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Preface 


Business Mathematics is a branch of applied mathematics that uses cal- 
culus, algebra and mathematical programming, to mention but a few. It finds 
applications in finance and economics. The mathematical analysis of stock 
markets, portfolios and the management of risks may also be covered by this 
subject. 

I began working on this book in October 2005 for the course Business 
Mathematics that I taught at Mahidol University at Kancanaburi. This vol- 
ume is the collection of hand-outs which were given to students before each 
lecture. Apart from examples and exercises, exam papers are also included. 
They are put in the appendix and represent in a way more examples. 

These hand-outs were used hand in hand with another set of printed 
pages whose contents are in larger letter, which were used on the camera 
projector. It would be good if in the future we could compile these hand-outs 
into lecture notes to be published as a book. 

There were 66 students in the class. From after the midterm exam my 
main concern has been for the students to do the exercises and problems by 
themselves. Since in practice this was hardly the case with the homework, we 
turned to holding within our class sessions of exercise practice and quiz. This 
has worked to a certain degree, all the time one main problem remain being 
the large size of the class, I do hope the students have learnt the importance 
of problem-oriented approach, and that each of them will find it useful in the 
future. I hope they have learnt from me as much as I have learnt for and from 
them. 


Kit Tyabandha, PhD 
Bangkok, 26%” April, 2006 
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Graph and derivative 
25" October 2005 


Definition 1. Let n be a positive integer. Then x” means that x is multi- 
plied to itself n times. 


8 


Definition 2. A monomial is an expression cosisting a real-numbered co- 
efficient times one or more variables each raised to the power of a positive 
integer. Adding and subtracting monomials to one another give us a polyno- 
mial. A univariate polynomial is a polynomial of one variable, of the form 
Anz” +...+a,x + ao, and its degree is the highest power of that variable, 
that is n. The degree of a polynomial is sometimes known as its order. A 
polynomial can be simplified into a product of two polynomials through the 
process called factoring. A polynomial equation is an equation of the form 
p(-) = 0, where p(-) is a polynomial. 

8 


Definition 3. A quadratic equation is a second-degree polynomial equation, 
that is one of the form ax? + bx +c = 0, where a, b and c are constants and 
a #0. 

8 


Example 1. The factors of mz? + nx + p are az + b and cx + d, where 
ac = m, bd = p and ad+ bc = n. Quadratic equations may be solved by 
factoring. 


Example 2. A quadratic polynomial of the form x?+bz can be transformed 
into the form (x+a)? by adding (b/2)? and then factor the result. Completing 
the square in the quadratic equation ax? + ba + c = 0 leads one to 


é b\? 2% —4ae 
x = 
2a 4a? 


which gives us the quadratic formula 


_ —b+ V0? = 4ac 


ie 2a 


Definition 4. At a point on the graph where the first derivative is zero, if 
the second derivative is positive, then that point is a local minimum; if the 
second derivative is negative, it is a local maximum; on the other hand, if the 
second derivative is either zero or is undefined, it is an inflection point. An 
asymptote is a straight line to which a non-linear curve smoothly approach as 
it goes towards infinity, never reaching it. At any point on a graph, if the first 
derivative is positive the function at that point is increasing, and if negative 
it is decreasing. At any point on a graph, if the second derivative is positive 
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it is said that the function is convex at that point, and if negative the latter 
is said to be concave. 


8 


global maximum 
a 


oe 


global minimum 


(a) (b) 


Figure 1 (a) global mazimum-, and (b) global minimum points of a 
function. 


y'=0 
y’< @) 


local maximum 


a 


local minimum 
y =0 
y’ 9 > @) 


Figure 2 Local minimum- and local maximum points of a function. 
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concave up 


curves down 


y’’> 0 


y’<0 


concave down 


curves up 


(a) (b) 
Figure 3 The two curvature types, namely (a) concave up (y” > 0) and, 
(b) concave down (y” <0). 


y’>0 
y’<0 y’>0 y>0 y<0 
y’<0 
(a) (b) 


Figure 4 The four possible curvatures in two dimensions, considering 


both y’ andy", (a) y” > 0 and, (b) y” <0. 

y’=0 
Wey, 
: () 


Figure 5 Inflection points, where y” = 0, (a) with y” increasing and, 
(b) with y” decreasing. 


y’=0 y =0 


y’<0 
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(a) (b) 


Figure 6 Stationary inflection points, where both y’ = 0 and y” = 0, (a) 
with y” increasing and, (b) with y” decreasing. 


Example 3. The cubic function of the form y = ax? + ba? + cx +d has 
local maximum and minimum, when these exist, at the point where the first 
derivative is zero, that is when 3ax? + 2b2 +c = 0. At such points, 


—-Or Jb = 3ac 
3a 


Depending on whether b? — 3ac is positive or negative, the local minimum and 
maximum either exist or do not exist. As examples, when a = 1, b= 2,c=3 
and d = 4 this value is —5, and therefore the value of « becomes complex. In 
the space of real number this means that there is neither a maximum- or a 
minimum point. Instead, in this case we have an inflection point, which is the 
point where the second derivative becomes zero, that is to say, 6ax + 2b = 0, 
or x =—. Figure 7 shows a plot of this case. 


y =ax® + ba? + ca +d, where a=1,b=2,c=3 andd=4 
10 T T T T T 


Bk 


6k 


~10 Ll Ll Ll Ll Ll 
3 


Figure 7 The cubic function y = x? + 2x7 + 3x44. 


Figure 8 shows the case where a= 1, b=5,c=4andd=3 
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y =aax® + bx? +cx +d, where a=1,b =5,c=4andd=3 
T T T T 


15 Ll L Ll Ll Ll Ll 
-5 -4 -3 -2 -1 0 1 


Figure 8 The cubic function y = x? + 5a? + 4x + 3. 


Example 4. The hyperbolic function of the form (a — a)(y — 6) = c has 
as its asymptotes the lines x = a and y = b. Figure 9 shows the graph when 
a=1,b=2andc=3. Here the asymptotes are x = 1 and y = 2. 


(x —a)(y —b) =c, where a= 1,b=2 andc=3 
10 T T T T T T T 


BE 


I 
I 
1 
\ 
\ 
6F i 
\ 
\ 
1 
\ 
\ 


| \ 

10 L L L L J 4 L 1 L 

10 -8 -6 -4 -2 0 2 4 6 8 10 
xc 


Figure 9 The hyperbolic function (x — 1)(y — 2) =3. 


Definition 5. The straight line has a form ax + by+c=0, or y=ma+d 
where m = —a/b and d = —c/b. The y-intercept is the point where x = 0, 
that is y = —c/b or y = d. The x-intercept is the point where y = 0, that is 
x = —c/a or © = —d/m. 


8 


Example 5. Given a sum of money, an isocost line represents the different 
combinations of two inputs of production that can be bought. The general 
formula is p.k + p,l = e, where k and | are capital and labour, px, and p; their 
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respective price, and e the allotted expenditures. In other words, p,k + pil — 
e = 0, or k = e/px — (pi/px)l. All the values concerned are positive, which 
means that the graph we are interested in is in the first quadrant only. 


Example 6. At the simplest case where there are only two activities to 
choose from, if x and y are the two activites, then trp, + typy =m, where tz 
and t, the number of hours spent on activities x and y respectively, pz and 
py are the price per unit hour of x and y, and m is the limit of income. 


Definition 6. The general demand function is of the form 


da = {(D,Y, Ps Pe; ta, 4, ---) 


where qq is the quantity demand of good z, p the price of x, y the income of 
the consumer, p, the price of substitute goods, p,. the price of complementary 
goods, t, the taste or fashion of the consumer, and a the advertisement level. 
In its simplest case where all other factors are constant, the demand equation 
takes the form p = c, — coqa, where p is the price-, while gq the quantity 
demanded of the good x, and c,; and cz are positive constants. The general 
supply function is of the form q, = f(p,c, po, te, 2, 0,-..), where qs is the quan- 
tity supplied of good 2, p the price of x, c the cost of production, po the price 
of other goods, t. the available technology, n the number of producers in the 
market, and o other factors, for example tax and subsidies. The simplified 
relation for the supply is p = c; + cogs, where q, is the quantity of x supplied, 
and c,; and cg are positive constants. 


8 


Example 7. When a tax of t per unit is imposed, the supply function 
becomes p — t = cy + Coq, where c; and c2 are positive constants, and the 
total cost function becomes c;, = ce +(k+t)q, where k is the cost of producing 
each unit. Here also cy = kq. 


Example 8. Revenue is the amount of money received when a firm sells its 
output. The relation is r; = pq, where r; is the total revenue, p the price and 
q the quantity. 


Definition 7. Elasticity « of x with regard to y means the ratio of the 
change in x to the change in y. Here y could be some economic variable, for 


example price or income. Hence, the price elasticity of demand is €g = meee 
the price elasticity of supply €, = cas a The point elasticity is the elasticity 


value calculated at a point, the arc elasticity or the midpoint elasticity is 
the same averaged over an interval. The arc price elasticity of demand or 


supply is then ég = eee The arc income elasticity of demand or supply 
: — dq yitye 
1S Ey = dy gitar’ 


8 


Elasticity measures the sensitivity or responsiveness of a certain quantity 
to changes in some variable. 
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Definition 8. A function of one independent variable is a relation in the 
form y = f(x) such that there exists one and only one value of y in the range 
of f for each real number x in the domain of f. The variable y is called the 
dependent variable. 


8 


Definition 9. A function is called multivalued function if the opposite to 
Definition 8 is true, that is there exist more than one values of y for some x. 


8 


Definition 10. An implicit function is a function in which both dependent- 
and independent variables appear on the same side. An explicit function is one 
where the dependent variable is on the left hand side-, and the independent 
variable on the right hand side of the equation. 


8 


Example 9. Examples of explicit functions are y = 5%, y = 2° + 2? —7 
and y = e* + (1 + 1)lInz. Examples of implicit functions are «+ y = 1, 
1 


x? +3cy-—yty?= ary and x? Iny = e¥(4#+ 2°). 


Definition 11. The exponential function has the general form y = a”, where 
the base a is aconstant and z is called the indez-, power-, or the exponent of the 
exponential function. The logarithmic function is then log, y = x, providing 
that both a and y are positive real numbers. We call a® the «> power of a, 
and call log, y the logarithm of y to the base a. When the logarithmic base 
is the natural number, we write log, n = Inn. 


8 


The logarithmic function is the inverse of the exponential function and 
vice versa. 


Theorem 1. Let p and q be real numbers, a and 6b positive numbers, and 
m and n positive integers. Then, 


a. ai= q?t4 
P 
a _ 
— =aqP-4 
qd 
(a?)4 — qP4 


a° = 1, provided that a 4 0 


ean 
aP 
(ab)? = aPbP 
a= at 
Vam =an 
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8 
Theorem 2. 
log, mn = log, m+ log, n 
log, ape log, m — log, n 
n 
log, m? = plog,m 
l 
log, n= cle 
log, a 
8 


Both Definition 12 and Definition 13 are a definition of limits in calculus. 
Definition 13 is a more formal and definitive one. 


Definition 12. If f(x) is a function which draws closer to a unique finite 
real number 7 for all values of x as the latter draws closer to a, but « 4 a, 
then / is called the limit of f(a) as x approaches a. In notation this is, 


lim f(x) =1 


na 


8 


Definition 13. For a function f(x), limp4 f(x) = 1 if and only if for every 
€ > 0, there exists 6 > 0 such that |f(a) — | < € whenever 0 < |x — al < 6. 


8 


By Definition 13 we mean that one can get f(a) to be as close to | as one 
wish, provided that x gets close enough to a. 


Theorem 3. Providing that both lim,-,, f(x) and limz_,, g(a) exist, and 
let k be a constant and n a positive integer, then the rules of limits are the 
following. 

lmk=k 


zLZ-a 
lim x2” = a” 
“za 


tim kf(x) =k jim f(a) 


lim [f(z) + g(x)] = lim f(x) + lim g(z) 


La La La 
Jim [f(e) - g(2)] = lim f(2) - lim g(c) 


A) = lims—ra ile) , provided that lim g(x) 4 0 
limz+a g(x) ra 
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lim [f(a)]” = [im f(x)| " , provided that n > 0 


«wa wa 


8 
Definition 14. Let y = f(x). Then, the derivative of y with respect to 2 is, 


dy See f(a + h) — f(a) 
dx = h-0 


The various notations for the derivative include df(«x)dz, 2 ‘(a), y’, Dy 


and D(f(a)). The process for obtaining this is called differentiation. 
8 


Theorem 4. Let y = 2”. Then the power rule for differentiation is the 
following. 

dy n—1 

— =nx 

dx 


8 


Theorem 5. The derivatives of exponential and logarithmic functions are 

respectively 92 = e* and 4@2 = 2. 

8 

It turns out that e” is the only function which is changed by neither dif- 

ferentiation nor integration. Leonhard Euler (1707-83) found that the natural 
number e = limn soo (1+ 4") = 2.71828... 


Exercise 1. Find the first derivative of the following: 


1. y = 2e” 2. y = 2In(z) 3. cy = M2 4+Ine 
4. p= 100(1 — e*) 5.q=VI+4 6. p=1+1lng 
7 q= "P -SR 8. ca =Ing+ 7 9. p= 98.55 


10.c=2?-—e*4+lnz 11.c= Mt +0.8leY 12. x = 34.3e° 


Theorem 6. Let y be a function of u, and u a function of x. Then the 
chain rule states that, 


dy _ dy du 
dx dudz 
8 
Exercise 2. Differentiate the following: 
1. p=5+ “y YS Tpoatt 
3. p=In(422) .y = (34 —5)7 


5. q = 0.95V/2I? + 31 c= V9 +4¢? 
7. y= (2—0.72)-8 . t = 100 + 30e9 84 
9. y=le th + sit —In(3z?—-274+5) 10. p= ai + 2In(q) 
11. s = 123(1 — 3e~° 4) 12. y= —+ 
(52+6) 3 


an » NY 
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Theorem 7. If y = u(x)v(x), then 


dy dw a du rn du 
=u u v 
dx "de de aa 
8 
Exercise 3. Find the first derivatives of the following: 
1. c=3qV2+¢4 2. p= Bare" 
3. y= a) 3/4a A. y= (#)e* +5 
5. s = 10ye—91Y 6. Ca = ji —qingt+ V4 
7 cr=yln(y®?)+V3yny 8. I= ; In(3y + 2) 
9. p= qe? In(q? — 8) 10. ¢ = & +aVq—1—q+ 100 
11. y=t?e' + 10lnt 12. g= olny 
Theorem 8. If y= nee then 
dy _ UGp — VGe 
dx v? 
8 


Exercise 4. Find the first derivative of the following and simplify your 
answer: 


— t-1 — 1-t a | 
1 Y= 2.¢t 14t 3.¢= "8 
— enz+1 = q = ate 
4.y=2mstl 5, p=100- 6. c= 500 (1- —*) 
3 26-4 f=I 
7 p= ia 8 C= Garang Oe Se 


Definition 15. The total cost c, comprises a fixed cost and a variable cost, 
that is ce = cf + cy. The total revenue r; is price times output, rz = pq. 


8 


Definition 16. The marginal cost cm is the change in total cost caused by 
the production of an additional unit. The marginal revenue rm, is the change 


in total revenue coming from the sale of an extra good. In other words, 
a dee _ dri 


Cm = Gq and rm, = aa where g is the output. 


8 


Definition 17. The average cost c, is the total cost per unit output, cg = ae 


The average revenue rq is the total revenue per unit output, rg = ae 


8 
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Example 10. From Definition 15, r, = pq, and from Definition 17, r, = 
raq, therefore p— rg. From Definition 17, cg = ai and from Definition 15, 
cH 


Ca = Cf +cy. Therefore cg = Caf + Cay, where the average fixed cost caf = a 


and the average variable cost cay = re 


Procedure 1 calculates and draws graphs of the total-, average-, and 
marginal costs, which are respectively c;, cg and c,,. Similar graphs of other 
quantities may be drawn likewise, for example ones for the total-, average-, 
and marginal products. 


Procedure 1 Total-, average-, and marginal graphs 


Given: c;(x) 
eae S 
Cm © C 
for each f € {c4,Ca,;Cm} do 
find the critical values x, for f’ = 0 
find f” 
n< |x,| 
for 1=1 to n do 
if f” (xo) > 0 then 
f(a£) is convex and is the relative minimum of f 
elseif f’(x¢) < 0 then 
f(a£) is concave and is the relative maximum of f 
endif 
endfor 
find inflection points x» from f” = 0 
endfor 


plot c:(x), then ca(x) and C(x) 
8 


Problem 1. Let the total cost function be c = q? — 40q? + 800. Sketch the 
graph to show the relationship between total-, average-, and marginal costs. 
If c = f(q), then c,, = dc;,/dq, and if r; = f(q), then r, = dr;/dq, where c; is 
the total cost, r; the total revenue, and gq the level of output. 


Problem 2. The demand function for a monopolist is q = 100 — 4p. Find 
Tt; Tm and rg, then evaluate these at q = 10 and explain the results. Further, 
find q when r, = 0, then find r,, when r, = 0, and then explain whether 
one should sell at this value of g. Draw graphs of r:, rm and rg on the same 
diagram. 


Problem 3. The demand function for a good is p = 200 — q!®. Find r;, 
Tm and rz, and then compare the slopes of the r,,- and rg curves. Evaluate 
Tt; Tm and ra at q= 11 and at gq = 26, and give an explanation of the values 
obtained. Find q when r,, = 0, and find gq when rz = 0. When does the sales 
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of more units begin to reduce the total revenue? Draw graphs of r;, r», and 
rq on the same diagram. 


Problem 4. The fixed costs of a firm are 2000 and the variable costs are 
4.7q. Write the equation for c and evaluate the ¢ at q = 35. Write the 
equation for c,, and evaluate the c,, at q = 35. Explain the meaning of c,, 
for this function. 


Problem 5. The average cost function of a firm is given by cg = q? —10q+ 
210 4 66. Write expressions for c; and c,,, and then calculate c; at q = 17. 


Write the equations for cr and cy. 


Problem 6. The average cost function of a firm is cg = OTe Differentiate 
Ca With respect to qg, then describe and explain how the former changes with 
the change in the latter. 


Definition 18. The relationship between input and output is called a pro- 
duction function, q = f(l,k,r,te, s,e,...), where | is labour, k phical capital 
such as buildings and machines, r raw materials, t. technology, s land, and 
e enterprise. Assuming a short period of time, then / becomes the only in- 
dependent variable and the other remaining factors are parameters, that is 
fixed, and therefore g = f(1). Then the marginal product of labour is pim = a4, 
and the average product of labour is pia = 4. 


8 


The value pjq in Definition 18 is a measure of productivity. 


Definition 19. The marginal propensity to consume is pem = ene The 
marginal propensity to save iS Pom = oe The average propensity to consume 
iS Dea = re The average propensity to save is Dsq = re Here y is the income, 
c the consumption, and s the saving. 


8 


Example 11. Since y = c+ 5s, it follows that peg + Psq = 1. Further, 


dy _ de ds —_ 
ay = ag + ag? and hence 1 = pem + Psm.- 


Example 12. The consumption function may be described by c = cg + by, 
where cg and 0 are positive constants. Then peg > Pem and Pgm > Dsa- 


Definition 20. The profit is 7 = r,— cq. At the break-even point x = 0, 
that is Te = Ct. 


§ 
Problem 7. Find the elasticity of demand when gq = re where a and c are 
constants. 
Problem 8. The demand function for train journeys is gq = a where q is 


the number of passengers in thousands. Find the elasticity of demand in this 
case. Determine the percentage change in demand when the fare increases or 
decreases by 10 per cent. 
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Problem 9. The demand for mineral water is given by g = 100 — 2p, where 
p is the price per bottle and q the number of bottles demanded. Find the 
expressions for rz, Tm and rg. Find the quantity and price that maximise the 
revenue. Find an equation for the price elasticity of demand firstly in terms 
of p, and then again in terms of g. Show that r; is maximum and r,, is zero 
when €g = —1, and find q at this eg. 


Problem 10. A demand for a wine is given by p = 1200e~°-°94, where 
p is the price per bottle and g the number of bottles demanded. Write the 
equations for rz, rm and rg. Find the price and quantity when the revenue is 
maximised. Show that ¢g = —1 when 7; is maximised. 


Problem 11. Derive the equation 


1 
rm =v(1+2) 


Problem 12. Show that the profit becomes maximised when 


Cm 
PFS 
1 
(+2) 


Problem 13. A shop selling shirts has a demand function p = 300 — 10q 
and a total cost function ¢ = 150+ 8q. Find the equations for r; and 7. 
Find the number of shirts which must be sold firstly in order to maximise 
the profit, and then again to maximise the total revenue. Show that ry, = re 
when the profit is maximised. 


Problem 14. The average cost function of a mobile phone is cg = 10+ 22%, 


and the average revenue function for the same is rg = 35. Find ry, c, Tm 
and c,,. How many mobile phones must be made and sold in order to break- 
even? Find the profit function, and show that neither profit nor revenue 
has a maximum. Explain using r,, and cm, why there is no maximum. Plot 
the graphs for rz, c, and 7 on one diagram, and plot r,, and c» on another 
diagram. Comment on the graphs. 
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Calculus of multivariable functions 
1° November 2005 


Definition 21 talks about functions of n independent variables, and Def- 
inition 22 the partial derivatives of these. In Theorem 9 we have the product 
rule, in Theorem 10 the quotient rule, and in Theorem 11 the generalised 
power function rule. 


Definition 21. A function y = f(x1,...,2,) is called a function of n inde- 
pendent variables if there exists one and only one value of y in the range of f 
for each tuple of real number (21,...,2n) in the domain of f. Here y is called 
the dependent variable while x;, 1 =1,...,n, the independent variables. 


8 


The word tuple in Definition 21 means an ordered list. It is also known 
as n-tuple, where n is the size of the list. 


Definition 22. Let a multivariable function be y = f(a1,...,@,). The 
partial derivative of y with respect to x;, where 1 <i <n, is a measure of the 
instantaneous rate of change of y with respect to 2; while x; is held constant 
for all 7 Ai, where 1 < 7 <n. This partial derivative is defined as 

Oy . £0..,0¢9 + Aa;,...) — f(a1,...,2n) 


= lim 
Ox; Axi 0 Ax; 


and can be written in either one of the following forms. 


Oy of 
oe ae fp (Bij: 24 Bri) Les OF Ya, 
8 
Theorem 9. Let z= g(z,y)-h(z,y). Then, 
Oz Oh Og 
=e. +h. 
dz © Ox Ox 
ae a dh, a 
Zz E 4.28 
Oy Oy Oy 
8 
Theorem 10. Let z= ee and h(x, y) #0. Then, 
dz __h- ge —g- oh 
Or h? 
pad a dh 
Oz Five om 8° By 
Oy h? 
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8 
Theorem 11. Let z = [g(x,y)]”. Then, 
Oz n—-1 Og 
dc Oe 
and 3 4 
ie = ng} . 88. 
Oy Oy 
8 
Exercise 5. Find the first-order partial derivatives of the following: 
1. c= 100(1 + e7}72) 2. q = 5l°-7hoS — 4] — 3k + 10 
8. q=In3x+4+ 2rIny 4.z2=0(1+yt+y") 
5. p = 150e°-74# 6.q=Inra+Iny 
Tz=e4a%?+2r42rytay 8 z2z=2x7y° 
9, gq S200 Ot 10. u = 8x3 y3 
11. z=ay+ 4% gS e 
8 


In Definition 23 we find the meaning of second-order partial derivatives. 
Theorem 12 is about critical points, and Procedure 2 is a procedure for de- 
termining critical points. 


Definition 23. Let z = f(x,y). Then, the second-order direct partial deriva- 
tives are 

DOA apa. eared 

Ox \ Ox Oy \ dy 


oO? o2 
5 and respectively fyy, (fy)y, oa 


These are also written 


fr, (far )ars 


The cross partial derivatives are 


2 (8) yg 2 (2 
Oy \ Ox Ox \ Oy 


These are also written as 


On oO? 
fy, (fe)y, Bia and respectively fyx, (fy)2, Tar 


Exercise 6. Find the second-order partial derivatives of the following: 
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1.z=23 +2y 
Pog Sey? 
5.q=Inax+lIny 
vA q= 99-77 f.0-1 


kit Tyabandha, PhD 


2. z= 19-2k9-3 — \(10 — 1 — 3k) 
4. c=100(1 + e714) 

6. gq = 519-449? + 50 — 71 — 5k 
& qgq=Ind5¢+alny 


9. p= 1109-87 
11. u=10zx*y4 


10.z=ayt+¢ 
12.z=a(1+y+y’) 
8 


When the second derivative is negative, the curve is concave towards the 
origin. 


Theorem 12. For a multivariable function z = f(x,y) to be a relative 
maximum at (a,b) necessarily f,,f, = 0, and fra, fyy <0 and fre: fyy > (fry)? 
at that point. For the same at the same to be a relative minimum, necessarily 
f,,fy = 0, and fra, fyy > 0 and fyx-fyy > (fry)? there. Moreover, an inflection 
point is a point (a,b) at which fr -fyy < (fry)?, and both f,. and fy, have 
the same sign. On the other hand, a saddle point is a point (a,b) at which 
foo «fyy < (fry)?, but fre and fy, are of different signs. 


8 


Procedure 2 Procedure for determining a critical point of a function 
with two independent variables 


Given z = f(x,y) and a point (a,b), at this point, 
if f, = 0 and f, =0 then 
(a,b) is a critical point 
if fix -fyy > (fey)? then 
if fr~ <0 and fy, <0 then 
(a,b) is a relative maximum of z 
elseif f,,, > 0 and f,, > 0 then 
(a,b) is a relative minimum of z 
else { 
endif 
elseif fix -fyy < (fey)? then 
if fr. -fyy > 0 then 
(a,b) is an inflection point 
elseif frz - fyy <0 then 
(a,b) is a saddle point 


else { 
endif 
else 
test inconclusive 
endif 
else 
(a,b) is no critical point 
endif 
8 
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Problem 15. There are two dead ends in Procedure 2. The first one (T) is 
the case where fxx-fyy > (fey)? and either (fer = 0, fyy = 0), (fre = 0, fyy < 0), 
(fax = 0,fyy > 0), (fea < 0,fyy = 0), (fea > 0, fyy = 9), (fea < 0, fy, > 0), 
or (fra > 0,fyy <0). The second one ({) is where fro - fy = 0. Find out 
what happen in these cases, and thus complete the missing lines of logic in 
Procedure 2. 


8 


Definition 24 gives the meaning of derivative and differential. Example 
13 looks at the derivative & differential of functions of one variable & two 
variables. 


Definition 24. By derivative du we mean an infinitesimally small change in 
y with respect to an infinitesimally small change in x. By differential dy and 
da we mean an infinitesimally small change in the values of y and respectively 
x. 


8 
Example 13. For a function of one variable y = f(z), the total derivative 
is 
dy 
dx 


and the differential of y is 


_(¥ 
dy = (+) dz 


For a function of two variables z = f(x,y) partial derivatives are, the first- 
order partial derivatives 


Oz a4 Oz 
eT aaidin 
Ox ro) 
and the second-order partial derivatives 
rz Orz o? wg Or z 
ee = 2a; = Zp, and —— = Zyz 
Ox? Oy? YW” OyOx A OxOy - 


The total differential of z is 


of Of 
dz= (=) da + (5) dy 


and for small changes which are not infinitesimal, dr becomes Az and the 
incremental change formula is 


Of Of 
Azx{[=—]A —|A 
=» (52) 4+ (5) % 
Definition 25 gives the general production function. In Example 14 we 
look at the Cobb-Douglas production function in more details. Theorem 13 


gives the law of diminishing returns to labour and the proof thereof, while 
similarly does Theorem 14 the law of diminishing returns to capital. 
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Definition 25. The general production function is q = f(l,k), where gq is 
output of the production, / labour and & capital. The Cobb-Douglas produc- 
tion function in its general form is 


q = al®kP (1) 


where a is a constant and0<a<1,0<8<1,l>0Oandk>0. 


8 


Example 14. With the Cobb-Douglas production function, the marginal 
product of labour is, 


0 - 
Pim = 1 = . nal? "he (2) 


and the marginal product of capital 


0 01 Be 
Pkm = {k = iE = aBl*ke (3) 


From this we see that pim > 0 and prem > 0. 


Theorem 13. From the Cobb-Douglas production function we have the 
law of diminishing returns to labour, which states that qy < 0. 


Proof. From Equation 1 in Definition 25, 


Pq 8 (Aq\ — OPim _ Cee 
maz aa) a 2 
Since 0 < a < 1, therefore qy < 0. q 


Theorem 14. Using the Cobb-Douglas production function, the law of 
diminishing returns to capital states that qz, < 0. 


Proof. From Equation 1 in Definition 25, 


_q _ O (aq — OPkm _ (9 184 
tke 9x2 ak \ Ok) Ok ke 
which, with 0 < @ <1, tells us that qzz < 0. q 


In Example 15 we see the changes in marginal product values. 
Example 15. Using the Cobb-Douglas production function, 


da =te=aopl "Ke 


Therefore, qx > 0 and gx: > 0. In other words, pj increases as capital input 
k increases, and respectively pz, increases as labour input / increases. 


Example 16 shows us the average functions of labour and capital, Ex- 
ample 17 the marginal functions of labour and capital, and Example 18 the 
comparison between marginal and average functions. 
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Example 16. For the Cobb-Douglas production function in Equation 1 the 
average product of labour is 


Pla = i Sal ke (4) 
and the average product of capital is 
Pha = 7 = al*hP? (5) 


Example 17. Again using the Cobb-Douglas production function of Equa- 
tion 1, the marginal product of labour is 


0 
Pim = Sp = aal*1k! (6) 


and the marginal product of capital is 


0 
Pkm = ok = apleke (7) 


Example 18. From the APL equation, Equation 4, and the MPL equation, 
Equation 6, and since 0 < a < 1, therefore py < Pig. Similarly from the APK 
equation, Equation 5, and the MPK equation, Equation 7, since 0 < 8 < 1, 
we have prem < Pra: 


In Example 19 one sees the conditions for using labour, Equation 8, and 
the conditions for using capital, Equation 9. 


Example 19. A producer likes to have a positive marginal function, which 
means that the productivity increases as the input increases. But the second 
derivative is negative, which means that this rate of increase slows down as 
time goes by. In practice, the conditions for using labour are, 


_ oq dpim _ 074q 
Pim = a > 0, qd 3p 0, and Pim < Pia (8) 


The conditions for using capital are similarly, 


_ oq dprm _ Oq 
OR ak OR 


Pkm < 0, and prem < Pka (9) 
Definition 26 and Theorem 15 deal respectively with production function 

graphs and slope of an isoquant. 

Definition 26. An isoquant is a graph in two dimensions, k = f(l), plotted 

to represent a production function g = f(l,k). The slope ae is called the 

marginal rate of technical substitution. The value of this slope at (lo, ko) is 

denoted by ar las 
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8 


Theorem 15. The slope of an isoquant is the ratio of the marginal products. 


Proof. The total differential of g = f(1,k) is 


_ (0 Oq 
aa= (22) ars (28) a 


Along any isoquant, dq = 0, therefore, 


0= (3) di + (3) dk (10) 


This directly yield, after some manipulation, 


dk 
dl qk 


Or, from Equation 10 together with Equation’s 6 and 7, it follows that, 


Definition 27. In the Cobb-Douglas production function equation, Equa- 
tion 1, let both inputs / and k change by the same proportion, and let \ be 
the constant of this proportionality. Then gz = a(Al)*(Ak)°, which leads to 
q2 = A°+8q,. When a+ 8 = 1, the case is described as constant returns to 
scale, when a+ 8 < 1 as decreasing returns to scale, and when a+ 6 > 1 as 
increasing returns to scale. 


8 


Definition 28. A homogeneous Cobb-Douglas production function of order 
r is, 

f(Al, Ak) = ATE(L, k) 
where r = (a, 8). 


8 


The utility function in its general form is given in Definition 29. Defi- 
nition 30 gives the Cobb-Douglas utility function, Definition 31 the marginal 
utility, and Definition 32 the meaning of indifference curves. 


Definition 29. A utility function expresses utility as a function of goods 
consumed. In its general form this is, 


u = f(x, y) 


where x and y are the quantities of goods X and respectively Y consumed. 


8 
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Definition 30. The Cobb-Douglas utility function is in its general form, 
u = ary? 


where a is a constant, and0O<a<1,0<6<1l,x>Oandy>0. 


8 
Definition 31. The marginal utility for a utility function with one variable, 
u = f(a), is se = Uz = Ugm. The marginal utility for a utility function with 
two variables, u = f(x,y), is ae = Uy = Ugm and a = Uy = Uym:- 
8 
Definition 32. The indifference curve is a graph y = f(a) drawn to repre- 
sent a utility function u = f(x,y). Its slope oy is called the marginal rate of 
substitution. Setting the total differential equal to zero, 
Ou Ou 
0=du= {| — ])dr+ {| — ]d 
(az) e+ (sp) 
we find 
dy Ug 
dz tly 
and 
dy _ _ Uxm 
dx Uym 
8 


Partial elasticities are describe in Definition 33, partial elasticities of 
demand, and Example’s 20 and 21, respectively the partial elasticity with 
respect to labour and the partial elasticity with respect to capital. 


Definition 33. Let a demand function be 


da = £(Da, Ys Pb) (11) 


where qq is the quantity demanded of good a, pg the price of a, y consumer’s 
income, and py, the price of another good b. Then, the price elasticity of 
demand is, 


Ae 04a Pa 
OPa Ia 
The income elasticity of demand is, 
0 
cy = Oe ¥ 
OY a 


And the cross-price elasticity of demand is, 


2, — 24a Po 
. Ope da 
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8 


Example 20. With the demand function as in Equation 11, the partial 
elasticity with respect to labour is, 


And from Equation’s 6 and 4, this leads to, 


Bey Pim 
ql — 
Pla 


For the Cobb-Douglas production function, Equation 1, then ég, = a. 


8 


Example 21. Again, with the demand function as in Equation 11, the 
partial elasticity with respect to capital is, 


Oq k 
Eqk => aka 
Then, from Equation’s 7 and 5, 
ei = Pkm 
i Pka 


For the Cobb-Douglas production function, Equation 1, we have €gx = £2. 


8 
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Exponential, logarithmic and nonlinear functions 
8*? November 2005 


Definition 34 gives a definition of a function. Example 22 gives some 
examples of this. In Definition 35 we talk about variables and parameters of 
a function. 


Definition 34. <A function is an operator or a procedure which accepts a 
permissible input and transforms it into a unique output. The input is some 
nonempty set. If a function is defined to be y = f(x), then x is the input 
vector, y the output, and f(-) the function itself. 


8 


Example 22. If f(-) is the function of dressing, then its input is possibly a 
person and its output a dressed person. If f(-) is the function of making up, 
then the input is perhaps a girl and the output a made-up girl. 


Definition 35. Let y = f(a,x) be a function, where a is a set of all its 
parameters, and x a set of all its variables. Then y is its dependent variable 
and 2;, for all i € x, are its independent variables. In other words, x; vary, 
y follows, and a; could assume any value within the range of its permissible 
ones, but its value must be constant. 


8 


Definition 36 and Example 23 address inverse function, and respectively 
operator and inverse operator. 


Definition 36. An inverse function is an expression of the independent 
variable in terms of the dependent variables. The inverse of the function 
f(-) is denoted by f~1(-). If f(-) is a function which admits one independent 
variable, namely x, then one could express it as, 


y =f(2) (12) 

Its inverse function is then, 
f'(y)=2 (13) 
8 


Example 23. Both the function and its inverse may be thought of as being 
an operator operating on an input to produce an output. The function, 


y = f(x) 


is understood diagrammatically as, 


Yor f(-) -— 2x 
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while its inverse function, 
«=f '(y) 


is displayed as a diagram as, 


y>| f1¢) |72% 


Theorem 16 is related to the domain and range of inverse functions. 
Example 24 gives some examples of inverse functions. 


Theorem 16. An inverse function must always be a one-to-one mapping. 


Proof. Let f(-) be a function. Then f(-) can be either one-to-one or many- 
to-one, and therefore f~1(-) could turn out to be either one-to-one or one-to- 
many. But since f~! is also a function, so for each of the values in its domain 
the corresponding value in its range must be unique. This means that in 
cases where f~! turns out to be one-to-many, some constraints must be put 
on its input in order to make the output one-to-one, which then makes all the 
outputs from f~1(-) one-to-one. q 


Example 24. Table 1 gives some of the functions and their corresponding 
inverse functions which are fundamental in mathematics. 


f(-) £~"() 


addition subtraction 
multiplication division 
power root 
exponential logarithm 


Table 1 Some of the functions and their corresponding inverse functions. 
Then, letting a be a constant, Table 1 becomes Table 2 
10 ee ae) 


ZT+a «e-a 


“a:-a a 
x Se 
a® log, x 


Table 2 The notational forms of functions and their inverses. 


in which division and logarithm are both undefined for a = 0. 
Definition 37 gives some of the basic building blocks of mathematics. 
Example 25 then shows how these are built one on top of another. 


Definition 37. The inverse of the addition, 
y=r+a 


is the subtraction, 
y-a=2 
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The inverse of the multiplication, 


y = ax 
is the division, 
yo 
[=f 
a 
The inverse of the power, 
y=u" 


is the root, 


The inverse of the exponential, 
y=a 
is the logarithm, 
loggy = # 
8 


Example 25. Figure 10 is a diagram which shows how addition makes 
multiplication makes power function. 


ia Pio 2 
t 
Chetek e 
ee ed c:a 
a 
ap 
ata2=2x 


Figure 10 The building blocks of mathematics, from addition to multi- 
plication to power function. 


In Figure 10 we start from considering a variable at the base. If in- 
stead of doing this we begin by considering addition of some constant a, then 
eventually a becomes a parameter in our more complicated functions. 


a---a 
SS a” 
av 
a 
a+...+@ 
Sa — a:n 
nm 
t 
at+a=2a 
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Figure 11 Starting from a constant to obtain in the end the exponential 
function. 


Our derivation in Figure 10 gives us x? when b is an integer, and similarly 
that in Figure 11 gives a” when z is an integer, but both the power- and the 
exponential functions can be extended to cover cases where the powers are 
noninteger, that is to say, when they are real or complex numbers. In these 
cases, however, the output may no longer be real. 

Next, we look at the exponential function, which is defined in Definition 
38. Example 26 discusses this further, and the exponential to the power of 
zero is looked at in Theorem 17. 


Definition 38. An exponential function is defined as y = a*, where a > 0 
anda #1. 
8 


Example 26. The domain of the exponential function y = a® is the set 
of all real numbers, while its range the set of all positive real numbers. The 
function is convex and increasing when a > 1, and convex and decreasing 
when 0<a< 1. At x =0, the value of the function is y = 1 for any a > 0. 


Theorem 17. For any a £0, 


lim a” = 1 
x—0 
8 
Problem 16. Try prove Theorem 17. 
8 


Theorem 18 gives some of the rules of exponential function and Example 
26 looks at some exponential function. 


Theorem 18. Three basic rules of the exponential function are, 


Proof. Write, say, a” as, 


Sa 
m 


and similarly for a”. Then all three equations above become obvious. q 
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Example 27. Figure 12 gives a graph of the exponential function when 
a> 1. Figure 13 gives a graph of the exponential function y = a” when 


0<a<l. 
y=2° 


0 Ll Ll Ll Ll Ll 
-3 -2 -1 0 4 2 3 


Figure 12 Example of the graph of the exponential function when a > 1. 
Here the graph is that of y = 2”. 


y = (0.5)* 
8 T T T 


0 
-3 -2 -1 0 1 2 3 


Figure 13 An example of graph of the exponential function y = a” when 
O0<a<1. Herea=0.5. 


x 


From Figure’s 12 and 13, one may see that the graph of y = a”, where 
0 <a <1, is the same as the graph of y = b-*, where 4 =6b> 1. This is 
obvious since by putting a = ; into y = a” one arrives at y = b~”, and if 
0<a<1,thend>1. 


One of the place we find use of an exponential function is the growth- 
and decay curves mentioned in Definition 39. Example 28 gives several basic 
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growth functions. 


Definition 39. Let a> 1. Then the graph of y = a® is called a growth 
curve, while that of y = a~* is called a decay curve. 


8 


Example 28. There are basically three laws of growth, namely unlimited, 
limited and logistic growth, all of which involve an exponential function. The 
model is for unlimited growth, 


for limited growth, 


and for logistic growth, 
y(t) 


where a, m and r are constants. 


ea +ae—rmt 


Interests, and future- and present values are discussed in Example’s 29 
and 30. 


Example 29. The value of a principal p compounded annually at an interest 
rate 2 for t years is, 
s=p(1+i)’ 


where 7 is expressed in decimal points. For compounding m times a year, 


then, 
; mt 
S=p (1 + ) 
m 


If the compounding is continuous, at 100 per cent interest for one year, then, 


1 m 
s=p lim (1+=) = pe 
m—- oo m 
where e is the Euler’s constant, e = 2.71828... 


Example 30. For multiple compounding, 


i mt 
pi +ie)t =p (1 4 “| 


m 


the effective annual rate of interest is, 
4 m 
e (1 Fe ) ij 
m 
The effective annual rate of interest for continuous compounding is, 
te =e —-1 
Definition 40 and Example 31 are about discounting. 
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Definition 40. Discounting is the process of finding the present value p of 
a future sum of money s. 


8 
Example 31. Discounting when under annual compounding is, 
s=p(1+i)’ 


when under multiple compounding, 


: —mt 
p=s(1+) 
m 


and when under continuous compounding, 


p= seve 


When discounting, the interest rate 2 is called the rate of discount. 


Example 32. A discrete growth s = p(1+i/m)™ can be converted to a 
continuous growth s = pe” thus, 


‘ mt 

Dp (1 + =) = per 
m 
j mt 

In (1+ =) =lne™ 
m 


r=min(1+ 4) 
m 


j mt ; 
s=p (1 = =) = pe™™ n(1+ mnt 


Example 33. Reversing the sign of x, that is replacing « by —x, has the 
effect of reflection of the original graph with respect to the y-axis. Reversing 
the sign of y, that is replacing y by —y, gives a reflection of the same with 
respect to the z-axis. The graphs of y = a** remain always above the x-axis, 
in other words the function y = a** maps —oco < « < co toy > 0. The two 
functions y = a” and y = a~” are the reflection of each other with respect to 
the y-axis. It can be easily seen that the functions y = —a** are the reflection 
with respect to the z-axis respectively of y = a=”. 


Therefore, 


Definition 41 introduces the logarithmic function, and Example 34 gives 
some elaboration regarding this. Some examples of natural logarithm are 
given in Example 35. Theorem 19 gives rules for logarithm. 
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Definition 41. The logarithmic function with base a is defined to be the 
inverse of the exponential function, and is written y = log, x, where a > 0 and 
a#1. The logarithmic function of base 10 is called the common logarithmic 
function, and one of base e, where e = limyn-+o0 (1 + 1)" is called the natural 
logarithmic function. By the notation y = log, « we mean that the logarithm 
base a of x is the power to which a must be raised to get «x. 


8 


Example 34. The domain of the logarithmic function y = log, x is the 
set of all positive real numbers, its range the set of all real numbers. The 
function is concave and increasing for a > 1, and is convex and decreasing for 
0<a< 1. Note also that log, x is the power which a must be raised to get 
x. 


Example 35. Note that el™* = a = Ine* where a > 0, ena g = ine 
where x > 0, and e!™(*) = f(x) = Ine where f(x) > 0. 


Theorem 19. Four basic rules for logarithm function are listed in the 
following. 
log, m + log, n = log, mn 


m 
log, m — log, n = log, = 


log, m* = zlog,m 


log, n 
] =— 
rea log,, b 
8 
Problem 17. Prove Theorem 19, the theorem for rules of logarithm. 
8 


Definition 42 sets out the meaning of the elasticity of substitution. Ex- 
ample 36 discuss the values of the elasticity of substitution. Definition 43 is 
on the constant elasticity of substitution production function. 


Definition 42. The elasticity of substitution o is defined as, 


where f is called the least-cost input ratio, and me the input-price ratio. 


8 


Example 36. The value o = 0 means there is no substitutability, that is the 
two inputs are complements of each other and both must be used together in a 
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fixed proportion. The value o = oo means that the two goods may substitute 
each other perfectly. Ultimately, 0 < a0 < oo. 


Definition 43. <A constant elasticity of substitution production function 
is a production function where, unlike the Cobb-Douglas function, has an 
elasticity of substitution whose value is constant but not necessarily 1. In its 
typical form, it is, 


BIH 


q=a(ak-?+(1—a)l-*) 

where a is called the efficiency parameter, a the distribution parameter, 3 the 

substitution parameter. Furthermore, 3 determines 0, anda >0,0<a< 1, 
and 6 > —1. 

8 


Example 37 discusses logarithmic transformation of nonlinear functions. 


Example 37. Some nonlinear functions can be converted to linear functions 
using logarithmic transformation, for example the Cobb-Douglas production 
function, 
q = ak*l? 
which becomes 
Ing=Ina+alnk+6lnl 


Other nonlinear functions can not be converted, for example the constant 
elasticity of substitution production function, 


ot 
q=a [ak-F +(1- a)l-F] e 
which becomes just another nonlinear function, 
1 
Ing=Ina- B In [ak +(1- a)l~F] 


others 


Example’s 38 and 39 give some examples of the use of nonlinear function 
in economics, namely the nonlinear total revenue and the nonlinear total cost. 


Example 38. Let the total revenue be r; = pq, and the demand function 
p =a-— bq, where q is the quantity sold. Then r; expressed as a function of 
q is nonlinear, for r; = (a — bq)q = aq — bq’. 


Example 39. A more realistic equation for the total cost instead of c = 
a+ bq is the nonlinear function ¢, = aq? — bq? +cq+d in which the production 
cost increases with quantity in at a decreasing rate (c, < 0) up to the inflection 
point at gq = 4, after which it increases at an increasing rate (c/ > 0). During 
the first stage the cost per unit decreases once the initial investment has been 
spent. During the second stage the cost per unit increases since more capital 
needs to be invested in order to allow more production capacity. 


Definition 44 talks about polynomial, and Example 40 about quadratic 
equation. 
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Definition 44. <A polynomial is an expression in the form ue age”. 
Here n is called the order of the polynomial. If n = 2 the polynomial is 
known as a quadratic polynomial, if n = 3 a cubic polynomial, ifn = 4a 
quartic, n = 5 a quintic and n = 6 a seztic. If we let p(a) be a polynomial, 
then a polynomial equation is the equation p(x) = 0. A polynomial function 
is a function of the form y = f(x) = p(x). 


8 
Example 40. The quadratic equation ax? + br + c = 0 has the solutions, 
fama C2 ES b2 aad 4 
pas ee (14) 
2a 


These solutions are called the roots of the quadratic equation. Equation 14 
is called the ‘minus-b formula’. The values of x obtained from the minus-b 
formula give the intersections of the graph of the quadratic function 

f(x) = p(x) = az”? + br +c 


on the z-axis. The value b? —4ac determines how the graph of f(x) lies relative 
to the x-axis, that is, 


> 0, there are two x-intersections 
b? —4ac< =0, the graph touches the x-axis at one point 
<0, the graph never touches the x-axis 


Furthermore, the graph reverses its direction with respect to the y-axis at 
the critical point where f’(a) = 0, that is when x = —<. Consequently the 


critical point is 
b b? i 
2a’ 4a © 


The graph of f(x) is symmetric with respect to the vertical line which passes 
through the turning point, that is to say, the line x = of The y-intercept is 
at the point (0,c). 


Example 41 is about another form of nonlinear function, the hyperbolic 
function. 


Example 41. A hyperbolic relation is an expression of the form, 


(px + q)(ry+s)=t 


From this we obtain, 


ee. 
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where p, q, r, s and t are constants, hence so are a f b=4e¢ ome bel 


pat 


economics we sometimes find hyperbolic functions of the form, 


4a 
~ bate 


y 


For example, a demand function of a good may be given by, 


m 
q+ra=— 
Pp 
which leads to, 
= m 
p= ga 


where p and q are respectively price and quantity demanded of a good, while 
m and a are constants. 


The graph of Equation 15 has the z-axis, that is the line y = 0, as its 
horizontal asymptote, and has the line z = —c/b as its vertical asymptote. If 
all the parameters are positive, then the curve in the first quadrant decreases 
with a decreasing rate. 


Bibliography 


Edward T Dowling. Introduction to mathematical economics. Schaum’s out- 
line series,2°¢ ed. 1992(1980) 


Teresa Bradley. Essential mathematics for economics and business. 2"¢ ed. 
2002 


God’s Ayudhya’s Defence 26" April, 2006 83 


Business Mathematics, notes and projections kit Tyabandha, PhD 


Matrix 
15*” November 2005 


Definition 45. Let A= {a;;}, B = {b;;} and C = {c,;} be three matrices. 
Then C = A+B is called the addition of the matrices A and B if ¢;; = aij +b;; 
for all 7 and 7. 


8 


Definition 46. Let A = (a;;) be an m x n matrix and B = (b,;) an n x p 
matrix. Then the product AB is an m x p matrix C = (cj) where, 


nm 
ci = s Aik del 
k=1 


where 1<i<mand1<l<p. 


8 


Definition 47. The expression obtained by eliminating the n variables 
X1,..-,% from n equations, 


441%, +... + A1nIyn = 0 
(16) 
An1X1 +... + adnntn = 0 


is called the determinant of this system of equations, Equation 16. The de- 
terminant of matrix A denoted by various different notations, for example 
det(A), |Al, Yo (£a1b2c3 rae ), D(a,b2c3 <i -), or |ayboc3 Pes |. 


8 


Example 42. For a linear system of three variables, Equation 16 can be 
written as, 

a,x + agy + a3z = 0 

byx oe boy Tr b3z =0 (17) 


cx + Coy + ¢3z = 0 


Eliminating x, y and z from Equation 17 gives us, 


a, b2c3 = a b3C9 + 0361 C9 ab, C3 t a2b3c, a3b2c1 =0 


Definition 48. <A minor Mj; of any matrix A is the determinant of a 
reduced matrix obtained by omitting the 7*® row and the j*® column of A. 


8 


Theorem 20. Determinant can be determined by, 


k 
| A| = S- aig Ciy 
i=l 
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where Cj; is called the cofactor of aj;. The cofactor Ci; can also be denoted 
as a’), and, i 

Cij = (-1) Mi 
where Mj; is a minor of A. 


8 


Problem 18. Prove Theorem 20,the theorem for finding the determinant 
of a matrix by Laplacian expansion. 


8 


Definition 49. Any pairwisely ordered pair in a permutation p is called a 
permutation inversion in p ifi > 7 and p; < pj. 


8 


Theorem 21. Determination of the determinant can also be determined 


by, 
|A| = Seay? II Qi r(i) 
i=1 


Tv 


where 7 is a permutation which ranges over all permutations of {1,...,n}, 
and I(7) is called the inversion number of 7. 


8 


Problem 19. Prove the theorem for the determination of determinant by 
permutation, Theorem 21. 


8 
Theorem 22. Let a be a constant and A an n x n matrix. Then, 
|aA| =a"|A| 
|=Al=(-1D)" | A] 
|AB|=|A||B| 
|Z] =|AA™*| =|A||A“*] =1 
1 
A= 
|A~*| 
8 
Problem 20. Prove Theorem 22, the theorem on properties of determinant. 
8 


Definition 50. A function in two or more variables is said to be multilinear 
if it is linear in each variable separately. 


8 


Theorem 23. Determinants of matrix are multilinear in rows and columns. 


8 
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Example 43. Consider an 3 x 3 matrix, 


a4 a5 a6 
a7 ag ag 


ay ag | 


What Theorem 23 says about multilinearity of determinants amounts to say- 
ing that, 


ay 0 0 0 ag 0 0 0 a3 
|A|=]a. a5 ag)/+]as a5 ag}+|as4 a5 ag 
a7 ag ag a7 ag ag a7 ag ag 
and 
a, ag a3 0 ag a3 0 ag a3 
|A|=| 0 as ag)/+]as a5 ag}+|] 0 as ag 
0 ag ag 0 ag ag az ag ag 


Problem 21. Prove the theorem on the multilinearity of determinants, 
Theorem 23. 
8 


Definition 51. A conformal mapping is a transformation that preserves 
local angle. The terms function, map and transformation are synonyms. 


8 


Definition 52. A similarity transformation is a conformal mapping the 
transformation matrix of which is, 


A’ = BAB"! 


Here A and A’ are similar matrices. 


8 


Theorem 24. Similarity transformation does not change the determinant. 
Proof. The proof for this is simply, 
1 


| BAB~'| =|B||A||B-*|=|B]| Al 
| B| 


= 4 


Example 44. 


|B-1AB — I | =|B-1AB— B-!)IB| 


=|B-\(A—\DB| 
=|B||A-AI||B| 
=|A—AlI| 
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Definition 53. Let A be a square, n x n matrix. Then the trace of A is, 


Tr(A) => 3 Agi 
i=1 


Definition 54. The transpose of a matrix A = {a;;} is A? = {a;i}. 
8 


Definition 55. The complex conjugate of a matrix A = {a;;} is A = {@;,}, 
whereG=p—qiifa=p+q. 
8 


Definition 56. Let ¢(n) or ¢(x) be a positive function, and let f(n) or 
f(x) be any function. Then f = O(¢) if |f| < Ad for some constant A and 
all values of n and x. Here O is called the big-O notation which denotes 
asymptoticity. The notation f = O(¢) is read, ‘f is of order ¢’. 


8 
Theorem 25. Some other properties of the determinant are, 
|A|=|A™| 
|A| =A] 
|I+¢€A| =1+Tr(A) + O(e?), for € small 
8 


Example 45. For a square matrix A, 

a. switching rows changes the sign of the determinant 

b. factoring out scalars from rows and columns leaves the value of the de- 
terminant unchanged 

c. adding rows and columns together leaves the determinant’s value un- 
changed 

d. multiplying a row by a constant c gives the same determinant multiplied 
by c 

e. if a row or a column is zero, then the determinant is zero 

f. if any two rows or columns are equal, then the determinant is zero 


Problem 22. Prove the properties of determinant given in Theorem 25. 


8 


Theorem 26. Some properties of matrix trace are, 
Tr(A) = Tr (A) 
Tr(A + B) = Tr(A) + Tr(B) 
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Tr(@A) = aTr(A) 


8 
Problem 23. Prove that, 
(AT) "= (A7t)" 
8 
Theorem 27. 
(AB)? = BT AT 
Proof. rie = ‘s 
(B A Ne aa ( ee (a ie 
= bKia5K 
= AjrbK: = (AB) 5; = (AB); 
q 


Definition 57. Let A be a square matrix. Then the inverse of A, if it 
exists, is 47! such that, 
AA a 


Furthermore, A is said to be nonsingular or invertible if its inverse exists, 
otherwise it is said to be singular. 


8 


Example 46. For a 2 x 2 matrix, 


the inverse of A is, 


1 d —b 
Atl= 
ad — be & a 


If A is a3 x 3 matrix, then the inverse of A is, 


At = a faet (ms) 


where m,; is a minor of A. 


If A is an n x n matrix, then A~! can be found by numerical methods, for 
example Gauss-Jordan elimination, Gaussian elimination, and LU decompo- 
sition. 
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Example 47. The Gaussian elimination procedure solves the matrix equa- 
tion Ax = b by first forming an augmented matrix equation [Ab] and then 
transform this into an upper triangular matrix [{ai,} b’], where aj, are all 
zero except when i < j. Then, 


k 
1 / 
oe b; y Ajj 
a jit 


The Gauss-Jordan elimination procedure finds matrix inverse by first forming 
a matrix [A J], and then use the Gaussian elimination to transform this matrix 
into [I B]. The result matrix B is in fact A7'. 

The LU decomposition forms from the matrix A a product LU of two matrices, 
one lower- while the other upper triangular. This gives us three types of 
equation to solve, namely when 7 < 7,7 = j andi > j, where 7 and j are the 
indices of row and respectively column of the matrix product. Then, 


Ax = (LU)x = L(Ux) =b 


Letting y = Ux we have Ly = b, and therefore, 


by 
Y= 7 
lia 
i i-l 
Y= 7 b; > lig Yj 
Lit aq 
j= 
where 7 = 2,...,n. Then since Ux =y, 
Yn 
t= 
Unn 
1 n 
y= —|y%- S- UijX; 
N44 re 
g=itl 


wherei=n-—1,...,1. 


Theorem 28. Let A and B be two square matrices of equal size. Then, 
(AB) = BoA 
Proof. Let C = AB. Then B = A~!C and A = CB“, therefore, 
C= AB=(CB™")(A'C) =CBAC 
Hence CB~'A~! =I, and thus B~'A~! = (AB)"1. q 
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Definition 58. The Einstein’s summation is the simplification of notation 
by omitting a summation sign, keeping in mind that repeated indices are 
implicitly summed over, for example 5°, ai,aij becomes ajpaij, and >>; aja; 
becomes a;a;. 


8 


Definition 59. The multiplication of two matrices A = {a,;} and B = {b;;} 
is the matrix C = AB such that cizx = aijbjr. 


8 
Theorem 29. The matrix multiplication is associative. 
Proof. 
[(ab)c],; = (ab) ances = (airbix) Cry = Ait (ines) = aa(be)iy = [a(be)],,; 
q 


Example 48. From Theorem 29, which shows us the associativity of matrix 
multiplication, we could write the multiplication of three matrices as [abc];;, 
which is the same as writing a,.bj,c,;. And this applies in a similar manner 
to the multiplication of four or more matrices. 


Theorem 30. If A and B are two square and diagonal matrices, then 
AB = BA. But in general matrix multiplication is not commutative. 


8 


Problem 24. Prove Theorem 30, which is a theorem about non-commutati- 
vity of matrix multiplication. 


8 


Definition 60. A block matrix is a matrix which is is made up of small 
matrices put together, for example, 


[é >] 


where A, B, C and D are matrices. 


8 


Theorem 31. Block matrices may be multiplied together in the usual 
manner, for example, 


Aj By Ag Bo 2. A, Ao le BiC, A, Bo as By, Do 
Ci D, C2 Do a COC Ag Sane Di Cz C7 Bo a D, D2 


provided that all the products involved are possible. 


Problem 25. Prove Theorem 31. 
§ 
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Definition 61. Let A = {a;;} be an n x n matrix. Then A is called a 
diagonal matrix if aj4j = 0 when i 4 j. Here 1 < i,7 <n. In other words, 
a diagonal matrix has its components in the form aj; = cidj;, where c; is a 
constant and 6;; is the Kronecker delta, 
ae 
li=j 
8 
Theorem 32. A square matrix A can be diagonalised by the transformation 
A = PDP™", where P is made up of the eigenvectors of A and D is the 
diagonal matrix desired. 
8 
Problem 26. Prove Theorem 32, the theorem on matrix diagonalisation. 


8 


Example 49. Matrix diagonalisation can greatly help reducing the number 
of parameters in a system of equations. For instance, the systems Ax = y 
when diagonalised becomes PDP~!x = y, that is Dx’ = y’, where x’ = P~!x 
and y’ = P~'y. In this case, if A is an n x n matrix, we say that our new 
system obtained through the process of diagonalisation has canonicalised from 
n xX nto n parameters. 


Definition 62. A symmetric matrix is a square matrix A which satisfies 
Ay A. 
8 


Example 50. If A is asymmetric matrix, then A~!A? = J. 


Definition 63. Let A be a square matrix. Then A is said to be orthogonal 
if AAT = IJ. 


8 
Example 51. Definition 63 is the same as saying that A~! = A’. 
Theorem 33. A matrix A is symmetric if it can be expressed as A = 
QDQ’, where Q is an orthogonal matrix and D is a diagonal matrix. 

8 
Problem 27. Prove Theorem 33, the problem on symmetric matrix. 

8 


Example 52. Any square matrix A may be decomposed into two terms 
added together, that is A, + A, where A, is a symmetric matrix and A, an 
antisymmetric matrix, called respectively a symmetric part and an antisym- 
metric part of A. Furthermore, 


As = =(A+ A’) 


Nile 


and 


Aa = = (A-A®) 


NlRe 
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Linear algebra 
224 November 2005 


Definition 64. Let A be a square, nonsingular matrix. Then the inverse 
matrix A~! of A is a unique matrix for which, 


AA =F SATA 


8 
Example 53. An inverse matrix may be found using the formula, 
Avs ae Adj A 
| A| 


Example 54. Matrix equations of the form Ax = b can be solved with the 
help of the inverse matrix A~! as x = A7'b, where A is an n x n matrix, x 
a vector of size n whose components are variables, and b a vector of size n 
containing constants. 


Theorem 34. Let A be the coefficient matrix and A; a matrix formed 
from A by replacing the column of coefficients of x; with the column vector 
of constants. Cramer’s rule solves a system of linear equations through the 
use of determinants as follows. 


oe, — Mil 
| A| 
8 
Problem 28. Prove Cramer’s rule, Theorem 34. 
8 


Definition 65. Let a system of n functions not necessarily linear be 


YW =fi (Giga ag ln) 


Uni = fg (Lip eign) 


Then a Jacobian determinant comprises all the first-order partial derivatives 
of the system arranged in ordered sequence, that is 


OULD Yaw, ue 
Ox, OL: 
| uiwdm |_| 
| J | ~~ | 0a4,...,02n | _ % 2 . 
Oyn .,,  OYn 
O24 Oxn 


8 


42 26" April, 2006 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Theorem 35. Let a system of n equations be y; = fi (%1,...,%n), t = 
1,...,n. If | J| = 0, then y; are functionally dependent. On the other hand 
if | J| #0, then y; are functionally independent. 


8 


Problem 29. Prove Theorem 35. 
§ 


Definition 66. A determinant | H | composed of all the second-order partial 
derivatives, with the direct partials on the principal diagonal and the cross 
partials off the same, is called a Hessian. In other words, let a multivariable 
function be z = f(x,y). Then the Hessian of z is 


Zz z 
zyx Zyy 

where Zz = Zyz- Moreover, the first principal minor is | Hi | = zzz and the 
second principal minor is 

Zz z 

| H2,=|] 0°" O° | = ZanZyy — (Zey) 
z z 
ay yy 


8 


Theorem 36. Let a multivariable function be z = f(x,y), and let the 
first-order conditions z, = z, = 0 are met. Then a sufficient condition for z 
to be at optimum is ZerZyy > (Zey)? together with zz2, Zyy < 0 in case of a 
maximum and Zz, Zyy > 0 in case of a minimum. 


8 


Problem 30. Prove Theorem 36, the theorem on the optimality of a mul- 
tivariable function. 


8 


Definition 67. From Definition 66, if | H,| > 0 and | H2| > 0 the Hessian 
| H| is said to be positive definite, and the second-order conditions for the 
minimum are met. If | H,| < 0 and | H2| > 0 it is said to be negative definite, 
and the second-order conditions for the maximum are met. 


8 


Algorithm 1 Procedure to test for the optimality of multivariable func- 
tions of two variables. 


a f(z, y) 
find z, and zy, 
if z, —0 and z, = 0 then 
find 222, Zry and Zyy 
find H, and H» 
if | H,| >0 and | H2| > 0 then 
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| H| is positive definite 
elseif | H, | < 0 and | Hz > 0| then 
| H| is negative definite 
endif 
endif 


Definition 68. Let y = f(a1,...,%,) be function of n variables. Then the 
n*>_order Hessian for this function is 


Yr oct° Yin 
|| = : 
Ynl1 oa Ynn 


Then the first principal minor | H | is simply x11, and the i*” principal minor 
is 
Yur Yui 
| Hil = 
Yil oct Vii 
8 


Theorem 37. Let y = f(#1,...,@n) be function of n variables. Let the 
Hessian of y be represented by | H |. Then if all the principal minors of | H | 
are positive, then | H'| is positive definite and the second-order conditions for 
a relative minimum are met. If the sign of the principal minors alternates 
between negagive and positive, then | H | is negative definite and the second- 
order conditions for a relative maximum are met. 


§ 
Example 55. For y = f(#1, 2,23) the third-order Hessian is 
Yl Yl2— 13 
| | = |¥21 Y22 Y23 
Y¥31 Y32 = Y33 
where 
O7y O7y 


and so on 


a Ox?’ = 0x2021" 


The first-, second- and third-order Hessian’s are respectively 


yi Ye Yur Yi2) 13 

|Ail=yu, |A2|= and | H3|=|Yyo1 Yoo Y23 
Yai -Y22 

¥31 «3233 


If |H,| > 0, | H2| > 0 and | H3| > 0, then A is positive definite and the 
second-order condition for minimum is fulfilled. If | Hi | < 0, | H2| > 0 and 
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| H3| < 0, then | H| is negative definite and the second-order condition for 
maximum is satisfied. 


Definition 69. A discriminant is a determinant of a quadratic form. Let 
the quadratic form be z = ax? + bry + cy”, which is in matrix form 


eft JE 


Then the discriminant is 


b 
a 2 
|D|= b : 
2 
The first principal minor of the discriminant is |D,| = a, and the second 
principal minor 
b 2 
a> b 
D2| = 2) =ac-— — 
| Da | g c A 


Theorem 38. Let a quadratic form be 
z= ax? + bry 4+ cy” 


and let the discriminant of z be | D|. If | Di | > 0 and | D2| > 0, then | D| is 
positive definite and z > 0 for all x,y 4 0. If | Di | < 0 and | D2| > 0, then 
| D| is negative definite and z < 0 for all x,y £0. 


8 


Theorem 39. Let f(x,y) be a function subject to a constraint g(x,y) =k, 
where & is a constant. Then the optimisation of f can be done by first 
transforming f together with g into a new function 


Ce ay, X) oF f(2,y) +A(k—g(2,y)) 


and then solve the following equations, 


F,(x,y,A) =0 
F,(z,y, A) =0 
F(a, y,A) =0 


to obtain the critical values x9, yo and Ap at which F' and hence f are opti- 
mised. 

8 
Problem 31. Prove Theorem 38, the definiteness of a function by the 
discriminant. Prove the theorem for constrained optimisation with Lagrange 
multipliers, Theorem 39 above. 


8 
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Definition 70. In the constrained optimisation with Lagrange multipliers 
in Theorem 39 above, f is called an objective or origin function and F the 
Lagrangian function. 


8 


Definition 71. Let f(a1,...,2,) be a function of n variables subject to 
constraints g(x1,...,2n,). Let 


F(a1,...,%n,A) = f(a,...,¢n) tA(k- g(a1,...,2n)) 


Then the bordered Hessian | H | is defined as either 


Fy, Fie Pin 1 
Fy, 92 
|H|=] : 7 : 
Pui Fran gn 
gi g2 +t gn O 
or 
0 1 gn 
2 n Fu Fin 
Ea ie 
This is simply the Hessian 
Pit! 434 ae 
Fri aes Fan 


bordered by the first derivatives of the constraint with zero on the principal 
diagonal. The order of a bordered principal minor being determined by the 
order of the principal minor being bordered, | H| = | H,, | since in this case 
an n X n principal minor is being bordered. 


8 
Theorem 40. Let f(z1,...,@,) be a function of n variables subject to con- 
straints g(@1,...,%n). Let | H| be the bordered Hessian defined in Definition 
71. Then if | H2|,...,|Hn| < 0, then the bordered Hessian | H | is positive 


definite, and therefore is a sufficient condition for a minimum. If | H2| >0 
| H3| <0, | H4| > 0, and so alternatingly on, then | H| is negative definite, 
which is a sufficient condition for a maximum. 


8 


Problem 32. Prove Theorem 40. 
§ 
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becomes 


F(2,y, r) = f(a,y) +A(k— g(x, y)) 


The first-order conditions for optimisation are F,, = 
second-order conditions for optimisation can be expressed together as a bor- 


dered Hessian 


Business Mathematics, notes and projections 


Let f(z,y) be a function to be optimised subject to a con- 
straint g(x,y) = k, where & is a constant. Then the Lagrangian function 


, = Fy = 0. The 


- Pex Poy Gx 
|H|=|Fye Fyy gy 
Ix GIy 0 
or 
7 0 Ge Dy 
| H | =19e Free Pry 
Jy Pye Fyy 
Note 1. Theorem 39 gives the first-order conditions for optimising a func- 


tion subject to some constraints. Theorem 40 gives the second-order condi- 
tions for optimising a function subject to some constraints. 


8 


Definition 72. A Marshallian demand function gives an expression of the 
amount of a good that a consumer will buy as a function of commodity prices 
and income available. It is derived by maximising the utility subjected to a 
budgetary constraint. 


8 


Example 57. Let a utility be wu = qq2 which is subject to a constraint 
Pigdi + p2g2 = b, where b is the amount of income available, that is to say, our 
budget. Then the Lagrangian function is 


U = qigg + A(b— pig — p24q2) 


The first partial derivatives are then 


U1 = q2—- Api = 0 (18) 
U2 = qd — Ape = 0 (19) 
uy =b—pigi — prq2 = 0 (20) 


where wu 1, uz are respectively u,, and ug,. Simultaneously solving Equation’s 


18, 19 and 20 leads us to 
q2 


Pl 


j= 
p2 


Hence g2 = qipi/p2 and qi = q2p2/pi and from Equation 20 we have, 


Pg p2q2 
b= pig + po— = pogo + pi1— 
P2 P1 
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which yield us, for q; and qz2, the Marshallian demand functions which max- 
imise satisfaction of the consumer subject to income and prices. 


Next, we test the second-order conditions by firstly finding u 41 = 0, w22 = 0, 
U12 = U21 = 1, gi = pi and go = pe, which give us 


7 0 1 P1 
|H|=)1 0 pe 
Pi poz O 


which gives | Hj | = 2p;p2 > 0 Hence | H| is negative definite and thus u is 
maximised. 


Definition 73. The production process of producing one good usually 
requires the input of many other intermediate goods. Let x; be the total 
demand for product 7, and let b be the final demand for the product from the 
ultimate users. Then, 


LH GILT eee Oi ey + b; 


fori =1,...,n, where aj; is a technical coefficient which represents the value 
of input 7 required to produce one monetary unit’s worth of product 7. If we 
consider the total demand for every one of the products, then 


x= Ax+b 


where 
Ly 411 "°° Gin by 


In an1 tt Ann bn 


It follows from this that x = (I — A)~'b. The matrix A is known as the 
matriz of technical coefficients. It is also known as the input-output table, the 
rows being the inputs and the columns the outputs.. The matrix I — A is 
known as the Leontief matric. 


8 


Example 58. Ina complete input-output table, labour and capital would 
also be included as inputs. These give the value added by the firm. They 
are normally put as an extra row at the bottom of the matrix of technical 
coefficients A. The vertical summation of each column of the table is then 
equal to 1. 


Definition 74. Let A be a square matrix. Then a scalar \ such that the 
equation 
Av = Av (21) 
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holds for some vector v ¥ 0 is called an eigenvalue t of A, and the vector v is 
called an eigenvector of A corresponding to the eigenvalue A. The eigenvalue 
is also known as the characteristic root, or the latent root, while the eigenvector 
is also known as the characteristic vector, or the latent vector. 


8 
Note 2. From Equation 21 it follows directly that 


(A—AD)v =0 (22) 


Then A—ALI is called the characteristic matrix of A. Since v assumes a unique 
value and by assumption v ¥ 0, it follows that A— AJ must be singular, which 
means that its rows must be a multiple of one another. Now A — XI is zero 
if and only if the characteristic determinant | A — AI| of A is zero. In other 
words 


| A—AIT| =0 (23) 


which is called the characteristic equation of A. With Equation 23 there will 
be an infinite solution for v in Equation 22. In particular, if v is a solution, 
that is if it is an eigenvector, so is kv for any k # 0. We force a unique 
solution by using the normalisation 


ye 


Then the sign-definiteness of A can be determined from the characteristic 
roots \’s. Thus if all \’s are positive, then A is positive definite; and if 
negative, negative definite. Let at least one A be zero, which is neither positive 
nor negative, if all the remaining ’s are nonnegative, then A is positive 
semidefinite; and if they are nonpositive, negative semidefinite. Lastly, if 
some of the ’s are positive while others are negative, then A is indefinite. 


8 
Problem 33. Prove all the necessary details need to be proved in Note 2. 
8 
Note 3. We have seen in Note 2 how, having found \;, where i = 1,...,n, 


we find through normalisation the corresponding, unique v;. On the other 
hand if we have found first the v,’s, their corresponding \;’s may be found 
by first forming a transformation matria 


T =[v1 ..- Val 


+ The word eigenvalue is a half-translation of the German word Eigen- 
wert. The latter means ‘appropriate value’ since Wert means ‘value’ and 
eigen means ‘proper or appropriate’. 
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and then the corresponding eigenvalues or the characteristic roots are ob- 
tained from 


Ay oO 2s: oi 
mare) 
0 Poker 
8 
Problem 34. Prove all the necessary details in Note 3. 
8 


Definition 75. The vector equation, Equation 21, has as its solutions the 
zero vector v = 0 together with all the corresponding eigenvalue-eigenvector 
pairs. The set of all the eigenvalues of A is called the spectrum of A. The spec- 
tral radius of A is then the largest of all the absolute values of the eigenvalues 
of A, that is to say, 

max |Az| 


The set of all eigenvectors v;;, together with 0, forms a vector space called 
the eigenspace of A corresponding to j;. 


8 
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Examples for linear algebra 
26" April, 2006 


1. Solve using the Gaussian elimination, 


321 pa 8x2 = 53 
621 + 2x2 = 50 


Solution. Write the equations in an augmented matrix, 


[A B]=|5 8 aa 


6 2 50 


Row 1| times 3, 
1 8 33 
3 3 
E 2 | 
Row 2 subtracted by 6 times of Row 1, 
14 4 
0 -14 —56 


Row 2 times -i. 


Row 1 subtracted by § times of Row 2, 


1 0 7 
0414 
Therefore %; = 7 and X2 = 4. 
# 
2. Let 
1 -3 3 
A= {3 —-5 3 
6 -6 4 


Find all the eigenvalues and a basis of each eigenspace. Can A be diago- 
nalised? Why? 


Solution. Form the characteristic matrix tf — A and find its determinant to 
obtain the characteristic polynomial A(t) of A; 


t-1 3 -8 
A(t)=|tl-A|=|] -3 t+5 —-3 | =(¢+2)?(t-4) 
-6 6 t-4 
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The roots of A(t) are —2 and 4, hence the eigenvalues of A are —2 and 4. 


# 
Next, find a basis of the eigenspace of the eigenvalue —2. Substitute 
t = —2 into the characteristic matrix tf — A, thus obtain the homogeneous 
system 
-3 3 -3 x 0 
-3 3 -3 y | =| 0 
-—6 6 -6 z 0 


In other words, 
—32 + 3y —3z=0 
—32 + 3y —3z =0 
—6z + by —6z =0 


that is x—y+z=0. The system has two independent solutions, for example 
ex=1,y=1,2=O0andz2=1,y=0, z=-1. Therefore u = (1,1,0) and 
v = (1,0, —1) are independent eigenvectors which generate the eigenspace of 
the eigenvalue -2. In other words, u and v form a basis of the eigenspace of 
-2. This means that every eigenvector belonging to -2 is a linear combination 
of u and v. 


Similarly, find a basis of the eigenspace of the eigenvalue 4. Substitute 
t = 4 into the characteristic matrix tf — A to obtain the homogeneous system 


3. ° 3) aa ae 0 
-3 9 -3|[y]=[0 
6 6 0 z 0 


This is the same as writing 


3x + 3y — 3z =0 
—32 + 9y —3z=0 


—6x + by = 0 
which can be reduced to 
xr+y—z=0 
2y—z=0 


Since the number of independent equations is less than the number of variables 
by one, this system has only one free variable. Therefore any particular non- 
zero solution, for example x = 1, y = 1, z = 2 generates its solution space. 
Hence w = (1, 1,2) is an eigenvector which generates, and so form a basis of 
the eigenspace of the eigenvalue 4. 


Since A has three linearly independent eigenvectors, A is diagonalisable. 


# 
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Let P be the matrix the columns of which are the three independent 
eigenvectors, 


ie. te» 2 
P={1 0 1 
0 2 


-1 
Then 
—2 0 0 
P-'AP=| 0 -2 0 
0 oO 4 


The diagonal elements of P~'AP are the eigenvalues of A corresponding to 
the columns of P. 
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Exercises for linear algebra 


26" April, 2006 


3. Find values of the variables x, y and z for each of the following systems of 
equations, using Gaussian elimination. 


2e—-2y-z = 


3 4a+3y+z = 5 

(i) x-ytzez = 2 (it) 2ea-y-z = 4 
rty+2z = 38 Py sg. “Se 3 
2e+2y—7z = 10 3a+2y+z2 = 16 

(iit) —“-y =" 5 (iv) ety+z2 = 5 
3a+2y+2 = -1 22 y -?Z SF 
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Linear programming 
29" November 2005 


Definition 76. <A problem of optimisation is one in which one tries to 
maximise or minimise a certain quantity called the objective, which depends 
on a finite number of variables. These may be either independent or related 
to one another through some constraints. 


8 


Definition 77. A mathematical programme is an optimisation problem in 
which the objective and the constraints are given as functions or mathematical 
relationship. In other words, 


optimise: z= f(x1,...,2n 


subject to: g;(@1,.--,2n) 


Some constraints are explicitly stated as requirements, others are hidden con- 
ditions. These latter need to be pin-pointed through the study and under- 
standing of the model and its inputs. 


8 


Definition 78. <A linear programme is a mathematical programme all the 
functions involved of which are linear. This means that, 


f(@1,---,0n) =Car +++ + enon 
gi(@1, ad sa) = 41%, +++ + GinIn 
where i = 1,...,m and c; and aj;, 7 = 1,...,n, are constants. If there is 


an additional restriction on the input variables that they be all integers, then 
the optimisation problem is called an integer programme. A mathematical 
programme which is not a linear programme is said to be nonlinear. 


8 


Definition 79. <A quadratic programme is a mathematical programme in 
which all the constraints are linear and the objective function is in quadratic 
form, which is in general, 


n n n 
f(x, cee ,Ln) = 2 > Cig Lily + S- dj Xj 
i=1 


i=1 j=1 


where c,; and d; are constants. 
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Definition 80. A linear programme is said to be in standard form if all 
the constraints are equalities and if one feasible solution is known. In other 
words, our problem is now 
optimise: z=c’x 
subject to: Ax=b 
with: x>0 


8 


Definition 81. One may change any linear programme into the standard 
form by adding a slack variable to the left-hand side of a constraint of the 
form >> a;;x; < b; to obtain 


n 
Yat; + Lp, = bi 
j=l 


where py > n and k = 1,2,.... Similarly one may add a surplus variable 
to the right-hand side of a constraint of the form )Ja;;7; > 6; to obtain 


D1 igtj = bi + tq, 
n 
DF digg — Bq, = 4 
j=l 


where q > n and / = 1,2,.... Next, all the slack and surplus variables are 
added to the objective function with zero coefficients. Then if we add an 
artificial variable to the left-hand side of each constraint where there is no 
slack variable, then the initial feasible solution is x) = b, where x is the 
vector of slack and artificial variables. The artificial variables are added to 
the objective function with a large negative coefficient —M. 


8 

Definition 82. A set of n vectors of m dimensions {vj,...,Vv,} is said to 

be linearly dependent if there exist some constants a1,...,Q@ , not all of which 
are zero, such that 

Q1V, +-:-+Q,Vpn = 0 (24) 


It is said to be linearly independent if the condition in Equation 24 implies 
QQ, = a9 =::: =a, = 0. 


8 


Theorem 41. Consider a set of n vectors of m dimensions. If n > m, then 
the set is linearly dependent. 


8 


Problem 35. Prove Theorem 41. 
§ 
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Definition 83. A vector v is called a convex combination of vectors v1,..., 
v,, if there exist some nonnegative constants 61,..., Gn, where 


Bits +Bu=1 


such that 
v= Bivit-+++BnVn 


8 


Definition 84. A set of m-dimensional vectors is said to be convez if for any 
two vectors belonging to the set the line segment between them also belongs 
to the set. 

8 


Theorem 42. All points on the line segment joining any two vectors may 
be expressed as a convex combination of the two vectors. 


8 


Problem 36. Prove Theorem 42. 
§ 


Definition 85. A vector v is called an extreme point of a convex set if it 
can not be expressed as a convex combination of two other vectors in the set. 


In other words, an extreme point of a convex set K is a point x in K that 
cannot be written as x = dy + (1— 0)z with 0 < 6 < 1, y and z in K, and 
y # z, that is to say, an extreme point is a point which is not an interior 
point of any line segment belonging to Kk. 

An equivalent definition of an extreme point is that x is an extreme point of 
a convex set K if K\ {x} is convex. 


8 


Definition 86. A metric space is a non-empty set X for which is defined a 
concept of distance. The distance d is called a metric on X, having such 
properties that, for any points x and y in X, we have d(z,y) > 0, and 
d(x, y) = 0 implies x = y; d(x, y) = d(y, x); and d(x, y) < d(x, z) + d(z, y). 


Let X be a metric space with metric d, let A be a subset of X and let x be 
any point of X. Then the distance from x to A is defined as 


d(x, A) = inf {d(x,a): a € A} 
whereas the diameter of A is defined as 
d(A) = sup {d(a1, a2) : a; and ag € A} 
Then a set is said to be bounded if its diameter is finite. 
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Further, let xg be a point in X and r a positive real number. Then the open 
sphere S,.(a9) with centre xo and radius r is the subset of X defined by 


S,-(ao) = {a : d(x, x0) < r} 


A point x in X is called a limit point of A if each open sphere centred on x 
contains at least one point of A different from x. A subset F of X is said to 
be closed if it contains all its limit points. 


8 


Definition 87. A linear space, aka a vector space, is a non-empty set L on 
which is defined two binary processes, say addition and scalar multiplication. 
Addition is defined such that for any x, y and z in L, then x+y is again in L; 
ety =ytau;e+(y+z) = (a@+y)4+2; there exists a unique identity element 
0, aka zero element or the origin, such that «+0 = «x for every x; and there 
exists a unique inverse element —x for every x, such that x+(—2x) = 0. Scalar 
multiplication is defined with regard to scalars, some instances of which are 
real and complex numbers, such that for any scalar a and any x and y in L, 
az is again in L; a(a+y) =axr+ay; (a+ B)x = ax + Bu; (aB)x = a(Bx); 
and lx = x, where 1 is the identity for scalar multiplication. A normed linear 
space is a linear space on which is defined a norm, that is a function which 
maps each element x in the space to a real number ||2|| in such a manner 
that ||z|| > 0, and ||a|| = 0 if and only if x = 0; ||a+y|| < |lz|| + |ly||; and 
lloa|| = Jo ||). 


8 
Theorem 43. A normed linear space is a metric space. 

8 
Problem 37. Prove Theorem 43. 

8 


Theorem 44. Any vector in a closed and bounded convex set with a finite 
number of extreme points can be expressed as a convex combinations of the 
extreme points. 


8 


Problem 38. Prove Theorem 44. 
§ 


Definition 88. For two vectors, that is points, « and y in R”, we write 
x > y if and only if 7; > y; for all 1 <i< n. A system of m weak 
linear inequalilties in n variables can be written as Ax > b, where A is an 
m xn matrix. A fundamental question concerning such system is whether it 
is consistent, that is to say, whether there exists some x such that Ax = b. 
A system may be inconsistent, or it may have a set of solutions which is 
unbounded. If we sketch our problem on a graph, we may see that it’s solution 
set is convex. 


58 26" April, 2006 God’s Ayudhya’s Defence 


kit Tyabandha, PhD Business Mathematics, notes and projections 


8 


Theorem 45. The solution space of a set of simultaneous linear equations 
is a convex set the number of extreme points of which is finite. 


8 


Problem 39. Prove Theorem 45. 
§ 


Theorem 46. Let S be the set of all feasible solutions to the linear pro- 
gramme in standard form in Definition 80, in other words, S is the set of all 
vectors x that satisfy Ax = b and x > 0, where A is an m x n matrix. Then 
S is a convex set, and the number of its extreme points is finite. The ob- 
jective function attains its optimum, provided that one exists, at an extreme 
point of S. If m <n, then the extreme points of S have at least n — m zero 
components. 


§ 
Problem 40. Prove Theorem 46. 


Algorithm 2 Procedure for finding basic feasible solutions. 


Input: Ax = b, A is an m x n matrix, m <n, rank A =m 
[a, -:- a,J<A 

(aja, +--+ + ana, = b) + (Ax = b) 

for i=m+1 ton do 

endfor 

(@1,-.--;Un) < solve x1a, +--+ +a2,a, =b 


Definition 89. The simplex method is a matrix procedure which solves 
linear programmes of the standard form as described in Definition 80 where 
b > 0. Starting from a basic feasible solution xo we locate successively other 
basic feasible solutions giving better values for our objective. For minimi- 
sation programmes the method uses Table 3, for maximisation programmes 
the same table is also used but with the sign of entries in the bottom row 
reversed. 


8 
Table 3 Table used for minimisation programming in simplex method. 
ee 
ar 
Xo Co A b 
cP —cdA —cib 
8 


Table 4 Description of the simplex method. 
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while negative number exists in d do 


Locate the most negative number in the bottom row of the simplex 
table, excluding the last column. The column in which we find 
this number is called the work column. If more than one such 
column exist, choose one of them. 

Find the smallest of the ratios between the elements in the last column 
and the elements in the work column of the same row, if these 
latter are positive. The element in the work column that yields 
this smallest ratio is called the pivot element. If there are more 
than one of these, choose one. If none of the elements in the work 
column is positive, the programme has no solution. 

Using elementary row operations, convert the pivot element to 1 and 
reduce all other elements in the work column to 0. 

Replace the x-variable in the pivot row and first column by the z- 
variable in the first row and pivot column. This new first column 
then becomes the current set of basic variables. 


endwhile 

The optimal solution is one in which all the basic variables assume the cor- 
responding values in the last column, while the remaining variables are zero. 
The optimal value of the objective function is then the value of the last row 
and last column for a maximisation programme, and the negative of this value 
if the programme is one of minimisation. 


8 


Algorithm 3 gives a procedure for the simplex method. Here M repre- 
sents a large positive integer, p a ratio, c a column, r a row, and xo contains 
all the basic variables. Also the last row in Table 3 is represented here by 
d=c? —cdA and e = —cdb. 

Algorithm 3 Algorithm for the simplex procedure. 

j <0 

while there exists a negative number in d do 

jgegtl 
for 1=1tondo 
{c’s} < (column number of the most negative number in the bot- 
tom row) 
(work column) «+ choose one of the {c}’s 
k < (work column) 
endfor 
Ppivot «+ M 
c+0 
soln + 0 
for 1= 1 to m do 
if (a;),, > 0 then 
soln + 1 
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(bj); 
(aj) in 
if p < Ppivor then 
ret 
endif 
endif 
endfor 
if soln = 0 then 
no solutions exist 
endif 
convert} A, such that (a;),,, = 1 and (a;),,=0,1<i<m,if#r 
(Xo), — ©k 
endwhile 
for i= 1 to m do 
(x0); — (b;); 
endfor 
for i= m-+1tondo 
x; +0 
endfor 
a e; 
if the programme is one of minimisation then 
ze —2* 


endif 


pe 


Definition 90. The two-phase method is a procedure modified from the 
simplex method to cope with cases when artificial variables exist in the initial 
solution xg, in order to minimise the round-off errors that occur in the cal- 
culation. The last row in Table 3 in this case is d = c? — cA A =d,+ Mdo, 
and consequently we have Table 5 which is used here. Algorithm 3 is then 
firstly applied to the last row, and then again to those elements directly above 
the zeros in that row. When an artificial variable is removed from the first 
column of the table, it ceases to be basic and may be removed from the top 
row of the table together with the entire column under it. When the last row 
contains only zeros, it may be deleted from the table. The programme has 
no solution if non-zero artificial variables are present in the final basic set. 


8 


Table 3 Table used for minimisation programming using the two-phase 
method. 


xl 
cl 

Xo Co A b 
d, —cpb 
d, 


{+ With the use of elementary row operations. 
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8 


Definition 91. Given a linear programme in the variables 71,...,2,,, there 
exists another linear programme associated with it, called its dual, which is in 
the variables w1,...,Wm. The original programme is called the primal. The 
primal completely determines the form of its dual. The symmetric dual of a 
primal linear programme in the matrix form 


minimise: z=c x 
subject to: Ax >b 
with: x>0 


is the linear programme 


maximise: z= b’w 
subject to: ATw<c 
with: w>0O 


The dual variables w,,...,wWm are known as shadow costs. The unsymmetric 
dual of the primal 
minimise: z7=c x 
subject to: Ax =b 
with: x>0 

is 

maximise: z= b’w 

subject to: A’w<c 


The unsymmetric dual of the primal 


maximise: z=c’ x 
subject to: Ax =b 
with: x>0 
is 
minimise: z= b’w 
subject to: A?Tw>c 


8 


Note 4. We may see from Definition 91 that the dual of a programme 
in standard form is not itself in standard form. These duals are said to be 
unsymmetric. 


8 


Theorem 47. If an optimal solution exists for either the primal or the dual 
programme, then the other programme also has an optimal solution. If the 
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duality is symmetric, then the two functions have the same optimal value. If 
the duality if unsymmetric, then the optimal value of each function can be 
derived from that of the other. 


8 
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Examples for linear programming 
26" April, 2006 


4. Use simplex method. 


maximise: z= 2, +9%2+2%3 
subject to: 2, + 2%. +323 <9 
321 + 2x2 + 2x3 < 15 


with: all variables non-negative 


Solution. 
Ly 1 
v2 9 
7_ fl 2 3 1 ~«0 _ {9 _ | &4 
aa a isc A=|5 22-0 1) B= [15] -% ial 
v4 0 
25 0 
cf. a 


optimise: z=c’ x 
subject to: Ax=b 


with: x>0 
apf 
on 
Xo Co A b 
+ (c? — cf A) Feb 


Note: (er - co A) and —cj b in case of a minimisation problem, whereas 
- (c? - co A) and cab in case of a maximisation problem. 


Tableau 1; 
U1 v2 v3 T4 Ls 
r4 1 2 3 1 0 9 
£5 3 2 2 0 1 15 
1 9 1 0 0 


The most negative number in the last row is —9. Therefore x2-column 


becomes the work column. And then, 


x2 
t 2 — positive — ee 
5 2 —+ positive > 3 =7.5 
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Since min(4.5,7.5) = 4.5, the value of x2 on the row corresponding to x4 
becomes our pivot element. Then carry out a series of elementary row op- 
erations, namely in order (1)2 © (I)1/2; UD)2 © (11) — 2(2)2; ID)2 — 
(IIT), + 9(D)2; 


Ly 22 x3 v4 XS 
1 3 1 9 
2 3 1 2 3 0 2 
as D 0 | 1 6 
7 25 9 81 
5 0 z 3 0 os 


Now the last row is all non-negative. Therefore x5 3, uz = 6, 27 
‘== Vent 2= = 
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Integer programming 
6°” December 2005 


Definition 92. Algorithms which change the boundary of the solution 
region in order to find the optimal solution of an integer programme are 
called cut algorithms. The branch-and-bound algorithm does this by splitting 
the solution region into two and then discard the one which does not contain 
the optimal solution. The Gomory algorithm, on the other hand, reduces the 
feasible region with the help of a new constraint without the region being 
splitted. 

8 


Definition 93. We call branching a process by which a programme whose 
solution contains a non-integral 7 < x; < k is made into two separate pro- 
grammes having the additional constraint x; < j in one, and 2; > k in the 
other, the objective together with all the constraints of the original problem 
of which remain the same. Here j and k are positive integers and j < k. 


8 


Definition 94. In the branch-and-bound algorithm, if the objective is max- 
imisation, the value of the objective obtained when the first integral approx- 
imation occurs is said to be the lower bound for the problem, and if the 
objective is minimisation it is said to be the upper bound of the same. 


8 


Algorithm 4 Branch-and-bound algorithm for integer programming. 


find first approximation 
while approximations not all integers do 

choose x; from all non-integral variables such that 

min (|; — |||, |e; — [a]|) is maximised 

branch 

choose the branch whose value of the objective is maximum 
endwhile 
solution¢last approximation 


Example 59. (Problem 6.9; Bronson, 1982) 


maximise: z= 2%, +2724 73 
subject to: 2x, + 349 + 3x3 < 11 


with: all variables non-negative and integral 
Solve by branch-and-bound algorithm. 


Solution. Draw a simplex table of Programme 1. 
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ZY Z2 x3 LA 
1 2 1 0 
v4 0 2 3 3 1 11 
-1 =2 -1 0 0 
Replace x4 for x2 as the basic variable. 
Ty x2 x3 LA 
2 a s 
1 2 22 
3 0 1 3 3 


Be > = 3.6, ef = 23 = 2, =0,2" = > Since 3-< 23 <4, branch 
into two programmes, namely Programme 1 where x2 < 3, and Programme 2 
where x2 > 4. Consider first Programme 2. 


maximise: z= 2%, +2724 73 
subject to: 2x, + 342 + 3x3 < 11 
rw <3 
with: all variables non-negative and integral 


Use the simplex method in a tabulated form. 


Ty r2 X3 TA U5 
1 2 1 0 0 
v4 0 2 3 3 1 0 11 
rs 0) 0 1 0 0 1 3 
-1 -2 -1 0 0 0 
Replace the basic variable x5 with x2 
T1 r2 x3 v4 v5 
LA 2 0 3 1 -3 
x2 0 1 0 0 1 3 
-1 0 -1 0 2 6 
Replace the basic variable x4 with 71. 
Cal 2 v3 v4 x5 
Ly 1 0 3 5 —3 1 
r2 0 1 0 0 1 3 
oo ¢ ¢G 7 
ej =1, 43 = 3, 73 = 2, = 23 = 0, 2* =7 Then consider Programme 3. 


maximise: z=2, +2724 73 
subject to: 2x, + 349 + 3x3 < 11 
o>) > 4 


with: all variables non-negative and integral 
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Draw a table for the two-phase method. 


v1 x2 x3 GA rs x6 
1 0 0 —M 
LA 0 2 3 1 0 0 11 
x6 —M 0 -1 4 
-1 -2 -1 0 0 0 0 
0 -1 0 0 1 -1 -15 
Change x4 for x2 in the basic variables. 
T1 r2 v3 v4 v5 v 
x2 2 1 1 3 0 0 a 
es ae ee ee a” | 1 
3 0 1 2 0 0 2 
3 0 1 3 1 -1 -% 
The coefficient parts of the row corresponding to the basic variable x¢ 
and the last row cancel each other. The optimal result is 75 = an LH a= 


| es eee ee * _ 22 
vy =@,=% =Oand 27 = =. 


(3) a > 4, 2" = %, (0,2) 


Algorithm 5 Gomory algorithm for integer programming. 


while solution not wholly all integers do 
choose one non-integral optimal approximation 
write a relation from the row where that variable is basic 
rewrite the relation to make all fractional coefficients 
some integer plus a proper fraction 
move all the fractions to LHS, and all the non-fractions to RHS 
write a new constraint as LHS > 0 
find the solution for the original problem together with 
the new constraint 
endwhile 


Example 60. (Problem 7.1; Bronson, 1982) 
maximise: z= 2%,+ 72 
subject to: 22, + 5% < 17 
321 + 2x2 < 10 


with: 21,22 non-negative and integral 
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Use cut algorithm. 


Solve 


Solution. Find the first approximation of Programme 1 normally using the 
simplex method. 


ey v2 x3 v4 
2 1 0 0 
a, 0 2 5 1 0 17 
v4 0 3 2 0 1 10 
-2 -1 0 0 0 


Since #29 < a we know that 3 is the pivot element, and therefore we 
replace the basic variable x4 with 71. 
Ly Ho) x3 L4 
11 2 31 
v3 0 a 1 3 3 
2 1 10 
1 2 20 
0 3 0 3 3 
x _ 10 eS 3) oS eS #2 520 s 
We have zj = 3, 13 = 3, ©} = 74 = 0 and z* = =. Since both 


xj and «#3 are non-integers, arbitrarily choose the former to generate a new 
constraint. Then our Programme 2 becomes, 


and our new programme becomes 


maximise: z= 2%, + %2+ 0x3 + 024 
11 2 31 


bject to: — 
subject to 3 2 gt 3 
es 10 
x sq - 5% = > 
Pt ee oe, 
2% +44>1 


with: all variables non-negative and integral 
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Uy x2 X3 vA x5 v% 
2 1 0 0 0 —M 
Gi 0 1 2 0 3 0 0 = 
v3 0 0 2 1 —3 0 0 a 
te —M 0 2 0 1 -1 1 1 
-2 -1 0 0 0 0 0 
0 -2 0 -1 1 -1 -1 


T1 t2 3 v4 U5 v6 
Ly 1 0 0 3 3 3 
3 0 0 1 . 7 i a 
m | o 1 of =} F 
1 1 1 1 
-2 0 0 3 ie 2 2 
0 0 0 0 0 0 0 
This becomes, 
Vy v2 x3 v4 x5 
Ly 1 0 0 0 3 3 
x3 0 0 1 -3 7 o 
2 a oe ar 1 
1 1 13 
0 0 0 5 é > 
Then our first approximation of Programme 2 is 7} = 3, x5 5, x3 a 
ri = vf = 0, and 2* = #8. Arbitrarily choose x3 to generate the new 
constraint. 
i 1 1 
t+ =%4- =%5 = = 
2 apts "5 
1 1 1 
L4 95 5) = —2@2 
1 1 1 1 
ee ee 
4— 5755s 
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Then our Programme 3 b 


maximise: 


subject to: 


with: 


We draw our table for thi 


Business Mathematics, notes and projections 


ecomes, 


z= 20, 4+ %2 + Org + 0x4 + Or5 


1 
co rd 
Bie: Ag. ke 
TIO Dee gee SD 
fe ee 
XQ gv gue 5 


%4—- 2X5 >1 


all variables non-negative and integral 


S programme. 


L1 x2 13 LA U5 x6 x7 
2 1 0 0 0 0 —M 
L1 0 1 0 0 0 5 0 0 3 
r2 0 0 1 0 4 4 0 0 4 
r3 0 1 0 0 -3 2 0 0 z 
U7 —M 0 0 0 1 -1 -1 1 1 
-2 -1 0 0 0 0 0 0 
0 0 0 -1 1 1 -1 -1 
Then x4 replaces the basic x7 to become the pivot element. 
Ly Zr v3 v4 5 v6 
Ly 1 0 0 0 ; 0 3 
re 0 1 0 0 0 3 0 
x3 1 0 0 0 — 3 - 3 11 
x4 0 0 0 1 -1 -1 1 
-2 -1 0 0 0 0 0 
Next, x, remains basic and becomes a pivot element. 
Ly x2 x3 x4 x5 6 
Ly 1 0 0 0 3 0 3 
re 0 1 0 0 0 3 0 
X3 0 0 0 0 -l - 3 8 
LA 0 0 0 1 -1 -1 1 
0 -1 0 0 3 6 
This becomes 
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Ty £2 x3 v4 v5 v6 
v1 1 0 0 0 $ 0 3 
re 0 1 0 0 3 0 
v3 0 0 0 0 -1 —3 8 
v4 0 0 0 1 -1 -1 1 
0 0 0 0 ; 5 6 


The optimum point for Programme 3 is then, xj = 3, 73 = 8, 73 = 1, 
v5 = x8 = x6 = 0 and z* = 6. Therefore the solution to the original problem 
Programme 1 is x} = 3, x3 = 0 at the objective value z* = 6. 


# 


Definition 95. A transportation problem involves m sources each of which 
supplies a;,i=1,...,m, units of a certain product, and n destinations each 
of which requires b;, 7 = 1,...,n, units of the same. The problem may be 
stated as following. 


maximise: z= 


subject to: So ti =a;, t=1,...,m 


with: all x,; non-negative and integral 


The total supply and the total demand are assumed to be equal. Were this 
not so, a fictitious destination or a fictitious source is added. 


8 


Definition 96. The north-west corner rule finds an initial basic solution for 
the transportation algorithm of the integer programming. It begins with the 
(1,1) cell in the m x n table, and allocates as many units as possible to 711 
violating neither the constraints of supply, that is the summation along each 
row, nor those of demand, that is the summation along each column. Then 
carry on moving for each step either right or downwards, until we reach the 
lower-right corner, 2mn.- 


8 


Definition 97. A loop, which is a sequence of cells in the table used for 
finding the solution in the transportation problem, has the following proper- 
ties. 
a. each pair of consecutive cells is on either the same row or the same column 
b. no three, or in fact any odd-numbered, consecutive cells lie in the same 
row or column 
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c. the first and the last cells are on the same row or column 
d. the path along the loop is self-avoiding, that is no cells appear more than 
once in the sequence 


8 


Algorithm 6 Transportation algorithm. 


while optimal solution not attained do 
find an initial, basic feasible solution using, for instance, the North- 
west corner rule 
let either uj = 0 or vj = 0 depending on whether the i**-row or the 
j*)-column 
has the maximum number of basic solutions 
find all u; and vj,7=1,...,m andj =1,...,n from uj+v,; = cj; for 
basic variables, 
and from ¢;; — uj — v; for non-basic variables 
improve the solution 
endwhile 


Note 5. In a transportation problem, optimal solution is achieved when 
Cij — Ui — Uj > O for all transportation costs per unit c;; of all non-basic 
variables. 


8 
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Financial mathematics 
13” December 2005 


Definition 98. The present value po or the principal is the amount initially 
borrowed or invested. The future value p, is the principal after a period of 
time ft. 


8 


Definition 99. Interest rates expressed per annum are called nominal rates, 
i. The annual percentage rate or effective annual rate ig is the equivalent 
annual rate of different interest rates variously compounded. 


8 


Definition 100. <A sequence is a list of numbers which follows a definite 
pattern. It is called an arithmetic sequence if each of its terms is obtained 
from the term immediately preceding it by an addition of a constant d, which 
is called the common difference. It is called a geometric sequence if each of 
its terms is obtained from the previous term by a multiplication of a constant 
r, the common ratio. 


8 


Definition 101. A series is the sum of the terms of sequence. It is called 
a finite series is one whose number of terms is finite, otherwise it is called an 
infinite series. An arithmetic series or arithmetic progression is the sum of 
the terms of an arithmetic sequence. Likewise a geometric series or geometric 
progression is the sum of the terms of geometric sequence. 


8 
Theorem 48. The value of the n‘® term of an arithmetic series is 
Tn =a+(n—1)d 
The sum of its first n terms is 
Ss, = 5 (2a+(n—1)d) 
8 
Problem 41. Prove Theorem 48. 
8 


Theorem 49. The n'" term of a geometric series is 
T, =ar™! 
The sum of the first n terms of it is 


Sn =a+ar-+s-++ar™! 


eb) 
~ [=r 
xed) 
~ pal 
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When the number of terms approaches infinity, the summation in cases where 
r < 1 becomes 


8 


Definition 102. <A simple interest is a fixed percentage of the principal 
paid to an investor each year. A compound interest is an interest paid on the 
principal plus any interest accumulated in previous years. 


8 
Theorem 50. The present value for simple interest is 
Pt = po(1 + tt) 
where 7 is the interest rate and ¢ the time in years. 
8 
Problem 42. Prove Theorem 50 
8 
Theorem 51. The present value in the case of compound interest is 
pt = po(l +i)! 
8 


Note 6. The interest may be compounded more than once a year, for 
example biannually, quarterly, monthly, weekly, daily, or continuously. Each 
time period is called a conversion period or interest period. The interest rate 
applied at each conversion is i/m, where m is the number of conversion periods 
per year. The number of conversion periods over t years is then n = mt. 


8 
Theorem 52. The present value at the end of n conversion periods is 
: n . mt 
i a 
Pt = Po (1+=) = Po (1+=) 
m m 
where all the variables and parameters are as previously defined. 
8 
Problem 43. Prove Theorem 52. 
8 


Theorem 53. When the number of compoundings per year becomes very 
large, the present value becomes 


Pt= poe" 


Proof. Since p; = po (1+ ye and limm-—oo (1 + a\™ = e’, we have the 
proof. q 
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Theorem 54. The annual percentage rate when compounding occurs m 


times per year is 
4 m 
coe 
m 


Problem 44. Prove Theorem 54. 
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Examples for financial mathematics 
26” April, 2006 


5. Given the demand function for a good as p = 1800—3q. Find the coefficient 
of point elasticity of demand when p is 300, 900 and 1200. Describe in words 
each of the results. Find the percentage change if the price of good increases 
by 10 per cent. 


Solution. The coefficient of point elasticity of demand is eg = —;*, and 


also €g = ae where Agg is percentage change in quantity demanded and Ap 


percentage change in price. 


Po 300 900 1200 
qo C800—Po) 500 = 300 200 
Ed = 3 = ; -1 “2 
lea| ce <1 1 al 
demand elastic unit elastic inelastic 
Ada eqAp -2% -10% -20% 

# 
Bibliography 
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Integral calculus 
10°” January 2006 


Definition 103. Let f(x) be a function, and let f’(x) be its derivative. The 
reverse process of differentiation is called antidifferentiation or integration. It 
gives us the original function, which is called the antiderivative or integral of 


f (2). 


8 
Theorem 55. Let c, n and & be constants. Then 
a. 
[ke =kxr+c 
b. 
i dx =x+c 

Cc ; 

n = n+1 sas 

fe ire ig +c n#-l 

d. 

[ota=mnzte, x>0 
e. 

[otar=tnlel +o, 0OAnr <0 
f ; 
[tae = re +c 

g. 

[rt@de =k f re)ae 
h. 

[U@sta@) = f reas f o(e)ae 

it 

[-te@dx=- [ taax 

8 


Definition 104. The approximation )7;"_, (f(x;) Az’) of the area under a 
continuous curve A is called a Riemann sum. That area under the curve is 
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8 


Theorem 56. Let F(x) be the integral of f(x). We call the fundamental 
theorem of calculus the expression. 


b 
/ f(a)dx = F(a)? = F() - F(a) 


§ 
Theorem 57. 
[ sea a ‘i ore 
b. ; 
f(e)de = F(a) — F(a) =0 
- c : b c 
| f(a)dx = ‘ f(a)da+ ] f(a)dz, ax<b<ec 
d. i, : ‘ 
[ teres [ oayde= f (f(@) + o(@)) de 
: 
[ ifade= kf Heoyae 
§ 


Property d of Theorem 57 is used to find the area between two curves. 


Theorem 58. The process of integration by parts is 


we have 


Bibliography 
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Examples for integral calculus 
26" April, 2006 


6. Find is x? dx, i (4 = im) dx and ar ada. 


Solution. 


# 

Consider f*, 2;dx = F(1) — F(-1). Here F(1) = 1 and F(-1) = -1 

But co = F(0+) # F(0~) = —oo, which means that f(x) = 1/zx? is not 

continuous at 0, which is a point in [—1,1]. Therefore the definite integral 
given does not exist. 


# 
Solution. (1) We know that if u is a differentiable function of «, then 
n+1 
pera ae +C 
n+1 
provided that n 4 —1. Let u = 2741, then du = 2zdz. And then, 
22 1 
Ve+1 | ae 
a pra 
3 
sue +C 
3 2 
as (2? a 1) 3 af C 
2 
# 
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Solution. (2) Let u = %/2? +1, then 3u?du = 2zdz. 


2z 3u? 
Se ype fp EG 
lyn i U 
=3 [udu 

ee 


+C 


8. Evaluate re 307/23 + Idx. 


Solution. (1) Let w= 2°+1. Then du = 3x7dz, u(—1) = 0 and u(1) = 2. 


And then, 


1 2 
i 327/23 + ldx = i Judu 
1 0 


Solution. (2) Again, substitute the u and du as above. 


indefinite integral, 


[seve vide = [ vadu 


2 
= Bui +e 


3 
2 


= 5 (8 +1) +C 


And then consider the corresponding definite integral, 


1 3 |t 
/ 3a? / a3 + 1dx = 2 (a3 +1)! 
—-1 


3 -1 
_ 4v2 
3 
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9. Find the average value of f(x) = V4 — «dz on the interval [—2, 2]. 


Solution. The average value required is, 


2 
os! /4— ede =" 
—2 


f(D) 2 
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Exercises for Integral calculus 


26%" April, 2006 


10. Find the definite integrals, 


J (202° + 324 — 6x3) dx, f sur de, f 8/2 —7dz, 
f 4e-3-5¢dt, and f 3073 da. 


11. Find the definite integral 


/ (2e%* _ 3e7**) dt 


given an initial, or a boundary condition F'(0) = 3. 


12. Find the values of the definite integrals, f, 5a dx, f 4e3 dt, Hes 6x~3, 
and f_,(—4)e~*dt. 


13. Find a firm’s total revenue r; function, given the marginal revenue func- 
tion rm = —.2%? — 1.32 + 500. 


14. Let the present value be p = a~" of the sum of money a to be received 


in the future when the interest is compounded continuously. Find the present 
value p,, of a stream of future income, that is to say, the money to be received 
each year for n years. 


15. Find the following integrals and definite integrals. 


[8x3 dx f? 803 dx f Vz + 2dx 
(jaxtde. “fat 2 3 as 
(0-227) de =f? (10-2e?) dx = ff 4(e@ — 3) 3 dx 
f?4(a—3)-3de ff 10e! dt JP 10e!-* de 
[Sat fi aa [ Fp42 

Sl az (tite) de IP (Jeti +2) ae 
f x(6x + 3)dx Ve x(6a + 3)da J (823 + 9x?) dx 

te (323 +907) dx f (x3 + 1) dx i (x3 + 1) dx 
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J 


S 
J 
J 
J 
S 
i) 
S 
J 
J 
J 
J 
J 
S 
J 
J 
J 
J 
J 


(5+.2)*dx 

: (308 + = +e") dx 
(1 _ et) dt 

2 (e® +e +2) de 
(ce?! + 34? + 4) dt 

* dt 

—5daz 

: a dau 

5e~ 7 dt 

* fda 

edt 

: (5x4 + + +e") dx 
(2ye-+2) a 

2 Jat ider 
(12Vz — 3) 
—V/ardx 

(4 + 2e*) dx 


2 30-lde 
2(2 —1)~*dax 


Cs SS ns Ra NS eR He 


*(5 +2)*dx 
2x(a + 3)?dx 

: (1 - e- 7) dt 

+ dx 

* (e* +30? + 4) dt 
furdu 

* _5dax 

xe* dx 

* Bet dt 

(6/x + 3) dx 

* esdt 

3e7! dt 

> (Qe +4) dt 
3./rdx 

> (12,/e — 1) de 
(2 +e”) dx 

‘ (+ + 2e) dx 


rip dz 
2, 


2(2 — 1)~*dx 
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‘ (30% t+ +e") dz 
He 2x(x + 3)? dx 
(Pets b) d 
fi tae 

f edt 

3 hg udu 

[daz 

* ret da 

side 

* (6Vr + 3) dx 
(5x4 + 4 + e) dx 
? 3e24 dt 

xa +7dx 

* 3,/xdx 

—Jfrdz 

: (2 + e”) dx 
3a-tdax 

2 


sa SR Se eS es Ss, RL eee 


2x 
ceri a 
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Integral calculus 
7” February 2006 


Definition 105. 
ac [Fede 0 
b. Hk f(z)dx =— fP f(a)dx 


Theorem 59. 
a. [Uk f(a)dx =k f° f(x)dx. And if k = —1, then 


[-1@e = — f feyae 


b. f? (f(x) £ g(x) da = f° f(x)da + f? g(a) dx 

c. If f(x) > g(x) on [a, bj, then rh f(a)dx > J gla)dx Let g(x) = 0. Then, 
f(x) > 0 on [a,b] implies ie f(x)dx >0 

d. If max f and min f are the maximum and minimum values of f on [a, b], 
then 


b 
min f-(6—a) < f f(x)dx < max f - (b— a) 


e. If f if integrable on the intervals between a, b and c, then 
b c c 
[feat [ pean =f pyar 
a b a 


Theorem 60 is called 60. 


Theorem 60. If f is continuous on the closed innterval [a,b], then at some 
point c in the interval [a, bJ, 


b 
fle) =p f Meade 


Theorem 61. The average- or mean value of an integrable function f on 


a, 5] is 
b 
— | f(x)dx 


Theorem 62. If f has a constant value c on [a, 6], then 


[ see = [cae =c0-«) 
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Theorem 63. If f is continuous on [a,b], then the function F(x) = 
J’ f(t)dt has a derivative at every point on [a,b] and 


dF 


d xz 
era (t= se) 


8 


Theorem 64. If f is continuous at every point of [a,b] and F is an an- 
tiderivative of f on [a,b], then 


b 
ih f(a)dx = F() - F(a) 


Bibliography 
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Simultaneous equations 
24" January 2005 


Definition 106. We call system of equations equations which together 
describe a mathematical model. A system of equations of n equations and vu 
variables is called an n x v system or a system with n x v dimensions. If n = vu 
the system of equations is called an exactly constrained system, ifn < v an 
under-constrained system, and if n > v an over-constrained system. 


8 


Note 7. A unique solution to a system exists only if there are as many 
equations as variables, that is to say, ifn > v. An under-constrained system 
may have an unlimited number of solutions or no solutions, but it may never 
have a unique solution. An exactly constrained or over-constrained system 
may have a unique solution, an infinite number of solutions, or no solution. 


8 


Definition 107. The graph of a (2 x 2) linear system of equations comprise 
two straight lines. If the two lines intersect, then the point of intersection 
(11, y1) satisfies both equations and therefore represents a unique solution of 
the system. If they do not intersect, then there are no solutions and the two 
corresponding equations are said to be inconsistent with each other. If the 
two equations have identical graph, then the system has an infinite number 
of solutions. Such equations are called dependent or equivalent equations. 


8 


Note 8. Consider a (2 x 2) system of linear equations in the slope-intercept 
form, 
y= mx + by 


y = max + be 


if my x msg then 
system has a unique solution 
else 
if by x bg then 
equations are inconsistent and the system has no solution 
else 
equations are equivalent and the system has infinitely many solu- 
tions 
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Differential equation 
31°* January 2005 


Definition 108. A differential equation (DE) is an equation which involves 
derivatives. An ordinary differential equation (ODE) is a differential equation 
in which there is exactly one independent variable. A partial differential 
equation (PDE) is one where there are at least two independent variables. 
The derivatives of an ODE are ordinary-, whereas those of a PDE are partial 
derivatives. 


8 


Definition 109. Consider a differential equation. The order of it is the 
order of the highest derivative appearing in it. Its degree is the degree of 
the highest ordered derivative therein. A primitive is a relation between the 
variables that involves n essential arbitrary constants, which gives rise to a 
differential equation of order n. The n constants are called essential if they 
cannot be replaced by a smaller number of constants. 


8 


Example 61. The differential equation y’” + 3(y”)? + 2y’ = sina is an 
ordinary differential equation of order 3 and degree one. The differential 
equation (y!")” + (y’)® + y = 2x is an ODE which has an order 2 and degree 
2. 


Problem 45. The problem of finding solutions of differential equations is 
essentially that of finding the primitive which gave rise to the equation. 


8 


Example 62. The differential equation y’” = 0 has a primitive y = Av? + 
Bet+C, y”—6y"+11y'—6y = 0 has y = Cye3* + Ce?" +C3e®, y? (y”) +y? = 
r? has («x -—C)? + y? =r’. 

8 
Definition 110. Existence theorems give conditions by which one could 
determine whether a differential equation is solvable. A particular solution 
of a differential equation is one obtained from the primitive by assigning 
definite values to the parameters, that is to say, the arbitrary constants. A 
singular solution is a solution which cannot be obtained from the primitive by 
any manipulation of the arbitrary constants. The primitive of a differential 
equation is usually called the general solution of the equation. 


8 


Definition 111. A differential equation is said to be variable separable if 
an integrating factor can be readily found. Such equation has the form 


fo(x) + g2(y)da + fo()- gi(y)dy = 0 
Through the use of the integrating factor 
1 
fa(x) - g2(y) 


8&8 26" April, 2006 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


the primitive of this is then 


A(z) 4... fH 5 
nytt | dy =C 


8 


Definition 112. A differential equation of the first order and first degree 
may be written in the form 


M(x,y)dz + N(x,y)dy =0 


If this such equation admits a solution f(z, y,C) = 0 where C is an arbitrary 
constant, then there exist infinitely many integrating factors xi(x, y) such that 
&(x, y) [M(a, y)da + N(x, y)dy] = 0 is exact, and there exist transformations 
of the variables which render the latter separated. But since no general rules 
exist for doing this, the use in practice is still somewhat limited. 


8 
Definition 113. A function f(z, y) is said to be homogeneous of degree n 
if 

f(Ax, Ay) = A" f(x,y) 

8 

Note 9. The equation 
(aru + bry + c1) dx + (agx + boy + co) dy = 0 
where a bz — agb; = 0, is reduced through the transformation 
dt —a,d 
ae+by=t and dy= eal iad a 
by 
to the form 
P(«,t)dxz + Q(az,t)dt =0 
8 


Note 10. The equation 


(aya + bry + C1) da + (anu + boy + cg) dy =0 
where a,b — agb; 4 0, is reduced through the transformation 
c=a'+h and y=y'+k 


in which x = h and y = k are the solutions of the equations a;2+b;y+c,; = 0 
and agx + boy + co = 0 into the homogeneous form 


(a,ax' + byy’) da’ + (agx’ + bay’) dy’ = 0 
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8 
Note 11. The equation of the form 
y: f(ry)dx + x - g(xy)dy = 0 
through the transformation 
z adz—zdx 
RUS ey. Sy 
x x 
into the form 
P(a,y)dx + Q(a, z)dz =0 
which is variable separable. 
8 
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Examples for differential equations 
26" April, 2006 


16. Solve x?(y + 1)dxz + y?(x — 1)dy =0. 


Solution. The integrating factor is ( by the multiplicative applica- 


1 
y+1)(x-1) 
tion of which the given equation becomes 


2 2 


y 
d 
a el 


dy =0 


And then from this, 


2 y? 
—+a+4In(x—-1)4 5 y+In(yt+1)=Ci 
x? + 2a —2y4+In(# —1)(y+1)=Cy 
(2 +1)? +(y—1)? +2ln(e-Dy+)=C 
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Difference equation 
20" February 2005 


Definition 114. A difference equation gives the relationship between an 
independent variable and a dependent variable, which changes at fixed, equally 
spaced intervals in time. The order of a difference equation is the number of 
time intervals spanned within the equation. 


8 


Example 63. In the difference equation y+; = 1.ly:, an independent 
variable is the time t while the dependent variable is the income y. The order 
is the span of t intervals within the equation, which in this case is (t+1)—t = 1. 


Theorem 65. ‘The general solution of a homogeneous first-order differ- 
ence equation is of the form y, = Aa‘, where t and y are respectively the 
independent and dependent variables. 


8 


Note 12. Given a homogeneous first-order difference equation of the form 
Yt+1 — by, = 0 and some conditions, we may find the parameters a and A of 
our general solution y, = Aa’ by first finding a by substitution of this general 
solution for t and t+ 1 in the difference equation, and then find A from the 
conditions given. 


8 


Theorem 66. The stability of the solution to a difference equation in 
general form y; = Aa’ is, 


range of a time path of y, solution time path 
-o<a<-l at's +0 unstable alternates 
-l<a<0 a’ +0 stable alternates 
0<a<l a’ +0 stable tends to zero 
l<a<aw a’ + co unstable tends to infinity 


8 


Theorem 67. The solution of a non-homogeneous difference equation is the 
sum of a complementary function and a particular integral, that is y = ye+Yp- 
The complementary function is the solution of the homogeneous part of the 
difference equation. The particular integral is a function which satisfies the 
full difference equation. 


8 


Note 13. The general form of the particular integral is deduced from the 
right-hand side of the difference equation. In particular, if c and b are con- 
stants, 
right-hand side general form of particular integral 
c Yp =k 
cbt i= hv 
8 
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Examples for difference equations 
26" April, 2006 


17. Solve the difference equation y:41 — 1.ly, = 0 by iteration for year 2, 3, 
4 and 5, given the income in year 1 is 20,000 Bahts.. 


Solution. We have y41 = 1.ly:, therefore 


yi = 20000 

yo = 1.1(20000) = 22000 
yg = 1.1(22000) = 24200 
y4 = 1.1(24200) = 26620 
y5 = 1.1(26620) = 29282 


# 


18. Write out the solution of the difference equation y+, — 1.ly, = 0 for t = 
1, 2, 3, 4 and 5 in terms of y,. Deduce the general expression for y; in terms 
of y;. Evaluate y4o given y, = 20000. 


Solution. 
Y2 = L.lyy 
y3 = Lye = 1.1 (1.1y1) = (1.1)?y1 
ya = 1.1y3 = 1.1(1.1)?y. = (1.1)? 
Ys = 1Llyg = 1.1(1.1)?y1 = (1.1)*y4 
# 
In general, 
ye = (1.1) *y2 
# 
t = 40, y: = 20000; 
yao = (1.1)3°20000 = 822895.56 
# 


19. Find a general solution of the difference equation y:4; — 0.9y, = 0. If 
yo = 100, find the particular solution. Evaluate y1, yo, y3, yao and yso. 


Solution. The general form of the solution is y, = Aa’. Therefore y41 = 
Aa't!, Then the difference equation becomes 


Aa‘ (a — 0.9) = 0 
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Since A # 0 and a! # 0, therefore a — 0.9 = 0, that is a = 0.9. Thus the 
general solution is 


Ut = A(0.9)! 
# 
From y2 = 100; 
Y2 = A(0.9)? 
100 = A(0.9)? 
A = 123.46 
The particular solution is then 
Yp = 123.46(0.9)* 
# 
Then we tabulate the required calculation, 
t Ut 
1 123.46(0.9) = 111.11 
2  123.46(0.9)? = 100.00 
3 123.46(0.9)? = 90.00 
20 123.46(0.9)?° = 15.01 
50  123.46(0.9)°° = 0.64 
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Exercises for difference equation 
26" April, 2006 


20. For each of the following difference equations, state 
i The order of the equation. 
ii Whether the equation is homogeneous or not. 
(a) Pepi — 0.8p; = 0 (b) ye+o2=8- yest (c) ye+o2 = 80+ yt 
21. Solve each of the following difference equations for the indicated variable 
by the iteration method. 
a. Yer — 0.8y: = 10, given y; = 1. Then find ys. 
b. pepo = 4pe41 — 8p4, given p, = 20, pp = 18. Then find ps. 
c. p, = 0.6p4_1 + 80, given p; = 100. Find ps. 
22. An amount of money A is invested at r% compounded annually. 
a. Show that the value of the investment at the end of any year ¢ is given 
by the difference equation p41 = (1+1r)p; 
b. Show that the general solution of the equation is py = po(1 +7)’. 
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Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Function 


A function of one independent variable is a relation in the form y = f(x) such 
that there exists one and only one value of y in the range of f for each real 
number x in the domain of f. The variable y is called the dependent variable. 
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Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


An implicit function is a function in which both dependent- and independent 
variables appear on the same side. 


An explicit function is one where the dependent variable is on the left hand 
side-, and the independent variable on the right hand side of the equation. 


Business mathematics, Graph and derivative, 25'” October 2005 -2- From 21° 
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Polynomial 


Anx” +... +a,% +49 
Polynomial equation 
Anx” +...+a;2+a9 =0 

Quadratic equation 


ax? +b2+c=0 


_ —b+ Vb? — 4ac 


2a 
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Definition of limits 


If f(a) is a function which draws closer to a unique finite real number / for all 
values of x as the latter draws closer to a, but « 4 a, then / is called the limit 
of f(a) as x approaches a. In notation this is, 


lim f(x) = 1 


na 
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Definition of limits 


For a function f(x), lim,-,, f(a) = 1 if and only if for every € > 0, there exists 
5 > 0 such that |f(a) — l| < € whenever 0 < |x —a| < 6. 
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Rules of limits 


lmk=k 
za 
lim 2” = a” 
«za 


lim kf(x) = k lim f(a) 


wa «wa 


lim [f(x) + g(x)] = lim f(x) + lim g(x) 


wa «wa wa 
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lim [f(z) - g(v)] = lim f(a) - lim g(a) 


Za za Za 
. f(x) lim,-+a f(x) F ; 
1 = ded that 1 
ea in) ee 
Jim (f(0))” = [Jim f(a)|" provided that n > 0 
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Definition of a derivative 


Let y = f(x). Then, the derivative of y with respect to = is, 


dy. 3... tae hy — Te) 
es h 


The various notations for the derivative include 


dx’ i a), y’, Dy and D(f(z)) 
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General rules of differentiation 


Let u, v and w are functions of x, and c is a constant. 


d 
dx 
d 

dx 


(c) =0 


(cz) =c 


d 
—(ca”) = nex” 


dx 


a ) du du 
ao da dx 


al 
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From 21° 


October 2005 , as of 26°" April, 2006 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 
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dy dydu 
dx dudz 
du 1 
Hy ~ 
dx = 
dy 
dy du 
d 
dx oe 
de" 
da 
dinzx 1 
dx  « 
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Laws of exponents 


Let p and q be real numbers, a and b positive numbers, and m and n positive 
integers. Then, 


a? - at = aPt4 
Pp 
a = 
— = qP-4 
aq 


(a?)4 = a4 
a° = 1, provided that a 0 
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pt 
aP 
(ab)? = aPbP 
Ya=an 
Yam =an 
aft _ Ya 
b Wb 
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Laws of logarithms 


log, mn = log, m+ log, n 


log, a log, m — log, n 
n 


log, m” = plog,m 


log, n 


] — 
eet log, a 
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global maximum 
uae 


7\ i|/ 


| 


Department of Mathematics, Mahidol University 


—e|- 


Lee 
global minimum 
(a) 


(b) 
function. 


Figure 1 (a) global maximum-, and (b) global minimum points of a 
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Figure 2 Local minimum- and local maximum points of a function. 
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y’=0 
y’<0 


local maximum 


re 


ee 


local minimum 
y’=0 
y’>0 
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concave up 


curves down 


y’> O 


y’<0 


concave down 


curves up 


(a) (b) 


Figure 3 The two curvature types, namely (a) concave up (y” > 0) and, 
(b) concave down (y” <0). 
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y’>0 
y’<0 y’>0 oOo > 
y’<0 
(b) 


(a) 


Figure 4 The four possible curvatures in two dimensions, considering 
both y’ and y", (a) y” > 0 and, (b) y” <0. 
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— a 
fon 2 \ 
. y’=0 4 


Figure 5 Inflection points, where y" = 0, (a) with y” increasing and, 
(b) with y"” decreasing. 
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y’=0 | y’=0 
y>0 
y’=0 


y’’<0 


(a) (b) 
Figure 6 Stationary inflection points, where both y’ = 0 and y” = 0, (a) 
with y” increasing and, (b) with y” decreasing. 
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Figure 7 shows a plot of the case where a = 1, b=2,c=3 andd=4. 


y =ar* + br? +er +d, where a=1,b=2,¢=3 and d=4 


Figure 7 The cubic function y = x3 + 2x7 4+ 3a 4-4. 
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Figure 8 shows the case where a= 1, b=5,c=4andd=3 


y=ax + bx? +cex-+d, where a= 1,b=5,¢=4 and d=3 
15 


toh E ff 


Figure 8 The cubic function y = x3 + 5a? + 4x 4+ 3. 
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Figure 9 The hyperbolic function (% — 1)(y — 2) =3. 
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The general demand function is of the form 


dd = Lp Oh, Dexa; tay he i) 


where qq is the quantity demand of good z, p the price of x, y the income of 
the consumer, p, the price of substitute goods, p, the price of complementary 
goods, t, the taste or fashion of the consumer, and a the advertisement level. 
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In its simplest case where all other factors are constant, the demand equation 
takes the form 


P=C1 — C2dd 


where p is the price-, while gg the quantity demanded of the good x, and c; 
and cz are positive constants. 
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The general supply function is of the form 


ds = f(p, C, Pos te; N,0,.. fl 


where qs is the quantity supplied of good 2, p the price of x, c the cost of 
production, po the price of other goods, t. the available technology, n the 
number of producers in the market, and o other factors, for example tax and 
subsidies. 
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The simplified relation for the supply is 


p=C1 + C2qs 


where q, is the quantity of x supplied, and c; and cg are positive constants. 
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Definition 15 Total cost and revenue 


The total cost c, comprises a fixed cost and a variable cost, that is 


Ce = Cf + Cy 


The total revenue r; is price times output, 


Tt = pd 
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Definition 16 Marginal cost and revenue 
The marginal cost cj, is the change in total cost caused by the production of 
an additional unit. 


The marginal revenue rm is the change in total revenue coming from the sale 
of an extra good. In other words, 


d 
Ci = a and 
Cs dr; 
eo dG 


where q is the output. 
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Definition 17 Average cost and revenue 


The average cost Cq is the total cost per unit output, 


C= 


q 
The average revenue rq is the total revenue per unit output, 


Le 
q 
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From Definition 15, r, = pq, and from Definition 17, rz = raq, therefore 
P—-Ta 
From Definition 17, cg = a and from Definition 15, cg = cf + cy. Therefore 
Ca = Caf + Cav 


where the average fixed cost caf = a and the average variable cost cg, = as 
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Procedure 1 Total-, average-, and marginal graphs 


Given: c;(2) 
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for each f € {cz,Ca,Cm} do 
find the critical values x. for f. =0 
find f” 
n © |Xe| 
for i=1tondo 
if f’ (af) > 0 then 
f(x) is convex and is the relative minimum of f 
elseif f’ (xf) <0 then 
f(x) is concave and is the relative maximum of f 
endif 
endfor 
find inflection points xy from f” =0 
endfor 
plot c:(x), then ca(x) and cm(x) 
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The relationship between input and output is called a production function, 


q=f(l,k,1,te, s,e,...) 


where / is labour, & phical capital such as buildings and machines, r raw 
materials, t. technology, s land, and e enterprise. 
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Assuming a short period of time, then / becomes the only independent variable 
and the other remaining factors are parameters, that is fixed, and therefore 
q=f(l). 


Then the marginal product of labour is 


dq 
Pim dl 
and the average product of labour is 
Pla = 4 
la 1 
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The marginal propensity to consume is Dem = ao 
The marginal propensity to save is Pom = ae 


The average propensity to consume is Dea = a 
The average propensity to save is Psa = me 
Here y is the income, c the consumption, and s the saving. 
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The profit is 
T=Tt— Ct 


At the break-even point 


that is 
Te = Ct 
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Calculus of multivariable functions 


Definition 21 functions of n independent variables 


A function y = f(x1,...,2n) is called a function of n independent variables 
if there exists one and only one value of y in the range of f for each tuple of 
real number (21,...,2n) in the domain of f. Here y is called the dependent 
variable while x;,1=1,...,n, the independent variables. 
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Definition 22 partial derivatives 
Let a multivariable function be 
y =f(a1,...,2n) 
The partial derivative of y with respect to x;, where 1 <i <n, is a measure 


of the instantaneous rate of change of y with respect to x; while x; is held 
constant for all 7 #71, where 1 < 7 <n. This partial derivative is defined as 


Oy i f(...,0; + Aa;,...) — f(a1,...,2n) 


= 1m 
Ox; Az;70 Ax; 


and can be written in either one of the following forms. 
dy at 
Ox,” Ox; 


» fx; (£1,---) Ln), fe; OF Ya, 
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Theorem 9 product rule 
Let z = g(x,y) - h(a, y). Then, 


and 
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Theorem 11 generalised power function rule 
Let z = [g(z,y)]". Then, 


Oz a: Sed Og 
On Ox 
and 5 3 
oe n-1, 98 
Te: Oy 
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Definition 23 second-order partial derivatives 


Let z = f(x,y). Then, the second-order direct partial derivatives are 
By (ced, emer Mm ee 
Ox \ Ox Oy \ Oy 


O? oe? 
fre, (fea, a and respectively fyy, (fy),, Ba 


These are also written 
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The cross partial derivatives are 


a (82) yy 2 (a 
Oy \ Ox MO \ By 


These are also written as 


oO? oO? 
legs (ea oe and respectively fyx, (fy)2; FED 
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Theorem 12 critical points 


For a multivariable function z = f(x,y) to be a relative maximum at (a,b) 
necessarily f,, fy = 0, and fer, fyy <0 and fre -fyy > (fry)? at that point. For 
the same at the same to be a relative minimum, necessarily fz,f, = 0, and 
fre, fyy > Oand frx+fyy > (fry)? there. Moreover, an inflection point is a point 
(a,b) at which fix -fyy < (fey)?, and both f, and fy, have the same sign. On 
the other hand, a saddle point is a point (a,b) at which fy, -fyy < (fry)?, but 
fr~ and fy, are of different signs. 
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Procedure 2 Procedure for determining a critical point of a function 
with two independent variables 


Given z = f(x,y) and a point (a,b), at this point, 
if f; = 0 and fy = 0 then 
(a, b) is a critical point 
if fox -fyy > (fey)? then 
if fr2 <0 and fy, <0 then 
(a, b) is a relative maximum of z 
elseif fr. > 0 and fy, > 0 then 
(a,b) is a relative minimum of z 
else + 
endif 
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elseif fx - fyy < (fry)? then 
if fr. -fyy > 0 then 
(a, b) is an inflection point 
elseif fx -fyy <0 then 
(a,b) is a saddle point 
else ¢ 
endif 
else 
test inconclusive 
endif 
else 
(a, b) is no critical point 
endif 
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Problem 15 detatls in the critical point procedure 


There are two dead ends in Procedure 2. The first one (fT) is the case where 
Ege pgp (eng) sand, either (ign: = O.  = 0) he OS Oy ea 
O,fyy > 0), (fer < 0,fyy = 0), (fer > O,fyy = 0), (fer < 0, fyy > 0), or 
(fra > O,fyy < 0). The second one ({) is where fez - fyy = 0. Find out 
what happen in these cases, and thus complete the missing lines of logic in 
Procedure 2. 
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Definition 24 derivative and differential 
By derivative dy we mean an infinitesimally small change in y with respect 
to an infinitesimally small change in x. By differential dy and dx we mean 


an infinitesimally small change in the values of y and respectively «x. 
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Example 13 derivative and differential of functions of one variable and 
two variables 


For a function of one variable y = f(x), the total derivative is 


dy 
dx 


dy 
dy={— Jd 
(33) : 


For a function of two variables z = f(x,y) partial derivatives are, the first- 
order partial derivatives 


and the differential of y is 
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and the second-order partial derivatives 
Or z rz o? Orz 


a2 = Zr, Oy ane Gok = Zpy and da5y = Zyex 


The total differential of z is 


of Of 
dz = (=) dx + (5) dy 


and for small changes which are not infinitesimal, dr becomes Az and the 
incremental change formula is 


Of of 
Azz ({—])Ar+{[—]A 
( Ox ( a ‘ 
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Definition 25 general production function 


The general production function is q = f(1,k), where q is output of the pro- 
duction, | labour and k capital. The Cobb-Douglas production function in its 
general form is 


q = al*k? (1) 


where a is a constant andO<a<1,0<6<1,l>0andk>0. 
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Example 14 Cobb-Douglas production function 


With the Cobb-Douglas production function, the marginal product of labour 
is, 


() 
Pim = 4 = A = alk (2) 


and the marginal product of capital 


From this we see that pim > 0 and ppm > 0. 
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Theorem 13 law of diminishing returns to labour 


From the Cobb-Douglas production function we have the law of diminishing 
returns to labour, which states that qu < 0. 
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Proof. From Equation 1 in Definition 25, 


2 
O*q  O (3) — Pim _ @x=1j%! 


A 3 BECO al P 


Since 0 < a < 1, therefore qy < 0. 
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Example 15 changes in marginal product values 


Using the Cobb-Douglas production function, 
Gh: = di = aepl-Tke* 


Therefore, qi, > 0 and qx; > 0. In other words, pj increases as capital input 
k increases, and respectively pym increases as labour input / increases. 
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Example 16 average functions of labour and capital 


For the Cobb-Douglas production function in Equation 1 the average product 
of labour is 
Pia = 5 = alk? (4) 


and the average product of capital is 


Pia = 7 = alk! (5) 
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Example 17 marginal functions of labour and capital 
Again using the Cobb-Douglas production function of Equation 1, the margin- 
al product of labour is 


0 
Pim = aI = aal®1kF (6) 


and the marginal product of capital is 


0 
Phim = 5 = apleke! (7) 
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Example 18 comparison between marginal and average functions 
From the APL equation, Equation 4, and the MPL equation, Equation 6, 
and since 0 < a < 1, therefore pm < pig. Similarly from the APK equation, 
Equation 5, and the MPK equation, Equation 7, since 0 < 6 < 1, we have 
Pkm < Pka- 
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Example 19 conditions for using labour and capital 


A producer likes to have a positive marginal function, which means that the 
productivity increases as the input increases. But the second derivative is 
negative, which means that this rate of increase slows down as time goes by. 
In practice, the conditions for using labour are, 


_ oq dpim _ 0°q 
Pim = al > 0, aa = ae << 0, and pim < Pla (8) 


The conditions for using capital are similarly, 


_ oq Cian: Og 
Phim = oe > 0 ag = Bp 


<0, and prm < Pra (9) 
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Definition 26 production function graphs 


An isoquant is a graph in two dimensions, k = f(l), plotted to represent a 
production function gq = f(1,k). The slope 


dk 
dl 
is called the marginal rate of technical substitution. The value of this slope at 
(Jo, ko) is denoted by 
dk 
dl |i ko 
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Theorem 15 slope of an isoquant 


The slope of an isoquant is the ratio of the marginal products. 
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Proof. The total differential of q = f(1,k) is 


_ (0q Oq 
dq = (3) dl + (32) dk 


Along any isoquant, dq = 0, therefore, 
oq oq 
= ({— }dl — ) dk 1 
0 (3) + (Fa (10) 
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This directly yield, after some manipulation, 


Or, from Equation 10 together with Equation’s 6 and 7, it follows that, 


dk = Plm 
dl = Pkm 


Business mathematics, Calculus of multivariable functions, 1** November 2005 —26- 
From 232 October 2005 , as of 26’ April, 2006 


128 26" April, 2006 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 27 returns to scale 


In the Cobb-Douglas production function equation, Equation 1, let both in- 
puts / and k change by the same proportion, and let be the constant of this 
proportionality. Then gz = a(Al)%(Ak)°, which leads to gg = A°+?q,. When 
a+ =1, the case is described as constant returns to scale, when a+ 3 <1 
as decreasing returns to scale, and when a+ 8 > 1 as increasing returns to 
scale. 
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Definition 28 homogeneous Cobb-Douglas production function 


The homogeneous Cobb-Douglas production function of order r is, 
f(Al, AK) = A’L(L, k) 


where r = (a, {). 
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Definition 29 utility function 


A utility function expresses utility as a function of goods consumed. In its 
general form this is, 


u= f(x,y) 


where x and y are the quantities of goods X and respectively Y consumed. 
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Definition 30 Cobb-Douglas utility function 
The Cobb-Douglas utility function is in its general form, 


u=axry 


where a is a constant, and0O<a<1,0<6<1l,x>Oandy>0. 
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Definition 31 marginal utility 
The marginal utility for a utility function with one variable, u = f(x), is 
qe = Ur = Urm- The marginal utility for a utility function with two variables, 
u = f(x,y), is de = Ur = Uem and ou = Uy = Uym.- 
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Definition 32 indifference curves 


The indifference curve is a graph y = f(a) drawn to represent a utility function 
u = f(x,y). Its slope dy is called the marginal rate of substitution. Setting 


the total differential equal to zero, 


Ou Ou 
0=du= (5) dx + (5) dy 


we find 
ia se 
dx Uy 
and 
dy Urm 
da Uym 
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Definition 33 partial elasticities of demand 
Let a demand function be 
da = f(a, Y, Po) (11) 


where qq is the quantity demanded of good a, pg the price of a, y consumer’s 
income, and py, the price of another good b. Then, the price elasticity of 
demand is, 
_ 94a Pa 

OPa Ia 


Ed 
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The income elasticity of demand is, 


And the cross-price elasticity of demand is, 


_ 94a Po 
Opp da 


c 
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Example 20 partial elasticity with respect to labour 


With the demand function as in Equation 11, the partial elasticity with re- 


spect to labour is, 
Oq | 
“a OL g 


And from Equation’s 6 and 4, this leads to, 


For the Cobb-Douglas production function, Equation 1, then ¢4; = a. 
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Example 21 partial elasticity with respect to capital 


Again, with the demand function as in Equation 11, the partial elasticity with 
respect to capital is, 


Oqk 
fa = Bie g 
Then, from Equation’s 7 and 5, 
pee Pkm 
k= 
. Pka 


For the Cobb-Douglas production function, Equation 1, we have égx = £. 
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Definition 34 function 


A function is an operator or a procedure which accepts a permissible input 
and transforms it into a unique output. The input is some nonempty set. If 
a function is defined to be 


y = fx) 


then x is the input vector, y the output, and f(-) the function itself. 
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Example 22 function 
If f(-) is the function of dressing, then its input is possibly a person and its 
output a dressed person. 


If f(-) is the function of making up, then the input is perhaps a girl and the 
output a made-up girl. 
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Definition 35 variables and parameters of a function 
Let 
y = f(a,x) 
be a function, where a is a set of all its parameters, and x a set of all its 


variables. Then y is its dependent variable and 2;, for all 7 € x, are its 
independent variables. 


In other words, x; vary, y follows, and a; could assume any value within the 
range of its permissible ones, but its value must be constant. 
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Definition 36 inverse function 
An inverse function is an expression of the independent variable in terms of 
the dependent variables. The inverse of the function f(-) is denoted by f~1(-). 
If f(-) is a function which admits one independent variable, namely x, then 
one could express it as, 
y = f(x) (12) 
Its inverse function is then, 
f(y) =a (13) 
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Example 23 operator and inverse operator 


Both the function and its inverse may be thought of as being an operator 
operating on an input to produce an output. The function, 


y = f(x) 


is understood diagrammatically as, 


yo] f(-) |—2@ 


while its inverse function, 
x =f '(y) 
is displayed as a diagram as, 


Gh) EAG). jae 
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Theorem 16 inverse function 


An inverse function must always be a one-to-one mapping. 
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Proof. Let f(-) be a function. Then f(-) can be either one-to-one or many- 
to-one, and therefore f~1(-) could turn out to be either one-to-one or one-to- 
many. But since f~! is also a function, so for each of the values in its domain 
the corresponding value in its range must be unique. This means that in 
cases where f~! turns out to be one-to-many, some constraints must be put 
on its input in order to make the output one-to-one, which then makes all the 
outputs from f~+(-) one-to-one. q 
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f(-) £*() 


addition subtraction 
multiplication division 
power root 
exponential logarithm 


Table 1 Some of the functions and their corresponding inverse functions. 
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v:a a 
x Ve 
a” log, x 


Table 2 The notational forms of functions and their inverses. 


in which division and logarithm are both undefined for a = 0. 
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Definition 37 some inverse functions 


The inverse of the addition, 


y=ut+a 
is the subtraction, 
y-a=2 
The inverse of the multiplication, 
y = ax 
is the division, 
y_ 
== 2 
a 
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The inverse of the power, 


is the root, 


The inverse of the exponential, 


is the logarithm, 


loggy =a 
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Example 25 building-blocks of mathematics 


Figure 10 addition makes multiplication makes power function 


LeL=2x 
t 
CPA SAS 
wes au-a 
a 
f 
ct+ue=2x 
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a:-:a 
—_=_ > a” 
x 
a 
at+...+a@ 
—_ SS a:n 
n 
t 
a+a= 2a 


Figure 11 Starting from a constant to obtain in the end the exponential 
function. 


Business mathematics, Exponential, log and nonlinear functions, 8°” November 2005 
-13- From 28" October 2005 , as of 26" April, 2006 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 38 exponential function 


An exponential function is defined as 
ya 


wherea>OandaFl. 
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Example 26 exponential function 


The domain of the exponential function y = a” is the set of all real numbers, 
while its range the set of all positive real numbers. The function is convex 
and increasing when a > 1, and convex and decreasing when 0 < a < 1. At 
x = 0, the value of the function is y = 1 for any a > 0. 
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Theorem 17 exponential to the power of zero 


For any a 4 0, 


lim a” = 1 
x0 
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Theorem 18 rules of exponential function 


Three basic rules of the exponential function are, 


a™a” pe qrtn 
m 

a =—qnr-n 
a” 


Business mathematics, Exponential, log and nonlinear functions, 8°” November 2005 
-17- From 28" October 2005 , as of 26" April, 2006 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Proof. Write, say, a” as, 


Sen 


m 


and similarly for a”. Then all three equations above become obvious. q 
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Example 27 graphs of the exponential function 


y = 28 
8 


0 L Hl i 
=3 =2 a 0 1 2 3 


Figure 12 Example of the graph of the exponential function when a > 1. 
Here the graph is that of y = 2”. 
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y = (0.5)" 


° i i i hi 
=3 =2 a 0 1 2 3 


Figure 13 An example of graph of the exponential function y = a® when 
O0<a<1. Herea=0.5. 
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Definition 39 growth- and decay curves 


Let a > 1. Then the graph of 
y=a" 


is called a growth curve, while that of 


x 


yra 


is called a decay curve. 
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Example 28 growth functions 


There are basically three laws of growth, namely unlimited, limited and lo- 
gistic growth, all of which involve an exponential function. The model is for 
unlimited growth, 


for limited growth, 


and for logistic growth, 
m 


t => 
y( ) 1+ qe—rmt 
where a, m and r are constants. 
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Example 29 interest compounding 


The value of a principal p compounded annually at an interest rate 7 for t 
years is, 
s=p(1+i)’ 


where 7 is expressed in decimal points. For compounding m times a year, 


then, 
j mt 
S=p (1 + ) 
m 


If the compounding is continuous, at 100 per cent interest for one year, then, 


1 m 
s=p lim (1+ =) = pe 


m—0oo m 


where e is the Euler’s constant, e = 2.71828... 
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Example 30 multiple compounding 


For multiple compounding, 


m 


i mt 
pil+ie)' =p (1 + =) 
the effective annual rate of interest is, 
i m 
= (1 de =) =f 
m 
The effective annual rate of interest for continuous compounding is, 


ie=e’—-1 
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Definition 40 discounting 


Discounting is the process of finding the present value p of a future sum of 
money Ss. 
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Example 31 discounting 
Discounting when under annual compounding is, 
s=p(1+i)® 


when under multiple compounding, 


Coe 
p=s|(1l+— 
m 


and when under continuous compounding, 


p=se™ 
When discounting, the interest rate 7 is called the rate of discount. 
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Example 32 converting exponential- to natural exponential functions 


A discrete growth s = p(1+i/m)™ can be converted to a continuous growth 


s = pe” thus, 
i mt 
s (1 ca =) 7 Pe 
m 


. mt 
In (1 + “ =Ine”™ 
m 


Therefore, 
eel In(1+4)t 
= 1 + min = min wl 
S=Dp ( =) pe 
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Example 33 reflections of graphs 

Reversing the sign of x, that is replacing x by —2, has the effect of reflection 
of the original graph with respect to the y-axis. Reversing the sign of y, that 
is replacing y by —y, gives a reflection of the same with respect to the x-axis. 
The graphs of y = a** remain always above the z-axis, in other words the 
function y = a+” maps —oco < x < oo toy > 0. The two functions y = a® 
and y = a~” are the reflection of each other with respect to the y-axis. It can 
be easily seen that the functions y = —a** are the reflection with respect to 
the x-axis respectively of y = a*®”. 
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Definition 41 logarithmic function 


The logarithmic function with base a is defined to be the inverse of the expo- 
nential function, and is written 


y = log, & 


where a > 0 anda 4 1. The logarithmic function of base 10 is called the 
common logarithmic function, and one of base e, where 


1 nm 
e= lim (1 + *) 
noo nr 


is called the natural logarithmic function. By the notation y = log, « we mean 
that the logarithm base a of x is the power to which a must be raised to get 
x. 
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Example 34 logarithmic function 


The domain of the logarithmic function 


y = log, x 


is the set of all positive real numbers, its range the set of all real numbers. 


The function is concave and increasing for a > 1, and is convex and decreasing 
for 0 <a< 1. Note also that log, x is the power which a must be raised to 
get x. 
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Example 35 examples of natural logarithm 
Note that 


e™* = qg=Ine% 


where a > 0, 


e?® — 7 =]ne” 


where x > 0, and 
elt fz) — f(x) = In el) 


where f(x) > 0. 
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Theorem 19 rules of logarithm 


Four basic rules for logarithm function are listed in the following. 
log, m + log, n = log, mn 
log, m — log, n = log, es 
n 


log, m* = zlog,m 
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Definition 42 elasticity of substitution 


The elasticity of substitution o is defined as, 


where f is called the least-cost input ratio, and ae the input-price ratio. 
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Example 36 values of the elasticity of substitution 
The value o = 0 means there is no substitutability, that is the two inputs 
are complements of each other and both must be used together in a fixed 
proportion. 


The value o = co means that the two goods may substitute each other per- 
fectly. Ultimately, 0 <0 < oo. 
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Definition 43 constant elasticity of substitution production function 


A constant elasticity of substitution production function is a production func- 
tion where, unlike the Cobb-Douglas function, has an elasticity of substitution 
whose value is constant but not necessarily 1. In its typical form, it is, 


lr 


q=a(ak* + (1—a)l-*) 


where a is called the efficiency parameter, a the distribution parameter, 3 the 
substitution parameter. Furthermore, 8 determines o, anda >0,0<a< 1, 
and 6 > —-1. 
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Example 37 logarithmic transformation of nonlinear functions 


Some nonlinear functions can be converted to linear functions using logarith- 
mic transformation, for example the Cobb-Douglas production function, 


gq = ak*l? 


which becomes 
Ing=Ina+alnk+6lnl 
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Other nonlinear functions can not be converted, for example the constant 
elasticity of substitution production function, 


lr 


q=a [ak +(1- al? ] 
which becomes just another nonlinear function, 
: -6 -p 
Ing=Ina— gm lak + (1—a)l~"} 
others 
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Example 38 nonlinear total revenue 
Let the total revenue be 
Tt = pq 
and the demand function 
p=a—bq 
where q is the quantity sold. Then r; expressed as a function of q is nonlinear, 
for 
ry = (a — bq)q = ag — bq? 
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Example 39 nonlinear total cost 
A more realistic equation for the total cost instead of 

Cc, = at bq 
is the nonlinear function 
c = aq? — bg? +ceq+d 
in which the production cost increases with quantity in at a decreasing rate 
(c, < 0) up to the inflection point at 
b 


~ 6a 
after which it increases at an increasing rate (c/ > 0). 
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Definition 44 polynomial 


A polynomial is an expression in the form 


n 
. ayv" 
i=0 


Here n is called the order of the polynomial. If n = 2 the polynomial is known 
as a quadratic polynomial, if n = 3 a cubic polynomial, if n = 4 a quartic, 
n = 5 a quintic and n = 6 a seztic. If we let p(x) be a polynomial, then 
a polynomial equation is the equation p(x) = 0. A polynomial function is a 
function of the form 


y = f(x) = p(2) 
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Example 40 quadratic equation 
The quadratic equation ax” + bx +c =0 has the solutions, 
—b+ Vb? — 4ac 


Oe (14) 


These solutions are called the roots of the quadratic equation. Equation 14 
is called the ‘minus-b formula’. The values of x obtained from the minus-b 
formula give the intersections of the graph of the quadratic function 


f(x) = p(x) = ax? + br +c 
on the z-axis. 
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The value b* — 4ac determines how the graph of f(z) lies relative to the x-axis, 
that is, 


> 0, there are two x-intersections 
b? —4dac’ =0, the graph touches the x-axis at one point 
<0, the graph never touches the x-axis 


Furthermore, the graph reverses its direction with respect to the y-axis at the 
critical point where f’(a7) = 0, that is when 


b 


~ 2a 
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Consequently the critical point is 


bb fy 
2a’ 4a * 
The graph of f(x) is symmetric with respect to the vertical line which passes 


through the turning point, that is to say, the line 


b 


~ 2a 


The y-intercept is at the point (0, c). 
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Example 41 hyperbolic function 
A hyperbolic relation has the form, 
(px + q)(ry+s)=t 


From this we obtain, 
S t 
(+2) (y+2)=+ 
Pp r pr 


_t 1 Ss 
oor r+4 r 
a 


—c 
x+b 
where p, g, Tr, s and t are constants, hence so are a = 


t 
pr? 
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In economics we sometimes find hyperbolic functions of the form, 


4 
~ br +e 


y (15) 
For example, a demand function of a good may be given by, 


m 
qta=— 
Pp 


which leads to, 
m 


~ gta 


Pp 


where p and q are respectively price and quantity demanded of a good, while 
m and a are constants. 
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The graph of Equation 15 has the z-axis, that is the line y = 0, as its hori- 
zontal asymptote, and has the line 


as its vertical asymptote. If all the parameters are positive, then the curve in 
the first quadrant decreases with a decreasing rate. 
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Definition 45 matrix addition 
Let A= {aij}, B = {b;;} and C = {c;;} be three matrices. Then 


C=A+B 
is called the addition of the matrices A and B if 
Cig = Gig + bi 


for all ¢ and j. 
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Definition 46 matrix multiplication 


Let A = (a;;) be an m x n matrix and B = (bg:) an n x p matrix. Then the 
product AB is an m x p matrix C = (cj) where, 


n 
C= S Gikdgl 
k=1 


where 1 <i<mand1<l<p. 
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Definition 47 determinant 


The expression obtained by eliminating the n variables x1,...,%, from n 
equations, 
441%, +... + A1nln = 0 


(16) 
AniX1 +... + Anntn = 0 


is called the determinant of this system of equations, Equation 16. The de- 
terminant of matrix A denoted by various different notations, for example 
det(A), |Al, Yo (ta1b2c3 “8 ‘), D(ay,b2c3 ve ‘), or |ayboc3 pace |. 
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Example 42 determinant of systems of three variables 


For a linear system of three variables, Equation 16 can be written as, 


a,x + aay +a3z =0 
bx + boy Tr b3z = 0 (17) 


12+ coy + ¢3z = 0 


Eliminating x, y and z from Equation 17 gives us, 


aybec3 a b3c2 t a3b1 C2 a2b1c3 + agb3c, = azboc1 =0 
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Definition 48 minor 


A minor Mj; of any matrix A is the determinant of a reduced matrix obtained 
by omitting the i*® row and the j*® column of A. 
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Theorem 20 Laplacian expansion 


Determinant can be determined by, 


k 
|A| = So aijsCi; 
i=l 


where Ci; is called the cofactor of a;;. The cofactor Ci; can also be denoted 
as a"), and, -_ 
Ci = (-1)'Y Mig 


where M;; is a minor of A. 
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Definition 49 permutation inversion 


Any pairwisely ordered pair in a permutation p is called a permutation inver- 
sion in pift > 7 and p; < pj. 
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Theorem 21 determination of determinant by permutation 


Determination of the determinant can also be determined by, 
|A| = yop) II Oy (4) 
wT w=1 


where z is a permutation which ranges over all permutations of {1,...,n}, 
and I(7) is called the inversion number of 7. 
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Theorem 22 properties of determinant 


Let a be a constant and A an n x n matrix. Then, 


|aA|=a"| Al 
[alee Al 
| AB| =| A||B| 
|Z] =|AA™*|=|A||A“*|=1 
1 
A — 
| | | A-"| 
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Definition 50 multilinearity 


A function in two or more variables is said to be multilinear if it is linear in 
each variable separately. 
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Theorem 23 multilinearity of determinants 


Determinants of matrix are multilinear in rows and columns. 
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Example 43 multilinearity of determinants 
Consider an 3 x 3 matrix, 
ay aq a3 
A= a4 a5 a6 
a7 ag ag 
What Theorem 23 says about multilinearity of determinants is the same as 
saying that, 


ay 0 0 0 ag 0 0 0 a3 
|A| =] a, as ag|]+]as, ads ag]+|as4 as ag 
a7 ag ag a7 ag ag a7 ag ag 
and 
a, ag a3 0 a2 a3 0 ag a3 
| A | =| 0 as ag}|+|a4 a5 ag|/+] 0 as ae 
0 ag ag 0 ag ag a7 ag ag 
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Definition 51 conformal mapping 


A conformal mapping is a transformation that preserves local angle. The 
terms function, map and transformation are synonyms. 


Business mathematics, Matrix, 15” November 2005 —13-From 5” Nov 2005 , as of 
26%" April, 2006 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 52 similarity transformation 


A similarity transformation is a conformal mapping the transformation matrix 
of which is, 
A’ = BAB" 


Here A and A’ are similar matrices. 
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Theorem 24 similarity transformation and determinant 


Similarity transformation does not change the determinant. 


Proof. The proof for this is simply, 


1 
ee en a hee Alay 1 
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Example 44 similarity transformation 


|B-1AB—\I| =|B-!AB— B-!)IB| 


=|B-1(A—dDB| 
=|B||A-AI||B| 
=|A-Al| 


Business mathematics, Matrix, 15” November 2005 —16-From 5” Nov 2005 , as of 
26" April, 2006 


164 26" April, 2006 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 53 matrix trace 


Let A be a square, n x n matrix. Then the trace of A is, 


Tr(A) = s Ay 
i=1 
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Definition 54 matrix transpose 
The transpose of a matrix 
A= {aij} 
is 


A? = {ays} 


Business mathematics, Matrix, 15'” November 2005 —18-From 5” Nov 2005 , as of 
26" April, 2006 


God’s Ayudhya’s Defence 26" April, 2006 165 


Business Mathematics, notes and projections kit Tyabandha, PhD 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 55 complex conjugate 


The complex conjugate of a matrix 


is a 
A= {aij} 
wherea=p-—qiifa=pt+di. 
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Definition 56 big-O notaton 
Let $(n) or g(x) be a positive function, and let f(n) or f(x) be any function. 
Then f = O(¢@) if |f| < Ad for some constant A and all values of n and 


x. Here O is called the big-O notation which denotes asymptoticity. The 
notation f = O(¢) is read, ‘f is of order @’. 
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Theorem 25 properties of determinant 
Some other properties of the determinant are, 
|A|=|A*| 


|A|=[A| 
| +eA| =1+Tr(A) + O(e?), for € small 
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Example 45 notes on determinants 


For a square matrix A, 

a. switching rows changes the sign of the determinant 

b. factoring out scalars from rows and columns leaves the value of the 
determinant unchanged 

c. adding rows and columns together leaves the determinant’s value 
unchanged 

d. multiplying a row by a constant c gives the same determinant mul- 
tiplied by c 

e. if a row or a column is zero, then the determinant is zero 

f. if any two rows or columns are equal, then the determinant is zero 
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Theorem 26 matriz trace 


Some properties of matrix trace are, 
Tr(A) = Tr (A*) 


Tr(A + B) = Tr(A) + Tr(B) 
Tr(@A) = aTr(A) 


Business mathematics, Matrix, 15” November 2005 —23-From 5” Nov 2005 , as of 
26%" April, 2006 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Problem 23 matriz transpose 
Prove that, 
(A?) = (4"1)" 
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Theorem 27 property of matrix transpose 


(AB)? = BT AT 
Proof. —— a . 
(B A he = (b Ve (a i 
= beiazK 
= a;k dK: = (AB): a (AB); 
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Definition 57 matrix inverse 
Let A be a square matrix. Then the inverse of A, if it exists, is 47! such 
that, 
AA + =I 
Furthermore, A is said to be nonsingular or invertible if its inverse exists, 
otherwise it is said to be singular. 
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Example 46 matrix inverse 


For a 2 x 2 matrix, 


the inverse of A is, 
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If A is a 3 x 3 matrix, then the inverse of A is, 


At = HF {det (mss)} 


where m,; is a minor of A. 

If A is an n x n matrix, then A~! can be found by numerical methods, for 
example Gauss-Jordan elimination, Gaussian elimination, and LU decompo- 
sition. 
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Example 47 finding matrix inverse 


The Gaussian elimination procedure solves the matrix equation Ax = b by 
first forming an augmented matrix equation [Ab] and then transform this 
into an upper triangular matrix [{ai,} b’ | , where aj, are all zero except when 
i<j. Then, 


k 
aed 1 b! / 
ran Aig U5 
ay. 9 
Ww j=itl 
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The Gauss-Jordan elimination procedure finds matrix inverse by first forming 
a matrix [A J], and then use the Gaussian elimination to transform this matrix 
into [I B]. The result matrix B is in fact A7?. 
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The LU decomposition forms from the matrix A a product LU of two matrices, 
one lower- while the other upper triangular. This gives us three types of 
equation to solve, namely when 7 < j,i = j andi > j, where i and 7 are the 
indices of row and respectively column of the matrix product. Then, 


Ax = (LU)x = L(Ux) =b 
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Letting y = Ux we have Ly = b, and therefore, 
by 


~ ha 


ea 
Yy 
w= 7 bi — Do its 


Y1 


where 7 = 2,...,n. 
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Then since Ux = y, 


Yn 
Ln = — 
Unn 
n 
1 
Y= [| Vi ; UjjL5 


j=i+l 


wherei=n-—1,..., 
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Theorem 28 matrix inverse 


Let A and B be two square matrices of equal size. Then, 
(AB) = BoA 
Proof. Let C = AB. Then B = A~!C and A =CB™", therefore, 
C= AB=(CB") (AC) =CB'A™C 


Hence CB~'A7! = I, and thus B~'A~' = (AB)7!. q 
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Definition 58 Einstein’s summation 


The Einstein’s summation is the simplification of notation by omitting a sum- 
mation sign, keeping in mind that repeated indices are implicitly summed 
over, for example }>, a;,a;; becomes 


Ajkaij 
and )*, a;a; becomes 


aya; 
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Definition 59 matrix multiplication 


The multiplication of two matrices A = {a;;} and B = {b;;} is the matrix 
C = AB such that 
Cik = Aizbdjr 
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Theorem 29 associativity of matrix multiplication 


The matrix multiplication is associative. 
Proof. 
[(ad)c],; = (ab) inch g = (@idie) Ces 


= ait (binCej) = ait(be)1z = [a(be)],; 
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Example 48 matric multiplication 


From Theorem 29, which shows us the associativity of matrix multiplication, 
we could write the multiplication of three matrices as [abc];;, which is the same 
as writing ajbi.cx%;. And this applies in a similar manner to the multiplication 
of four or more matrices. 
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Theorem 30 non-commutativity of matriz multiplication 


If A and B are two square and diagonal matrices, then 
AB=BA 


But in general matrix multiplication is not commutative. 


Business mathematics, Matrix, 15” November 2005 -39-From 5” Nov 2005 , as of 
26%" April, 2006 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 60 block matrix 


A block matrix is a matrix which is is made up of small matrices put together, 


for example, 
A B 
C D 


where A, B, C and D are matrices. 
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Theorem 31 block matrix multiplication 


Block matrices may be multiplied together in the usual manner, for example, 


Aj By Apo Bg = A, Ag + Bi Co A, Bo iF By, De 
Cy D, C2 Dz - COC Ao + Di Cz C, Bo aa D, Dz 


provided that all the products involved are possible. 
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Definition 61 diagonal matrix 


Let A = {aj;} be an n x n matrix. Then A is called a diagonal matrix if 
aij = 0 when i 4 j. Here 1 < i,j <n. In other words, a diagonal matrix 
has its components in the form aj; = cjd;;, where c; is a constant and 4;; is 
the Kronecker delta, 
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Theorem 32 matrix diagonalisation 


A square matrix A can be diagonalised by the transformation 
AS PDR: 


where P is made up of the eigenvectors of A and D is the diagonal matrix 
desired. 
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Example 49 matrix diagonalisation 


Matrix diagonalisation can greatly help reducing the number of parameters in 
a system of equations. For instance, the systems Ax = y when diagonalised 
becomes 

PDP-!x=y 


that is Dx’ = y’, where x’ = P~'x and y’ = P~'y. In this case, if A is an 
n Xn matrix, we say that our new system obtained through the process of 
diagonalisation has canonicalised from n x n to n parameters. 
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Definition 62 symmetric matrix 
A symmetric matrix is a square matrix A which satisfies 
AT=A 
Example 50 symmetric matrix 
If A is a symmetric matrix, then 
AMATHT 
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Definition 63 orthogonal matrix 
Let A be a square matrix. Then A is said to be orthogonal if 
AAT =I 
Example 51 orthogonal matrix 
Definition 63 is the same as saying that 
Al =A? 
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Theorem 33 symmetric matriz 


A matrix A is symmetric if it can be expressed as 
A= QDQ? 


where Q is an orthogonal matrix and D is a diagonal matrix. 
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Example 52 symmetric matrix 


Any square matrix A may be decomposed into two terms added together, 
that is A, + Aq where A, is a symmetric matrix and A, an antisymmetric 
matrix, called respectively a symmetric part and an antisymmetric part of A. 
Furthermore, 


As = = (A+ A’) 


Nl re 


and, 


Aq = = (A- AT) 


Nl eRe 
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Definition 64 inverse matrix 


Let A be a square, nonsingular matrix. Then the inverse matrix A~' of A is 
a unique matrix for which, 


AAS =f =A-1A 
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Example 53 finding the inverse of a matrix 


An inverse matrix may be found using the formula, 


A = ra Adi 
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Example 54 solving linear equations with the inverse 


Matrix equations of the form 
Ax=b 


can be solved with the help of the inverse matrix A~! as 
x=A7'b 
where A is an n X n matrix, x a vector of size n whose components are 


variables, and b a vector of size n containing constants. 


Business mathematics, Linear algebra, 22" November 2005 —3- From 5°” 
November 2005 , as of 26'" April, 2006 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Theorem 34 Cramer’s rule 


Let A be the coefficient matrix and A; a matrix formed from A by replacing 
the column of coefficients of x; with the column vector of constants. Cramer’s 
rule solves a system of linear equations through the use of determinants as 
follows. 


oe, = il 
i= 
| A| 
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Definition 65 Jacobian determinant 


Let a system of n functions not necessarily linear be 


Y1 =fi (Dice rane) 


Yn Sf (Biya Ord) 
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Then a Jacobian determinant comprises all the first-order partial derivatives 
of the system arranged in ordered sequence, that is 


Oy. ON 
Ox, OLn 
| OYyts Oyun | : 
| J | ~~ | Ox1,.:.;00% | , 
dyn Oyn 
0x1 OLn 
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Theorem 35 functional dependence from Jacobian determinant 
Let a system of n equations be 
Y= fi (ip ices ey) 


7=1,...,n. 
If | J| = 0, then y; are functionally dependent. 
On the other hand if | J| 4 0, then y; are functionally independent. 
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Definition 66 Hessian 


A determinant | H | composed of all the second-order partial derivatives, with 
the direct partials on the principal diagonal and the cross partials off the 
same, is called a Hessian. In other words, let a multivariable function be 


z= f(x,y) 
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Then the Hessian of z is 
zyx Zyy 


where Zzy = Zyz. Moreover, the first principal minor is 
| veut | = £u2u 
and the second principal minor is 


Zz z 
| H, | _ LL cy | _ Cian ns Care 


Zey yy 
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Theorem 36 optimality of a multivariable function 
Let a multivariable function be 
z= f(x,y) 
and let the first-order conditions 


Zn = Zy =0 


are met. Then a sufficient condition for z to be at optimum is 


a (Beale 


together with 
ZrasZyy <0 
in case of a maximum and 
Zan, Zyy > 0 
in case of a minimum. 
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Definition 67 definiteness 
From Definition 66, 


if |H,| > 0 and | H2| > 0 the Hessian | H | is said to be positive definite, and 
the second-order conditions for the minimum are met. 


If | H,| < 0 and | H2| > 0 it is said to be negative definite, and the second- 
order conditions for the maximum are met. 
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Algorithm 1 Procedure to test for the optimality of multivariable functions of 
two variables. 


z= f(x,y) 
find zz and zy 
if z, =0 and z, = 0 then 
find Zz2, Zry and Zyy 
find H; and Ho 
if |H,| >0 and |H.|>0 then 
| H| is positive definite 
elseif | H; | < 0 and | H2 > 0| then 
| H| is negative definite 
endif 
endif 
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Definition 68 general Hessian 


Let y = f(x1,-..,%n) be function of n variables. Then the n**-order Hessian 
for this function is 

Yl oct" Yin 

[al : 


Ynl1 aie Ynn 


Then the first principal minor | H, | is simply x11, and the i*” principal minor 
is 


Yr Yui 
| Hi| = 
Yr oc Vii 
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Theorem 37 positive definiteness through Hessian 


Let y = f(a1,...,%n) be function of n variables. Let the Hessian of y be 
represented by | H |. 

Then if all the principal minors of | H| are positive, then | H| is positive 
definite and the second-order conditions for a relative minimum are met. 

If the sign of the principal minors alternates between negagive and positive, 
then | H | is negative definite and the second-order conditions for a relative 
maximum are met. 
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Example 55 positive definiteness for three dimensions 
For 
y = f (x1, 22,23) 


the third-order Hessian is 


Yur Y12 13 
|H|=|yo1 yee Yyo3 

Y¥31  -Y32 «33 

where 9 7 
y ¥ 
=.=, = ———,, and so on 
te One ate 022021 
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The first-, second- and third-order Hessian’s are respectively 


|i | = yu, |H| = [Bs ik 


Ya. Y22 
and 
Yr Y12 13 
| H3 | = Y21 Y22 Y23 
Y31 ~Y32 = -Y33 
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If |H,| > 0, | H2| > 0 and | H3| > 0, then H is positive definite and the 
second-order condition for minimum is fulfilled. 


If | H,| < 0, | H2| > 0 and | H3| < 0, then | H| is negative definite and the 
second-order condition for maximum is satisfied. 
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Definition 69 discriminant 


A discriminant is a determinant of a quadratic form. Let the quadratic form 
be 
edt 2 
z=ax* + bry +cy 


which is in matrix form 


a 2) [¢ 
ral alt lls] 
Then the discriminant is " 
ie 
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The first principal minor of the discriminant is 


| Dy, | =a 
and the second principal minor 
b 2 
a 3 b 
D — 2 =ac- — 
| Po | B C 4 
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Theorem 38 definiteness of a function by the discriminant 


Let a quadratic form be 
z= ax" + bay + cy? 


and let the discriminant of z be | D|. 


If |D,| > 0 and | D2| > 0, then | D| is positive definite and z > 0 for all 
Ly # 0. 


If |D,| < 0 and | D2| > 0, then | D| is negative definite and z < 0 for all 
£,y # 0. 
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Theorem 39 constrained optimisation with Lagrange multipliers 
Let 
f(xy) 


be a function subject to a constraint 


g(x,y) =k 


where & is a constant. Then the optimisation of f can be done by first 
transforming f together with g into a new function 


F(z,y, X) = f(x,y) + A(k— g(x,y) 
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and then solve the following equations, 


F,(x,y, A) =0 
Fy(x,y, A) =0 
Fy (x,y, A) =0 


to obtain the critical values xo, yo and Ao at which F and hence f are opti- 
mised. 
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Definition 70 constrained optimisation 


In the constrained optimisation with Lagrange multipliers in Theorem 39 


above, f is called an 
objective or origin function 
and F the 


Lagrangian function 


Business mathematics, Linear algebra, 22"¢ November 2005 
November 2005 , as of 26" April, 2006 


~23— From 5"” 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 71 bordered Hessian 
Let 
P(Bigsivs LR) 


be a function of n variables subject to constraints 


g(@1,---,;%n) 


Let 


F(a1,..-,@n,A) = f(a,...,%n) +A(k- g(a1,... 
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Then the bordered Hessian | H | is defined as either 


Fu Fig ss) Fin 1 
Foi 92 
|H|=| : : : 
Eni Fan Gn 
J. g2 gn 0 
or 
0 91 Gn 
ie gn Fu Fin 
|H| = 8 
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This is simply the Hessian 


Fics Fin 
Fut Fi Finn 
bordered by the first derivatives of the constraint with zero on the principal 


diagonal. The order of a bordered principal minor being determined by the 
order of the principal minor being bordered, 


|H|=|H,| 


since in this case an n X n principal minor is being bordered. 
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Theorem 40 bordered Hessian 


Let f(a1,...,%n) be a function of n variables subject to constraints 


g(@1,---,;%n) 


Let | H| be the bordered Hessian defined in Definition 71. 


Then if | H2|,...,|Hn| <0, then the bordered Hessian | H | is positive defi- 
nite, and therefore is a sufficient condition for a minimum. 


If | H2| > 0| H3| <0, | H4| > 0, and so alternatingly on, then | H | is negative 


definite, which is a sufficient condition for a maximum. 
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Example 56 constrained optimisation using bordered Hessian 


Let f(x,y) be a function to be optimised subject to a constraint g(x,y) =k, 
where k is a constant. Then the Lagrangian function becomes 


F(a,y, X) = f(x,y) + A(k— g(x,y) 


The first-order conditions for optimisation are 


Business mathematics, Linear algebra, 22" November 2005 —28-— From 5°” 
November 2005 , as of 26°” April, 2006 


194 26" April, 2006 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


The second-order conditions for optimisation can be expressed together as a 
bordered Hessian 


= Pox Foy Gx 
|H| =| Fyc Fyy Gy 
Ix Dy 0 
or 
_ 0 Gx Dy 
| A | =]9e Pra Frey 
Iy Pye Fyy 
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Note 1 first- and second-order conditions for optimisation of a function 
subject to some constraints 
Theorem 39 gives the first-order conditions for optimising a function sub- 
ject to some constraints. Theorem 40 gives the second-order conditions for 
optimising a function subject to some constraints. 
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Definition 72 Marshallian demand function 


A Marshallian demand function gives an expression of the amount of a good 
that a consumer will buy as a function of commodity prices and income avail- 
able. It is derived by maximising the utility subjected to a budgetary con- 
straint. 
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Example 57 derivation of a Marshallian demand function 


Let a utility be 
U= 41492 


which is subject to a constraint 


pig + page = b 


where b is the amount of income available, that is to say, our budget. Then 
the Lagrangian function is 


U = qigg + A(b— pigi — P2492) 
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The first partial derivatives are then 


U1 = q2—- Api =0 (18) 
U2 = q1— Ape =0 (19) 
uy =b—piqi — pag = 0 (20) 


where u;, ug are respectively ug, and ug. 
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Simultaneously solving Equation’s 18, 19 and 20 leads us to 


Boge ae 
PL p2 


Hence g2 = qipi/p2 and q = q2p2/pi and from Equation 20 we have, 


Pd p2q2 
b= pig, + po— = pogo + pir — 
P2 Pl 


which yield us, for qj and gg, the Marshallian demand functions which max- 
imise satisfaction of the consumer subject to income and prices. 


Business mathematics, Linear algebra, 22" November 2005 —34— From 5°” 
November 2005 , as of 26°” April, 2006 


God’s Ayudhya’s Defence 26" April, 2006 197 


Business Mathematics, notes and projections kit Tyabandha, PhD 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Next, we test the second-order conditions by firstly finding uy, = 0, u22 = 0, 
uy2 = U21 = 1, gi = pi and go = po, which give us 


0 1 pi 
|H|=/1 0 pe 
pi po O 


which gives | H2| = 2p1p2 > 0 Hence | H| is negative definite and thus u is 


maximised. 
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Definition 73 input-output analysis 


The production process of producing one good usually requires the input of 
many other intermediate goods. Let x; be the total demand for product i, 
and let b be the final demand for the product from the ultimate users. Then, 


Uj, = A141 +... + Ajinkn + b; 
fori =1,...,n, where a;; is a technical coefficient which represents the value 


of input 7 required to produce one monetary unit’s worth of product j. 
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If we consider the total demand for every one of the products, then 


x= Ax+b 
where 
T1 G11 *°*) Gin by 
Ln Ani *** Gann bn 
It follows from this that 
x=(I-—A)-'b 
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The matrix A is known as the matrix of technical coefficients. It is also 
known as the input-output table, the rows being the inputs and the columns 
the outputs.. The matrix J — A is known as the Leontief matriz. 
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Example 58 complete input-output table 


In a complete input-output table, labour and capital would also be included 
as inputs. These give the value added by the firm. They are normally put 
as an extra row at the bottom of the matrix of technical coefficients A. The 
vertical summation of each column of the table is then equal to 1. 
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Definition 74 eigenvalue and eigenvector 


Let A be a square matrix. Then a scalar A such that the equation 
Av =v (21) 


holds for some vector v ¥ 0 is called an eigenvalue ¢ of A, and the vector v is 
called an eigenvector of A corresponding to the eigenvalue A. The eigenvalue A 
is also known as the characteristic root, or the latent root, while the eigenvector 
is also known as the characteristic vector, or the latent vector. 
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Note 2 eigenvalue 


From Equation 21 it follows directly that 
(A—Al)v =0 (22) 


Then A—ALI is called the characteristic matrix of A. Since v assumes a unique 
value and by assumption v 4 0, it follows that A— AJ must be singular, which 
means that its rows must be a multiple of one another. Now A — AI is zero 
if and only if the characteristic determinant | A — AI | of A is zero. 
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In other words 
|A—AI| =0 (23) 


which is called the characteristic equation of A. With Equation 23 there will 
be an infinite solution for v in Equation 22. In particular, if v is a solution, 
that is if it is an eigenvector, so is kv for any k # 0. We force a unique 
solution by using the normalisation 


We =a 
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Then the sign-definiteness of A can be determined from the characteristic 
roots A’s. 


Thus if all \’s are positive, then A is positive definite; and if negative, negative 
definite. 

Let at least one » be zero, which is neither positive nor negative, if all the 
remaining \’s are nonnegative, then A is positive semidefinite; and if they are 
nonpositive, negative semidefinite. 

Lastly, if some of the \’s are positive while others are negative, then A is 
indefinite. 
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Note 3 
We have seen in Note 2 how, having found ;, where 7 = 1,...,n, we find 


through normalisation the corresponding, unique v;. On the other hand if 
we have found first the v;’s, their corresponding ;’s may be found by first 
forming a transformation matria 


T=[v1 .-- Vn] 


and then the corresponding eigenvalues or the characteristic roots are ob- 
tained from 


dy De 286 1H 
A2 
TTA S 
0 a 
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Definition 75 


The vector equation, Equation 21, has as its solutions the zero vector v = 0 
together with all the corresponding eigenvalue-eigenvector pairs. The set of 
all the eigenvalues of A is called the spectrum of A. The spectral radius of A 
is then the largest of all the absolute values of the eigenvalues of A, that is to 
say, 

max |X; 


The set of all eigenvectors v;;, together with 0, forms a vector space called 


the eigenspace of A corresponding to j. 
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Definition 76 optimisation 


A problem of optimisation is one in which one tries to maximise or minimise 
a certain quantity called the objective, which depends on a finite number of 
variables. These may be either independent or related to one another through 
some constraints. 
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Definition 77 mathematical programme 


A mathematical programme is an optimisation problem in which the objective 
and the constraints are given as functions or mathematical relationship. In 
other words, 


optimise: z= f(x1,...,2n 


subject to: g;(a1,...,2n) 


Some constraints are explicitly stated as requirements, others are hidden con- 
ditions. These latter need to be pin-pointed through the study and under- 
standing of the model and its inputs. 
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Definition 78 linear programme 


A linear programme is a mathematical programme all the functions involved 
of which are linear. This means that, 


f(@1,---,0n) = C11 +++ + Cn en 
Gi(@1,---,Un) = a2, + +++ + Gintn 
where 7 = 1,...,m and c; and aj;, 7 = 1,...,n, are constants. If there is 


an additional restriction on the input variables that they be all integers, then 
the optimisation problem is called an integer programme. A mathematical 
programme which is not a linear programme is said to be nonlinear. 
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Definition 79 quadratic programme 


A quadratic programme is a mathematical programme in which all the con- 
straints are linear and the objective function is in quadratic form, which is in 


general, 
n n n 
f(x, are | = S- S- Cig Lily + > dx; 
w=1 


i=1 j=1 


where c;; and d; are constants. 
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Definition 80 standard form 
A linear programme is said to be in standard form if all the constraints are 
equalities and if one feasible solution is known. In other words, our problem 
is now 
optimise: z= c’x 
subject to: Ax=b 


with: x>0 
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Definition 81 initial feasible solution 


One may change any linear programme into the standard form by adding a 
slack variable to the left-hand side of a constraint of the form })a,;;7; < b; 
to obtain 


n 
YF aijay + tp, = bi 
j=l 


where pp > n and k = 1,2,.... 
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Similarly one may add a surplus variable to the right-hand side of a constraint 
of the form > Ajj L5 = b; to obtain > Aijiti = bj + La 


n 
DT aijty — tq = bi 
j=l 


where q > n and /=1,2,.... Next, all the slack and surplus variables are 
added to the objective function with zero coefficients. 
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Then if we add an artificial variable to the left-hand side of each constraint 
where there is no slack variable, then the initial feasible solution is xo = b, 
where x is the vector of slack and artificial variables. The artificial variables 
are added to the objective function with a large negative coefficient —M. 
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Definition 82 linear dependence 


A set of n vectors of m dimensions {v1,...,Vn} is said to be linearly dependent 
if there exist some constants ay,...,@, not all of which are zero, such that 
Q1Vv, +--+ +anvn = 0 (24) 


It is said to be linearly independent if the condition in Equation 24 implies 
ay = ag =+:'=An = 0. 
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Theorem 41 


Consider a set of n vectors of m dimensions. If n > m, then the set is linearly 
dependent. 
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Definition 83 conver combination 


A vector v is called a convex combination of vectors v1,...,Vn if there exist 
some nonnegative constants 31,...,8,, where 


Brt--+By=1 


such that 
v= ivi +--+ Brvn 
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Definition 84 convex set 


A set of m-dimensional vectors is said to be convex if for any two vectors 
belonging to the set the line segment between them also belongs to the set. 
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Theorem 42 


All points on the line segment joining any two vectors may be expressed as a 
convex combination of the two vectors. 
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Definition 85 extreme point 


A vector v is called an extreme point of a convex set if it can not be expressed 
as a convex combination of two other vectors in the set. 
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In other words, an extreme point of a convex set K is a point x in K that 
cannot be written as x = y+ (1—0)z with 0 < 6 < 1, y and z in K, and 
y # z, that is to say, an extreme point is a point which is not an interior 
point of any line segment belonging to Kk. 
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An equivalent definition of an extreme point is that x is an extreme point of 
a convex set K if K\ {a} is convex. 
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Definition 86 bounded set 


A metric space is a non-empty set X for which is defined a concept of distance. 
The distance d is called a metric on X, having such properties that, for any 
points x and y in X, we have d(x, y) > 0, and d(x, y) = 0 implies 


r=y 
Furthermore, 

d(x, y) = dy, x) 
and 


U(x, y) < d(x, z) + d(z,y) 
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Let X be a metric space with metric d, let A be a subset of X and let x be 
any point of X. Then the distance from x to A is defined as 


d(x, A) = inf {d(z,a): ae A} 
whereas the diameter of A is defined as 


d(A) = sup {d(a1, a2) : a; and ag € A} 
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Then a set is said to be bounded if its diameter is finite. 
Further, let zp be a point in X and r a positive real number. Then the open 
sphere S;(%o9) with centre xp and radius r is the subset of X defined by 


S,(ao) = {@ : d(x, x0) <r} 
A point x in X is called a limit point of A if each open sphere centred on x 


contains at least one point of A different from a. A subset F of X is said to 
be closed if it contains all its limit points. 
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Definition 87 linear space 
A linear space, aka a vector space, is a non-empty set L on which is defined two 
binary processes, say addition and scalar multiplication. Addition is defined 
such that for any x, y and z in L, then «+ y is again in L; 


E+y=yrtsr 


at+(ytz)=(@+y)+2 


there exists a unique identity element 0, aka zero element or the origin, such 
that « +0 =z for every x; and there exists a unique inverse element —x for 
every x, such that a + (—x) = 0. 
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Scalar multiplication is defined with regard to scalars, some instances of which 
are real and complex numbers, such that for any scalar a and any x and y in 
L, ax is again in L 

a(a+y) =ar+ay 


(a+ B)x =axn4+ Bu 
(a8)x = a(Gx) 


and lx = x, where 1 is the identity for scalar multiplication. 


Business mathematics, Linear programming, 29°” November 2005 -21-— From 5‘” 
November 2005 , as of 26'" April, 2006 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


A normed linear space is a linear space on which is defined a norm, that is 
a function which maps each element x in the space to a real number ||2|| in 
such a manner that ||2|| > 0, and ||2|| = 0 if and only if « = 0 


Ila + yl < lel] + IlyIl 


and 
lax|| = |a |x|] 
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Theorem 43 


A normed linear space is a metric space. 


Theorem 44 


Any vector in a closed and bounded convex set with a finite number of extreme 
points can be expressed as a convex combinations of the extreme points. 
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Definition 88 definition of the problem 


For two vectors, that is points, x and y in R”, we write x > y if and only 
if 2; > y; for all 1 <i< n. A system of m weak linear inequalilties in n 
variables can be written as Ax > b, where A is an m x n matrix. 


A fundamental question concerning such system is whether it is consistent, 
that is to say, whether there exists some x such that Ax = b. 


A system may be inconsistent, or it may have a set of solutions which is 
unbounded. If we sketch our problem on a graph, we may see that it’s solution 
set is convex. 
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Theorem 45 solution space 


The solution space of a set of simultaneous linear equations is a convex set 
the number of extreme points of which is finite. 
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Theorem 46 extreme-point solution 


Let S be the set of all feasible solutions to the linear programme in standard 
form in Definition 80, in other words, S is the set of all vectors x that satisfy 
Ax = b and x > 0, where A is an m x n matrix. Then S is a convex set, and 
the number of its extreme points is finite. The objective function attains its 
optimum, provided that one exists, at an extreme point of S. If m <n, then 
the extreme points of S have at least n — m zero components. 
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Algorithm 2 
Procedure for finding basic feasible solutions.: 
Input: Ax = b, A is an m x n matrix, m <n, rank A =m 
[ar --: anJ<A 
(aja, +--+ +a4n;a, = b) + (Ax =b) 
for i=m+1 ton do 
endfor 
(@1,.-.,Un) <— solve xa; + +--+ 2,a, =b 
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Definition 89 simplex method 


The simplex method is a matrix procedure which solves linear programmes of 
the standard form as described in Definition 80 where b > 0. Starting from a 
basic feasible solution xg we locate successively other basic feasible solutions 
giving better values for our objective. For minimisation programmes the 
method uses Table 3, for maximisation programmes the same table is also 
used but with the sign of entries in the bottom row reversed. 
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Table 3 
Table used for minimisation programming in simplex method.: 
xP 
cl 
Xo Co A b 
ch cf A —cpb 
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Table 4 Description of the simplex method. 
while negative number exists in d do 


Locate the most negative number in the bottom row of the simplex 
table, excluding the last column. The column in which we find 
this number is called the work column. If more than one such 
column exist, choose one of them. 
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Find the smallest of the ratios between the elements in the last column 
and the elements in the work column of the same row, if these 
latter are positive. The element in the work column that yields 
this smallest ratio is called the pivot element. If there are more 
than one of these, choose one. If none of the elements in the work 
column is positive, the programme has no solution. 
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Using elementary row operations, convert the pivot element to 1 and 
reduce all other elements in the work column to 0. 

Replace the z-variable in the pivot row and first column by the z- 
variable in the first row and pivot column. This new first column 
then becomes the current set of basic variables. 


endwhile 
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The optimal solution is one in which all the basic variables assume the cor- 
responding values in the last column, while the remaining variables are zero. 
The optimal value of the objective function is then the value of the last row 
and last column for a maximisation programme, and the negative of this value 
if the programme is one of minimisation. 
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Algorithm 3 Algorithm for the simplex procedure. 
j<0 
while there exists a negative number in d do 
jejtl 
for i=1 ton do 
{c’s} < (column no. of the most negative no. in the bottom row) 
(work column) ¢ choose one of the {c}’s 
k © (work column) 
endfor 
Ppivot — M 
co 0 
soln < 0 
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for i=1 tom do 
if (a;),, > 0 then 
soln —1 


if p < Ppivor then 
rei 
endif 
endif 
endfor 
if soln = 0 then 
no solutions exist 
endif 
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convert} A, such that (a;),,, = 1 and (a;),, =0,1<i<m,i¥Fr 
(Xo), — Te 

endwhile 

for i=1 tom do 
(x0); — (by); 

endfor 

for i=m-+1 ton do 
xr; +0 

endfor 

i 

if the programme is one of minimisation then 
a ae 


endif 
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Definition 90 two-phase method 


The two-phase method is a procedure modified from the simplex method to 
cope with cases when artificial variables exist in the initial solution xo, in 
order to minimise the round-off errors that occur in the calculation. The last 
row in Table 3 in this case is 


d=c? —cfA=d,+ Md, 
and consequently we have Table 5 which is used here. Algorithm 3 is then 
firstly applied to the last row, and then again to those elements directly above 


the zeros in that row. 
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When an artificial variable is removed from the first column of the table, it 
ceases to be basic and may be removed from the top row of the table together 
with the entire column under it. When the last row contains only zeros, it 
may be deleted from the table. The programme has no solution if non-zero 
artificial variables are present in the final basic set. 
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Table 3 Table used for minimisation programming using the two-phase 
method. 


x? 
cl 

Xo Co A b 
d, —cjb 
dy 
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Definition 91 duality in linear programming 


Given a linear programme in the variables 21,...,2,, there exists another 
linear programme associated with it, called its dual, which is in the vari- 
ables wi,...,Wm-. The original programme is called the primal. The primal 
completely determines the form of its dual. 
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The symmetric dual of a primal linear programme in the matrix form 


minimise: z=c!x 
subject to: Ax >b 
with: x>0O 


is the linear programme 


maximise: z= b’w 
subject to: A’w<c 
with: w>0O 
The dual variables w1,...,Wm are known as shadow costs. 
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The unsymmetric dual of the primal 


minimise: z=c x 
subject to: Ax=b 
with: x>0 
is 
maximise: z= b!w 


subject to: A’w<c 
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The unsymmetric dual of the primal 


maximise: z=c x 
subject to: Ax=b 
with: x>0 
is 
minimise: z= b!w 


subject to: Aw >c 
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Note 4 


We may see from Definition 91 that the dual of a programme in standard 
form is not itself in standard form. These duals are said to be unsymmetric. 
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Theorem 47 


If an optimal solution exists for either the primal or the dual programme, 
then the other programme also has an optimal solution. If the duality is 
symmetric, then the two functions have the same optimal value. If the duality 
if unsymmetric, then the optimal value of each function can be derived from 
that of the other. 
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Definition 92 cut algorithms 
Algorithms which change the boundary of the solution region in order to find 
the optimal solution of an integer programme are called cut algorithms. 


The branch-and-bound algorithm does this by splitting the solution region 
into two and then discard the one which does not contain the optimal solution. 


The Gomory algorithm, on the other hand, reduce the feasible region with 
the help of a new constraint without the region being splitted. 
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Definition 93 branching 


We call branching a process by which a programme whose solution contains 
a non-integral 

JS UK k 
is made into two separate programmes having the additional constraint 


Li j 


in one, and 


in the other, the objective together with all the constraints of the original 
problem of which remain the same. Here 7 and k are positive integers and 
g<k. 
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Definition 94 bounding 


In the branch-and-bound algorithm, if the objective is maximisation, the value 
of the objective obtained when the first integral approximation occurs is said 
to be the lower bound for the problem, and if the objective is minimisation it 
is said to be the upper bound of the same. 
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Algorithm 4 Branch-and-bound algorithm for integer programming. 


find first approximation 
while approximations not all integers do 
choose x; from all non-integral variables such that 
min (|x; — |x;||, |"; — [av;]|) is maximised 
branch 
choose the branch whose value of the objective 
is maximum 
endwhile 
solutiontlast approximation 
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Example 59 branch-and-bound example 
(Problem 6.9; Bronson, 1982) 


maximise: z= 2%, +2%2+ 23 
subject to: 27, + 3272 + 3x73 < 11 


with: all variables non-negative and integral 


Solve by branch-and-bound algorithm. 
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Draw a simplex table of Programme 1. 


Ty no) v3 v4 
1 0 
L4 0) 2 3 3 1 11 
-1 -2 -1 0 0 
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Replace x4 for x2 as the basic variable. 


Ly v2 x3 L4 
m a a a 

1 2 22 

3 0 1 3 3 
=] =36,.0) S25 = 0h 
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Since 3 < #3 < 4, branch into two programmes, namely Programme 1 where 
v2 <3, and Programme 2 where x2 > 4. Consider first Programme 2. 


maximise: z= 2%, +2%2+ 23 
subject to: 2271; + 372 + 323 < 11 
rw <3 
with: all variables non-negative and integral 
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Use the simplex method in a tabulated form. 


Ty r2 x3 v4 x5 
1 2 1 0 0 
wa 0 2 3 3 1 0 11 
U5 0 0 1 0 0 1 3 
-1 -2 -1 0 0 0 
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Replace the basic variable x5 with x. 


Ly op) 23 LA U5 
v4 2 0 3 1 -3 2 
Lr 0 1 0 0 3 
-1 0 -l 0 2 6 
-10— From 5‘” 
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Replace the basic variable x4 with x1. 


Ly v2 x3 v4 x5 
vi 0 5 3 1 
xg 0 1 0 0 1 3 
0 0 3 3 3 7 
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Then consider Programme 3. 


maximise: z= 2, +2%2+ 273 
subject to: 2271; + 372 + 3x3 < 11 
v2 > 4 


with: all variables non-negative and integral 
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Draw a table for the two-phase method. 


@ 2 £3 @ 4X5 Xr 
0 0 —-M 
LA 0 2 3 3 1 0 0 11 
te —M 0 1 0 0 -1 1 4 
-1 -2 -1 0 0 0 0 
0 -1 0 0 1 -1 -15 
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Change x4 for x2 in the basic variables. 
T1 v2 v3 v4 U5 v6 
v2 3 1 1 3 0 0 + 
LE —3 0 aa 3 -1 1 $ 
1 2 22 
5 0 1 5 0 0 > 
2 1 34 
5 0 1 5 1 -1 => 
The coefficient parts of the row corresponding to the basic variable xg and 
the last row cancel each other. The optimal result is 23 = Ht, Li = 2% = 
vy = af = 2§ =0 and z* = 2%. 
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Therefore the solution is <j = 1, 73 = 3, x3 = 2, = x5 =0, and z* =7. 
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Algorithm 5 Gomory algorithm for integer programming. 


while solution not wholly all integers do 
choose one non-integral optimal approximation 
write a relation from the row where that variable is basic 
rewrite the relation to make all fractional coefficients 
some integer plus a proper fraction 
move all the fractions to LHS, and all the non-fractions 
to RHS 
write a new constraint as LHS > 0 
find the solution for the original problem together 
with the new constraint 
endwhile 
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Example 60 
(Problem 7.1; Bronson, 1982) 


maximise: z= 227,+ 22 
subject to: 2%; +529 < 17 
321 + 229 < 10 


with: 21,22 non-negative and integral 


Use cut algorithm. 
Solve 


Business mathematics, Integer programming, 6°” December 2005 -17—- From 5‘” 
November 2005 , as of 26" April, 2006 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Find the first approximation of Programme 1 normally using the simplex 
method. 


Ly x2 x3 LA 

0 0 
x3 0 2 5 1 0 17 
LA 0 3 2 0 1 10 
-2 -l1 0 0 0 
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Since 2» < zZ, we know that 3 is the pivot element, and therefore we replace 
the basic variable x4 with 21. 
Vy v2 x3 v4 
11 2 31 
3 0 3 1 -3 3 
m1 a es m 
1 2 20 
0 3 0 3 3 
We have aj = 2, 23 = 4, 23 = cj =O and 2 = 2. 
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Since both xj and 23 are non-integers, arbitrarily choose the former to gen- 
erate a new constraint. Then our Programme 2 becomes, 


| — =3 
Ly Dated 3 ak 
2 1 1 3 
to + ~2x — 
372 34 3 Ly 
2 1 D4 
Ss ab ihe Se 
en aie og 
2 Fe i 
ae air z 
aa care 
2% +44 >1 


and our new programme becomes 
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maximise: 


subject to: 


with: 


November 2005 , as of 26¢” 


Business mathematics, Integer programming, 6°” December 2005 


Department of Mathematics, Mahidol University 


z= 24, + 2724+ 0273 + 024 


11 2 31 
Bee ree a 

_2 1 10 
ey 372 34 = 3 
2% +44 >1 


all variables non-negative and integral 
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ry v2 23. v4 a5 rE 
2 1 0 0 0 —M 
a, 0 > 0 | O 0 y 
wz (OO 0 #2 1 -2 0 0 a 
te —M 0 2 0 dl -l 1 1 
-2 -1 0 0 0 0 0 
0 -2 0 -1 1 -1 -1 
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Now £2 replaces zg in the basic variables and becomes the pivot element. 


Ly x2 v3 LA x5 v 
v1 1 0 0 0 $ -3 3 
x3 0 0 1 % * % i 
xr 0 1 0 5 5 5 $ 
2 0 0 2 2 8 2 
0 0 0 0 0 0 0 
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This becomes, 


T1 r2 x3 LA x5 
vy 1 0 0 0 ; 3 
13 DS ur 
re 0 1 0 4 —4 4 
. ot © @ 4 8 
Then our first approximation of Programme 2 is xj = 3, x5 $, xe a 


ee See * _ 13 
ry = 25 =0, and 2*= >. 
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Arbitrarily choose x35 to generate the new constraint. 


1 1 1 
24 = 
2 T4 U5 2 
1 1 Lr 
5 a a ca aa 
1 1 1 1 
94 Chas a= 
t—-%5>1 
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Then our Programme 3 becomes, 


maximise: z= 22%; +%2+02r3 + 024+ Or5 


subject to: 21+ =%5 =3 


with: all variables non-negative and integral 
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Ly x2 x3 L4 25 LE “7 
2 1 0 0) 0 0 —-M 
a 600 1 0 0 0 ¢ O °O 3 
© 0 0 1 0 $ — 4 0 0 $ 
tz «(0 1 0 0 -3 # 0 0 a 
zz —-M 0 0 0 -1 -1 1 1 
-2. -1 0 0 0 0 0 0 
0 0 0 -1 1 1 -1 -1 
Business mathematics, Integer programming, 6°” December 2005 -27- From 5‘” 
November 2005 , as of 26" April, 2006 
Kit Tyabandha, PhD Department of Mathematics, Mahidol University 
Then x4 replaces the basic x7 to become the pivot element. 
ay v2 X3 v4 a5 Ba 
Ly 1 0 0 0 3 0 3 
x2 0 1 0 0 0 4 0 
3 1 0 0 0 —3 —3 11 
LA 0 0 0 1 -1 -1 
-2 -1 0 0 0 0 0 
Business mathematics, Integer programming, 6°” December 2005 -—28- From 5‘” 
November 2005 , as of 26°” April, 2006 


240 26" April, 2006 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Next, x, remains basic and becomes a pivot element. 


v1 T2 v3 v4 v5 v6 
ry i 0 0 0 $ 0 3 
re 0 1 0 0 0 3 0 
r3 0 0 0 0 -1 —3 8 
r4 0 0 0 1 -1 -1 1 
0 -1 0 0 ; 0 6 
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This becomes 


v1 r2 v3 v4 v5 v6 
ry i 0 0 0 ; 0 3 
2 0 i 0 0 0 $ 0 
v3 0 0 0 0 -1 —3 8 
v4 0 0 0 1 -1 -1 1 
0 0 0 0 5 5 6 
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The optimum point for Programme 3 is then, x] = 3, 73 = 8, 74 = 1, 
v5 = @E = x6 =O and z* = 6. 


Therefore the solution to the original problem Programme 1 is x} = 3, x5 = 0 
at the objective value z* = 6. 


# 
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Definition 95 transportation problem 
A transportation problem involves m sources each of which supplies a;, i = 
1,...,m, units of a certain product, and n destinations each of which requires 
bj, i =1,...,n, units of the same. The problem may be stated as following. 
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m n 
maximise: z= S- CijLij 


subject to: Say = 0, “4S ee em 


with: all x;; non-negative and integral 


The total supply and the total demand are assumed to be equal. Were this 
not so, a fictitious destination or a fictitious source is added. 
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Definition 96 north-west corner rule 


The north-west corner rule finds an initial basic solution for the transportation 
algorithm of the integer programming. It begins with the (1,1) cell in the 
m xn table, and allocates as many units as possible to x1; violating neither 
the constraints of supply, that is the summation along each row, nor those of 
demand, that is the summation along each column. Then carry on moving 
for each step either right or downwards, until we reach the lower-right corner, 


XLmn-: 
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Definition 97 loop 


A loop, which is a sequence of cells in the table used for finding the solution 
in the transportation problem, has the following properties. 
a. each pair of consecutive cells is on either the same row or the same 
column 
b. no three, or in fact any odd-numbered, consecutive cells lie in the 
same row or column 
c. the first and the last cells are on the same row or column 
d. the path along the loop is self-avoiding, that is no cells appear more 
than once in the sequence 
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Algorithm 6 Transportation algorithm. 


while optimal solution not attained do 
find an initial, basic feasible solution using, for instance, 
the North-west corner rule 
let either u; = 0 or vj; =0 depending on whether 
the i**-row or the j**-column 
has the maximum number of basic solutions 
find all u; and vj, =1,...,mandj=1,...,n 
from u; +; = cj for basic variables, and 
from cj; — u; — v; for non-basic variables 
improve the solution 
endwhile 
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Note 5 optimum for transportation problem 


In a transportation problem, optimal solution is achieved when 


Cij — Uji — Uj > O 


for all transportation costs per unit c,; of all non-basic variables. 
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Definition 98 Present and future values 


The present value po or the principal is the amount initially borrowed or 
invested. The future value p; is the principal after a period of time t. 
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Definition 99 Annual percentage rate 


Interest rates expressed per annum are called nominal rates, i. The annual 
percentage rate or effective annual rate ig is the equivalent annual rate of 
different interest rates variously compounded. 
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Definition 100 Sequence 


A sequence is a list of numbers which follows a definite pattern. It is called an 
arithmetic sequence if each of its terms is obtained from the term immediately 
preceding it by an addition of a constant d, which is called the common 
difference. It is called a geometric sequence if each of its term is obtained 
from the previous term by a multiplication of a constant r, the common 
ratio. 
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Definition 101 Series 


A series is the sum of the terms of sequence. It is called a finite series is 
one whose number of terms is finite, otherwise it is called an infinite series. 
An arithmetic series or arithmetic progression is the sum of the terms of an 
arithmetic sequence. Likewise a geometric series or geometric progression is 
the sum of the terms of geometric sequence. 
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Theorem 48 Arithmetic series 


The value of the n* term of an arithmetic series is 
Tn =at(n—1)d 
The sum of its first n terms is 


Sn = 5 2a+ (n—1)4) 
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Theorem 49 Geometric series 
The n*" term of a geometric series is 
T,, =ar"+ 
The sum of the first n terms of it is 
Sa =atar+---+ar™ 


_a(l—r”) 

— Jer 

_ a(r”™—1) 
7A 


When the number of terms approaches infinity, the summation in cases where 


r <1 becomes é 


l-r 
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Definition 102 Simple and compound interests 
A simple interest is a fixed percentage of the principal paid to an investor 
each year. A compound interest is an interest paid on the principal plus any 
interest accumulated in previous years. 
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Theorem 50 Present value for simple interest 


The present value for simple interest is 
Pt = po(1 + tt) 


where 7 is the interest rate and t the time in years. 
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Theorem 51 Present value, compound interest 


The present value in the case of compound interest is 


pe = po(1 + i)* 
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Note 6 Compounding more than once a year 


The interest may be compounded more than once a year, for example bian- 
nually, quarterly, monthly, weekly, daily, or continuously. Each time period 
is called a conversion period or interest period. The interest rate applied at 
each conversion is /m, where m is the number of conversion periods per year. 
The number of conversion periods over t years is then n = mt. 
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Theorem 52 Present value, compounding several times per year 


The present value at the end of n conversion periods is 


‘ n : mt 
i i 

p= m (1+) =m (1+5) 
m m 


where all the variables and parameters are as previously defined. 
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Theorem 53 Present value, continuous compounding 


When the number of compoundings per year becomes very large, the present 
value becomes 


Pt= poe” 


Proof. Since p; = po (1+ aye and limtm-—soo (1+ a\™ = e’, we have the 
proof. q 
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Theorem 54 Annual percentage rate 


The annual percentage rate when compounding occurs m times per year is 


ig= (145) —1 
m 
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Date 
25 Oct 2005 
1 Nov 2005 
8 Nov 2005 
15 Nov 2005 
22 Nov 2005 
29 Nov 2005 
6 Dec 2005 
20 Dec 2005 
10 Jan 2006 
17 Jan 2006 
24 Jan 2006 
31 Jan 2006 
7 Feb 2006 
14 Feb 2006 
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Appendix 


Course Outline 


Topic of lecture 
Graph and derivative 
Calculus of multivariable functions 
Exponential, log and nonlinear functions 
Matrix 
Linear algebra 
Linear programming 
Integer programming 
Financial mathematics 
Integral calculus 
Integral calculus 
Simultaneous equations 
Differential equation 
Trigonometric functions and power series 
Differential equation 


Materials 


26" April, 2006 


= 
S 
Ps 


WWWWWWwWwWwWwWwwwww iw 
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Content Date begun Dated Last updated 
Graph and derivative 20 Oct 05 25 Oct 05 10 Dec 05 
Calculus of multivariable 

functions 20 Oct 05 1 Nov 05 18 Nov 05 


Exponential-, logarithmic 
and nonlinear functions 28 Oct 05 8 Nov 05 18 Nov 05 


Matrix 5 Nov 05 15 Nov 05 10 Dec 05 
Linear algebra 5 Nov 05 22 Nov 05 10 Dec 05 
Examples for 

linear algebra 17 Jan 06 31 Jan 06 31 Jan 06 
Exercises for 

linear algebra 7 Feb 06 7 Feb 06 7 Feb 06 
Linear programming 5 Nov 05 29 Nov 05 10 Dec 05 
Examples for 

linear programming 5 Nov 05 6 Dec 06 10 Dec 06 
Integer programming 5 Nov 05 6 Dec 05 11 Dec 05 
Financial mathematics 5 Nov 05 13 Dec 05 20 Dec 05 


Examples for 
financial mathematics 18 Feb 06 20 Feb 06 20 Feb 06 
Simultaneous equations 5 Nov 05 24 Jan 06 26 Jan 06 


Differential equation 5 Nov 05 31 Jan 06 31 Jan 06 
Integral calculus 5 Nov 05 10 Jan 06 13 Feb 06 
Integral calculus 6 Feb 06 7 Feb 06 7 Feb 06 
Examples for 

integral calculus 17 Jan 06 7 Feb 06 7 Feb 06 
Exercises for 

integral calculus 10 Jan 06 7 Feb 06 7 Feb 06 
Difference equation 5 Nov 05 20 Feb 06 19 Feb 06 
Examples for 

difference equation 19 Feb 06 20 Feb 06 20 Feb 06 

Appendix 

Midterm examination 19 Dec 05 20 Dec 06 20 Dec 06 

solution 5 Jan 06 
Quiz 1 30 Jan 06 31 Jan 06 13 Feb 06 

solution 16 Feb 06 
Quiz 2 7 Feb 06 7 Feb 06 7 Feb 06 

solution 23 Feb 06 
Quiz 3 13 Feb 06 14 Feb 06 13 Feb 06 

solution 14 Feb 06 
Final examination 19 Feb 06 20 Feb 06 20 Feb 06 

solution 22 Feb 06 


Teaching method 


There were lectures and practices in class. All practice exercises, quizzes 
and exams were done in an opened-book manner. Practice of problem and 
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exercise in class were preferred to assignment and homework, since in the 
latter many students did not do the work themselves. In some of these practice 
the questions for the students were all different in order that they may start 
to think for themselves. 


Teaching media 


A camera projector and a microphone were the hardware media used. The 
material used as a media were lecture projections in this collection and lecture 
hand-outs, which will go into another collection. Because of the distance be- 
tween the class and its lecturer, feedbacks of homework were sent to students 
by email. Textbooks on Calculus from the library were used in some of the 
exercise sessions, for the students to learn how to use books as a resource to 
help them solve problems and answer questions. 


Evaluation methods 


means per cent 
Attendance 10 
Homework 10 
Quiz 1 10 
Quiz 2 10 
Quiz 3 10 
Midterm exam 20 
Final exam 30 


Evaluation will be based on the distribution plot, as relative performance 
among students. 


Bibliography 
Frank Ayres, Jr. Theory and problems of Differential Equations. Schaum’s Outline 
Series, 1981(1952) 


Teresa Bradley and Paul Patton. Essential mathematics for economics and business. 
"4 ed. 2002 

Richard Bronson. Theory and problems of operations research. Schaum’s outline 
series, McGraw-Hill, Singapore, 1982 (1983) 

Edward T Dowling. Introduction to mathematical economics. Schaum’s outline 
series,2°4 ed. 1992(1980) 


Edward T Dowling. Mathematical methods for business and economics. Schaum’s 
outline series, 1993 


David Kincaid and Ward Cheney. Numerical analysis. Brook/Cole, 1991 
Erwin Kreyszig. Advanced engineering mathematics. i ed, 1993 


G F Simmons. Introduction to topology and modern analysis. McGraw-Hill, Singa- 
pore, 1963 


George B Thomas, Jr and Ross L Finney. Calculus and analytic geometry. he ed, 
1992 


God’s Ayudhya’s Defence 26"" April, 2006 255 


Business Mathematics, notes and projections kit Tyabandha, PhD 


Quiz 1 
Business Mathematics 
31°* January 2006 
Time: 1 hour (10-11pm) 


1. Let 
3 SL 
A=|[|2 4 2 
1 1 3 
Find all the eigenvalues and a basis for each eigenspace. If possible, find 
invertible matrices P such that P~1AP is diagonal. 
Solution. Form the characteristic matrix, 


tl-A=|-2 t-4 -2 
—1 -1 t-2 


t — 3) ((t — 4)(¢ — 3) — 2) + ((t — 3)(—2) — 2) - (2+t-4) 


( 
t 
= (t — 2) (t? — 8t + 12) 
(¢- 2-6) 


Therefore the eigenvalues are 2 and 6. 


# 
For the eigenvalue 2; 
2-3 -1 —1 x -1 -1 -1 x 0 
-2 2-4 -2 y}=|-2 -2 -2 y}= {0 
1 1 2-3 z -1 -1 -1 z 0 


which gives x+y+z=0. The system has thus two free variables, for example 
{x =1,y=0,z=-—1} and {x =1,y=-—1,z=0}. Let u = (1,0,—1) and 
v = (1,-1,0). Then all linear combinations of u and v forms a set of all 
possible eigenvectors, which together with zero vector comprise the eigenspace 
corresponding to the eigenvalue 2. In other words, u and v form a basis of 
the eigenspace of the eigenvalue 2. 


# 


256 26" April, 2006 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


For the eigenvalue 6; 


te ue x 8. Sh -t\ fe 0 
a: DGiacd 29 y)=[-2 2 -2]}[/y]= [0 
1 toa ke af on 2 z 0 


Use Gaussian elimination to find the solution of this system of equation. First 
form an augmented matrix 


: es aes 8 
me, 5 ae 
Bhs ede. 25 
(esta G Sa eet ad 3.0): 
a4 ik oe 5G) 
5 Or ae 20 
i ig 
aa); 
Ie. Aihe SSS SG 
ao Oe EG 
hy ah - Be 


Ci Oak OF Se AO Te 282s ae 


1 1 -3 0 
0 4 -8 0 
0 -4 -8 0 
z(11), (1 —2 0); (III) + 4(ID), (-4 8 0)+4(1 2: <0); 
1 1 -3 0 
0 1 -2 0 
0-0 20°. 0 


Therefore rank of A is 2. The system has only one free variable, that is to 
say, one independent solution. Any particular non-zero solution generates its 
solution space, that is the eigenspace, for example x = 1, y = 2 and z = 1. 
So w = (1, 2,1) forms a basis of the eigenspace of the eigenvalue 6. 


# 
1 1 1 

Let Pbe | 0 -1 2 | Since det P = 1(—-1) —1(24+ 1) = —4 £0, hence 
-l1 OoO 1 


P has an inverse. Therefore P is the required matrix such that P~!AP is 
diagonal. 


# 
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Quiz 2 
Business Mathematics 
7” February 2006 
Time: 1 hour (10-11pm) 


1. Evaluate the following integrals. [10] 


1 2 14 2 1 
- ada, | (a _ 3) dz, / —dz, ‘ (ar + Vz) dx, i xe” dx 
0 0 -14 0 0 


Solution. 


# 
In 4 tdz the limits of integration pass through 0, where there is disconti- 
nuity and where 1/z is undefined. 


Further, lim,_,9- f(z) = —oo and lim,_,9+ f(z) = co. We look at 


1 

1 
| “dz = lim (Inz|j = lim (In1 — Ina) = 0— (—co) = 
O22 


a0t aot 
and : , , 
fh -de= f vde+ f —dz 
-1% -14 0 
where 
fre Hi, | eg) a ines heii) 
-azZ = lm = = hm n(—Z = i nho— in 
Ag 8 pe} ele) BO ~1 650- 
= lim (Inb 
oe 


the latter of which is undefined. Therefore ier tdz is also undefined. 


# 
. ee et ee ee 
2 d = Eat = 3/2 = } = 
| (x + \/z) z 3 Tige ay 3 3 
# 
Let u = x and dv = e*dx. Then du = dz and v = e*. And then, 
1 1 
[ serac= (eer - [ e"dx = (xe —e*|, = 1 
0 0 
# 
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Quiz 3 


Business M 


athematics 


14‘” February 2006 
Time: 1 hour (10-11pm) 


Choose only one problem[10], either 
1. Solve 


y+3z= 
24+ 2y—z=8 
—“2£+5z=8 


by Gaussian elimination. 


Solution. Write an augmented matrix, 


Or. 
2: 2 
=f) 0 
(1) (11D, (0 1 3 9)a(-1 0 
21620 
2 2 
0 1 


“1D, 1(2t 02.8.8). —20): 


iPS. TE, (0 2-9: 24) 
2(1 3 9); 


Therefore, directly we have z = 2, y= 


or 
2. Solve 


3 9 
-1 8 
5 8 
5 8); 
5 8 
= 
3 9 


Qo =f B= 0 


—5 —-8 
9 24 
3. (CO@d 
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3 9); (III) — (1D), (2 9 24)- 


—5 —8 
3 39 
3. (6 
9 


2% +y-—2z=10 


3a + 2y 4 
ox + 4y 4 
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+ 3z = 4 
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3(2) =3 and z = —8 + 5(2) =2. 


# 
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by any method. 


Solution. Form an augmented matrix, 


kit Tyabandha, PhD 


=1 5); (i= 


ae eee) 
a. 2> 9 
oe 

31), 5 —-1 5);UD-3(), 3 2 

5(1),(5 4 3 4)-5(1 $ -1 5); 
i ; =I 
0 3 5 
03 8 

2(1I), (1 10 —28); (111) — 3(Z1), (3 
1 4 -1 
0 1 10 
0.0 ee 

Directly, z y 3, Y 28 
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—28 
21 


10(—3) = 2 and x = 5— $(2)+(-3) =1. 


# 
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Midterm examination 
Business mathematics 


20°" December 2005 
Time: 3 hours (14pm) 


1. 
. Given 
f(z) = ax? + ba+e 


where a = c= land b= 2. Find all the z-intercepts, y-intercepts, and critical 


points. 
. Given 

Iny=1-<2 
Find y. 


. Given 


Find 0. 


. Given 
logy = 3loga+logb — loge 


Find y. 


Solution. 


f(v)=2? +2241 


x-intercept, y = 0; 
x? +22+1=0 


(c+1)(x+1)=0 


xr=-l 
f(x) touches x-axis at x = —1 
Critical point; f’(x) = 0; 
2x+2=0 
g£=-1 
y-intercept, x = 0; 
y=l 
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y= el 
# 
qgtty-2)n = a? 
b=(a+y—2z)n 
# 
log a® + log b — loge = logy 
3b 
log a log y 
c 
a°b 
yr 
c 
# 
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2. 


Find the first- and second-order partial derivatives. 


Solution. 
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q=Inax+lIny 
gaa tar?4+e4 Qeyt ay’ 


p= 150e9:74t 


| 
Sle 


2 
Og 0 
OyOu 


— = 8e7 +26 +1+ 2y+y? 
Ox 


Oz 


Tae | 
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oe 
Oz 
gE er 
# 
072 
=2+2 
OxOy Be! 
# 
O2z 
ay 
oS 
Oz 
=2+2 
OyOu wad 
# 
p = 150(0.74)e°"* = 11169 
sia 
Op 0.74t 0.74t 
na = 111(0.74)e = 82.14e 
Se 
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3. Find the determinant of the following matrices. 


a. 
a b 
c d 
[2] 
b. 
a be 
de f 
g h 2 
[3] 
Cc. 
a bed 
e f g h 
ij kil 
mn o p 
[5] 
Solution. 
a. 
a b 
: | = ad — be 
# 
b. 
a be 
de fl=al, i|-2|2 Ae hae 
ee g g 
=a(ci — fh) — b(di — fg) + c(dh — eg) 
=aci —afh — bdi + bfg + cdh — ceg 
# 
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senoe 
ore a 


a 
e 
a 

m 


é JF 9 
ij k 
m n O 


foh 
ae alle 
n p 


= a(f(kp — lo) — g(jp — In) + h(jo — kn)) 
b (e(kp — lo) — g(ip — lm) + h(to — km)) 
t e(e(jp — In) — f (ip — Im) + h(in — jm)) 
d(e(jo — kn) — f(io — km) + g(in — jm)) 
= afkp — aflo— agjp + agln + ahjo — ahkn 
bekp + belo + bgip — bglm — bhio + bhkm 
+ cejp — celn —cfip+cflm+ chin — chjm 
dejo + dekn + dfio — dfkm — dgin + dgjm 
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4. The relationship between the total revenue r;, the price p, and the output 
quantity q is 


Tt = pq 


The demand function is p = a — bq, where a and 6 are positive constants. 


Find r,;, the marginal revenue r,,, and the average revenue rz. Then find q 
at the maximum r;. And then sketch the graphs of r;, ry, and rg. (assume 
ag > 4b) 


[10] 
Solution. 
re = pq = (a — bq)q = aq — bq? 
oS 
dr; 
m= =a— 2b 
r a5 a a 
oe 
Ta —- rt — Pp =a—- bq 
qd 
# 
At maximum r;; 
r, =a—2bq=0 
_ a 
I 3 
# 
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Tt, Pm and rq (a2 > 4b,a=5,b=1) 
8 T T T T T 


2 
maxr;, = G 


1 0 1 2 3 4 5 6 
qd 
Figure 14 
a (<) b a? as 1 b i? 
a a O) 402) ~ 462 2) 4b 
q=0; m=0 
r, = 0; 
q(a — bq) = 
a 
q=0, b 
q=0; Tm=a 
a 
m = 0; ey 
I 35 
q = 0; Ta =a 
a 
a = 0; > a 
r q 5 


This can be summarised as Figure 14. 
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5. 


. Let A oS B ee) and C GS) represent the second-order condi- 


9 


tions of critical point of function. Suppose the graph of a function has the 
shape as shown in Figure 15. 


Figure 15 
Which of these conditions is satisfied at points J, P and T’?[3] Explain. [2] 


Find f’(z) and f”(«).[1] Show that 1, —1 and 2 are the critical points. [2] 
Which of these are maximum, minimum or inflection point? [2] 


Solution. 
I ¢B 
P ¢C 
T ¢A 


At I, P and T, f’(-) =0. 


At I, PandT, f’(-)=0, <0, and > 0 respectively. 


fi(z) =32+-— 925 a 


f(z) = = 1243 — 2727 + 6r +9 
# 
f'(x) = (3x? — 6x + 3) (w +1)(a — 2) = 3(@ — 1)?(x + 1)(a — 2) 
Critical points, 1, -1 and 2. At these points f’(-) = 0. 
# 
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f"-) =0 +1. inflection point 
f’(-1) <0 —-I1 maximum point 
f"(2)>0 +2 minimum point 
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6. Solve the following programme by the simplex method. 


maximise: z= 327, + 4%2 + 5x3 


subject to: 2, +2%2+2%3 <2 


@+%24+ 3243 <1 
321 + 2%. +23 <4 


with: all the variables non-negative 


[10] 
Solution. Draw simplex tables. 
Ly v2 x3 v4 XS X6 
3 4 5 0 0 0 
x4 0 1 1 1 1 0 0 2 
Xs 0 1 1 3 0 1 0 1 
x6 0 3 2 1 0 0 1 4 
-3 -4 -5 0 0 0 0 
Hine v2 x3 v4 “5 rE 
2 2 1 5 
4 3 3 0 1 —3 0 3 
1 1 1 1 
a 2 = se Se i 
v6 3 3 3 3 
4 7 5 5 
a | 0 0 3 0 3 
Ly r2 v3 v4 U5 v6 
L4 0 0 -2 1 - 0 1 
x2 1 1 3 0 1 0 1 
x6 1 0 -5 0 5 1 2 
1 0 7 0 4 0 4 
eH 1,73 =1, 7% =2, c7 = 273 = 25 = 0 and 2* = 4. 
# 
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Final Examination 
Business Mathematics 


20° February 2006 
Time: 3 hours (9-12am) 


1. 
. Give the definition of a homogeneous function. [2] 
. What are order and degree of a differential equation? [2] Given the differential 
equation 
yl +4(y") + (y')? = sine 
Give its degree and order. [2] Is this equation homogeneous? [1] 


. A primitive gives rise to a differential equation. Explain how the primitive 
3 
y = Ax? + Bx +C gives rise to the differential equation +$ = 0. [3] 


Solution. 
. A function f(z,y) is said to be homogeneous of degree n if f(Ax, Ay) 
A” f(xy). 


I 


. The order of a differential equation is the order of its highest derivative. 


# 


The degree of a differential equation is the degree of its highest ordered deriva- 
tive. 


# 
The degree is 1, and the order 4. 
# 
The equation is not homogeneous. 
# 
d 
“Y ~ oA +B 
dx 
d*y 
—.~=2A 
dx? 
d°y 
ae ae) 
dx? 
# 


2. 

. Explain what sequence and series are. When are they said to be arithmetic?, 
when geometric? [6] 

. Given an arithmetic series, 1+3+5+7+--:. What is the value of the n'® 
term? Find the value of the 50°" term. What is the sum of the first n term? 
Find the sum of the first 50 terms. [4] 
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. Given a geometric series, 1 + 3 + + + 3 +--+, What is the value of the n** 
term? Find the value of the 50°" term. What is the sum of the first n terms? 
Find the sum of the first 50 terms. What is the sum to infinity, that is to say, 
Soo? [5] 


Solution. 
. A sequence is a list of number following a definite pattern. A series is the sum 
of the terms of sequence. They are said to be arithmetic if each of their terms 
is obtained from the term immediately preceding it by an addition of some 
constant. They are said to be geometric if each of their terms is obtained 
from the previous term by a multiplication of some constant. 

# 
. From what is given, a = 1 and d = 2. The value of the n*" term is then T,, = 
a+(n—1)d = 14+(n—1)2d, while that of the 50° term Tso = 1+(50—1)2 = 99 
The sum of the first n terms is S, = } (2a+ (n—1)d) = $(2+(n-—1)2), 
and the sum of the first 50 terms is S59 = 32 (2(1) + (50 — 1)2) = 2500. 


# 
. From the given geometric series, a = 1 and r = $. Then J, =ar™ 1) = 
vag _4y50-1 ae WEG) i = 
(3)" 1 Tyo = (3) = ga, Sn = SS = A = 2 (1 - oe), St0 = 
50 
2(1— 950) = 2yas* and Sxo = 72; = Gy = 2. 
# 


3. 

. Give the formula for finding the present value of compound interest. Let the 
annual interest rate be 10 per cent. What is the present value of 1,000 Bahts 
in ten years’ time? [4] 

. What is the present value at the end of n conversion periods in t years? Here 
m. is the number of conversion periods per year. Find the present value of 
1,000 Bahts in five years’ time, when m = 5 and i = 10 per cent. [4] 

. Compare the present values obtained from (a) and (b), and discuss the dif- 
ference between them. [2] 


Solution. 
. The formula for the present value of compound interest is p, = po(1+i)*. From 
what is given, i = 0.1, pp = 1000 and t = 5. Therefore, ps = 1000(1+0.1)° = 
1610.51. 

# 


. The present value of compound interest at the end of n = mt conversion 
period is py = po @ + oe From what is given, 7 = 0.1, m = 5, po = 1000 
and t = 5. Therefore ps; = 1000 (fi + ot) = 1485.95 Bahts. 


- pPe> Pt,m 
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A, 

. What is the average or mean value of an integrable function f(x) on [a,b]? 
| 

. What is an integrable function? [2] Give an example of a function that is not 
integrable on some range. [2] 

. Consider a (2 x 2) system of linear equations in the slope-intercept form, 


y=myaxt+ by 


y = mga + bg 


Give conditions for the system to have a unique solution, no solutions, and 
infinitely many solutions. [4] 


Solution. 
. The average or mean value is 


b 
[fae 
# 


. An integrable function is a function whose value never becomes infinitely large 
on the range of interest. 


# 


The function f(x) = 4 is not integrable over the range [—1, 1]. 


. The system has a unique solution when m, 4 mga, no solutions when b; ¥ be, 
and an infinite number of solutions when m; = m2 and b; = bo. 


# 


5. 
. Explain revenue and elasticity mathematically and verbally. [1] Explain 
marginal and average costs both mathematically and also in words. [1] 


b. Give the condition for a° = 1 to be true. [1] 


. The following is a formal definition of limit. 
For a function f(x), lim,_,, f(z) =1 if and only if for every « > 0, there 
exists 6 > 0 such that |f(a) — 1] < « whenever 0 < |x — al < 6. 


Explain in words what we mean by this. Give your answer both in English 
and in Thai. [1] 

. What is a cut algorithm? [1] In the following solved example, identify all pivot 
elements and explain in detail how the problem was solved. [5] 


maximise: z= 2%,+ 72 
subject to: 22, + 5% < 17 
321, + 2x2 < 10 


with: 21,22 non-negative and integral 
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Use cut algorithm. 
Solution 


Find the first approximation of Programme 1 normally using the simplex 
method. 


Ty v2 v3 v4 
2 0 0 
23 0 2 5 1 0 17 
wa 0 3 2 0 1 10 
-2 -1 0) 0 0 


Since 2 << q, we know that 3 is the pivot element, and therefore we replace 
the basic variable x4 with x1. 


71 r2 v3 v4 
a er 
bs 1 2 0 3 8 
a z 
We have a7 = 2, L3 = a v5 = 2, = 0 and z* = 2. Since both xj and «3 


are non-integers, arbitrarily choose the former to generate a new constraint. 
Then our Programme 2 becomes, 


52 wi = AO ai 
Ly ge gts 3 
2 1 1 
gre + 304 373 Ly 
2 1 a, 

372 34 3- 
ee 1 

x —2x = 
aie eae 


2%2+%4>1 
and our new programme becomes 


maximise: z= 2%, +%2+ 0x3 + 024 
11 2 31 


bject to: a 
subject to 3 72 gt 3 
* 10 
x =X. - -@4 = 
sea iar ac 
2% +44>1 


with: all variables non-negative and integral 
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ry r2 r3 v4 U5 v% 
2 1 0 0 0 —M 
v1 0 1 2 0 $ 0 * 
ts 0 i a mh: a 0 3 
XG —-M 0 2 0 1 -1 1 1 
-2 -1 0 0 0 0 0 
0 -2 0 -1 1 -1 -1 
Now 22 replaces zg in the basic variables and becomes the pivot element 
Ly ip) X3 XA a5 XE 
X1 1 0 0 0 3 3 3 
£3 0 0 1 i tt a Md 
m | o 1 0 f=} F 
1 1 1 1 
-2 0 0 3 3 2 2 
0 0 0 0 0 0 0 
This becomes, 
ry v2 x3 X4 x5 
Ly 1 0 0 0 3 3 
v3 0 0 1 -3 i 4 
wo Ge Sh. Bt 1 
jf ge — a= B 
Then our first approximation of Programme 2 is xj = 3, x3 T L3 w 
ri = vf = 0, and 2* = #8. Arbitrarily choose x3 to generate the new 
constraint. 


1 1 
BET 5 C8 Dts 5 
i 1 
v4 v5 2 v2 
d A sy 
2 pe ea 
94 975 a= 


4-45 >1 
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Then our Programme 3 becomes, 


maximise: 


subject to: 


with: 


z= 20, 4+ %2 + Org + 0x4 + Or5 


1 
co rd 
Bie: Ag. ke 
TIO Dee gee SD 
fe ee 
XQ gv gue 5 


%4—- 2X5 >1 


all variables non-negative and integral 


We draw our table for this programme. 


Ly v2 L3 vA a5 X6 v7 
2 1 0 0 0 0 —M 
x4 0 1 0 0 0 5 0 0 3 
r2 0 0 1 0 4 4 0 0 4 
r3 0 1 0 0 -3 2 0 0 z 
U7 —M 0 0 0 1 -1 -1 1 1 
-2 -1 0 0 0 0 0 0 
0 0 0 -1 1 -1 -1 
Then x4 replaces the basic +7 to become the pivot element. 
Ly Zr v3 v4 a) v6 
X1 1 0 0 0 ; 0 3 
re 0 1 0 0 0 3 0 
x3 1 0 0 0 —3 - 11 
v4 0 0 0 1 -1 -1 1 
-2 -1 0 0 0 0 0 
Next, x, remains basic and becomes a pivot element. 
Ly x2 x3 x4 x5 6 
X1 1 0 0 0 3 0 3 
Le 0 1 0 0 0 3 0 
£3 0 0 0 0 -l - 8 
x4 0 0 0 1 -1 - 1 
0 -1 0 0 3 6 
This becomes 
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Ty v2 x3 L4 Xs X6 
Ly 1 0 0 0 3 0 3 
re 0 1 0 0 0 3 0 
x3 0 0 0 0 -1 —3 8 
L4 0 0 0 1 -1 -1 1 
0 0 0 0 3 $ 6 
The optimum point for Programme 3 is then, zj = 3, 73 = 8, x7 = 1, 


v5 = x8 = x4 = 0 and z* = 6. Therefore the solution to the original problem 
Programme 1 is x} = 3, x3 = 0 at the objective value z* = 6. 
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Students’ scores 
Business Mathematics 
2005-6 


Kit Tyabandha, PhD 
26%” April, 2006 


Introduction 


This is a report of study and examination results for the course Business 
Mathematics taught by me during the second semester, 2005-6 academic 
year. It has been my concern from the beginning of 2006 that students had 
not done their homework themselves, therefore I have had since then practices 
and quizzes in class, with the hope that they would obtain some problem- 
solving skill. Everything, which includes exams, quizzes and practices, was 
done in an opened-book manner because that is how the real world works and 
that is where the progress lies. 


Midterm exam 


There were six questions in the midterm exam, namely from M1 to M6. MT 
is the total marks and MM the actual contribution toward the students’ final 
score. Table 6 shows the six marks, M@1—M6, for each student ID, together 
with the total mark and the same in per cent for each. 


ID M1 M2 M8 M4 M5 M6 Total (60) Per cent (20) 


4661 0.5 83 52 0O 2.2 1.5 17.7 5.9 
4801 6 7.5 5.2 0 2.7 3 24.4 8.1 
4802 7 64 5.2 0 2 3 23.6 7.9 
4803 6.5 76 52 35 3 1.50 27.3 9.1 
4804 

4805 0 5 5.2 0 0 4.7 14.9 5 
4806 7 3.3 5.2 0 1 2 18.5 6.2 
4807 7 0.8 5.2 1 1 0 15 5 
4808 1 42 5.2 1 1 0 12.4 4.1 
4809 1 5 5.2 1 1 2 15.2 5.1 
4810 1.5 5 5.2 3.5 1 0 16.2 5.4 


Table 6 Midterm marks, Business Mathematics, 2005-6, contribution of 
20 per cent toward the overall points 
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ID 

4811 
4812 
4813 
4814 
4815 
4816 
4817 
4818 
4819 
4820 
4821 
4822 
4823 
4824 
4825 
4826 
4827 
4828 
4829 
4830 
4831 
4832 
4833 
4834 
4835 
4836 
4837 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 
4848 
4849 
4850 
4851 
4852 


Table 6 (continued) Midterm marks. 
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M1 
6 


M2 
7.6 


7.6 
9 
8.7 
6.8 
1.7 
10 


M8 
5.2 


5.2 
5.2 


OAM MAMAS 
oR WY) Dow wh bw 


i) 


OOMANMOMOOOSO 
NS) 


hehe 
ww on 


ww w 
to 


“I 
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Total (60) Per cent (20) 
19.8 6.6 
21.3 7.1 
32.4 10.8 
25.3 8.4 
10.9 3.6 
24.7 8.2 
16.8 5.3 
10.1 3.4 
12.9 4.3 
12.8 4.3 
30.1 10 
29.6 9.9 
20.1 6.7 
32.3 10.8 
3.2 10.8 
30.2 10.1 
11.3 3.8 
13.9 4.6 
12 4 
13 4.3 
12.8 4.3 
14.5 4.8 
30.4 10.1 
31.1 10.4 
31.6 10.5 
33.1 11 
38.1 12.7 
26.6 8.9 
37.6 12.5 
11.5 3.8 
17.3 5.8 
18.7 6.2 
32.3 10.8 
36 12 
31.4 10.5 
25.6 8.5 
25.3 8.4 
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ID M1 M2 M3 M{ M5 M6 Total (60) Per cent (20) 
4853 


4854 6 3.9 5.2 3 3 1.5 22.6 7.9 
4855 8.5 7.5 52 3 3 1.2 28.4 9.5 
4856 8 7.5 7 05 05 15 25 8.3 
4857 7 8.3 8 05 2.7 1.5 28 9.3 
4858 85 7.9 8 1.5 2.7 1.5 30.1 10 
4859 8 83 7.8 0.5 2.7 1.5 28.8 9.6 
4860 

4861 5 8.3 8 0 2 3 26.3 8.8 
4862 85 29 52 3 2.5 1.5 23.6 7.9 
4863 8 2 5.2 3.5 2.5 1.5 22.7 7.6 
4864 

4865 65 66 52 35 3 1.5 26.3 8.8 
4866 6 5.1 5.2 3.5 3 1.5 24.3 8.1 
4867 7 2.1 5.2 3 3 1.5 21.8 7.3 
4868 7 2.2 65.2 3 2.5 1.5 21.4 7.1 
4869 

4870 8 79 5.2 3 0.7 1.5 26.3 8.8 
4871 6 76 8 3 5 1.5 31.1 10.4 
4872 7 7 8 3 2.7 1.5 29.2 9.7 


Table 6 (continued) Midterm marks. 


Table 7 shows mark and rank of each student. 
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ID 

4661 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
4811 
4813 
4814 
4815 
4816 
4817 


5.07 
5.4 
6.6 
7.1 
9.4 
10.8 
8.43 
3.63 


ID 

4818 
4819 
4820 
4821 
4822 
4823 
4824 
4825 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 


53 
46 
51 
AT 
49 


ID 

4836 
4837 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 
4848 
4850 
4851 
4852 


Table 7 Mark and rank of students’ midterm scores. 


SCOTE 
4.83 
10.13 
10.37 
10.53 
11.03 
12.7 
8.87 
12.53 
3.83 
5.77 
6.23 
10.77 
12 
10.47 
8.53 
8.43 
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ID 

4854 
4855 
4856 
4857 
4858 
4859 
4861 
4862 
4863 
4865 
4866 
4867 
4868 
4870 
4871 
4872 
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10 12 14 16 
score 


Figure 16 Distribution of students’ midterm scores. 


20 


The mean score for the midterm exam is 7.78, median 8.18, minimum 3.37 
and maximum 12.7. The standard deviation is 2.53. 
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ID 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
4811 
4813 
4814 
4815 
4817 
4818 
4819 
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Quiz 1 score distribution, Business Maths, 2005-6 
T T 


number of students 


5 
score 


Figure 17 Distribution of students’ scores from Quiz 1. 
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ID 
4661 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
A811 
4813 
4814 
4815 
4816 
4817 


Table 8 Quiz IT marks, Business Mathematics, 2005-6, contribution of 10 per cent toward the 
overall points 


ID 
4818 
4819 
4820 
4821 
4822 
4823 
4824 
4825 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 


ID 
4836 
4837 
4838 
4839 
4840 
4841 
4842 
4843 
4845 
4846 
4847 
4848 
4849 
4850 
4851 
4852 


19 
16 


ID 
4854 
4855 
4856 
4857 
4858 
4859 
4861 
4862 
4863 
4865 
4866 
4867 
4868 
4870 
4871 
4872 
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Quiz 2 score distribution, Business Maths, 2005-6 
i) T 


number of students 


5 
score 


Figure 18 Distribution of students’ quiz 2’s scores. 
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ID Score 
4661 10 
4801 10 
4802 10 
4803 10 
4805 10 
4806 10 
4807 10 
4808 2 
4809 2 
4810 10 
4811 10 
4813 2 
4814 10 
4815 10 
4816 10 
4817 2 


Table 9 Quiz 3 marks, Business Mathematics, 2005-6, contribution of 10 per cent toward the 


overall points 


Rank 


a 


ID 
4818 
4819 
4820 
4822 
4823 
4824 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 
4836 
4837 


10 


R 
1 
1 
5 
1 
1 
1 
1 
i 
il 
1 
5 
1 
1 
3 
2 
1 


ID 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 
4848 
4849 
4850 
4852 
4854 
4855 


10 
4 
10 
5 
10 
10 
10 
5 


Rank 
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ID 
4856 
4857 
4858 
4859 
4861 
4862 
4863 
4865 
4866 
4867 
4868 
4870 
4871 
4872 
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Figure 19 Distribution of students’ quiz 3’s scores. 
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Final Exam 


The final examination was held on 20 February 2006. There were five ques- 
tions, with marks respectively from the first to the last, 10, 15, 10, 10 and 10, 
totalling 55. 


ID Q1 Q2 Q8 Q4 Q5 Total (55) Per cent Scaled (30) 
4661 2 4.9 3 0 0.5 10.4 18.91 5.67 
4801 5 0 0 4 2.8 11.8 21.45 6.44 
4802 2 6.4 2 4 3 17.4 31.64 9.49 
4803 2 0 11 4 1 8.1 14.73 4.42 
4805 4 6.5 0 6.1 4 20.6 37.45 11.24 
4806 1 12.5 0 3 2 18.5 33.64 10.09 
4807 5 6.2 0 2.1 58 19.1 34.73 10.42 
4808 3 0 0 0) 18 4.8 8.73 2.62 
4809 2 6.4 0 0.1 O 8.5 15.45 4.64 
4810 3 10.4 O 3.1 2 18.5 33.64 10.09 
4811 0 1.7 0 41 1 6.8 12.36 3.71 
4813 2 6 0 0 2 10 18.18 5.45 
4814 5 6.5 0 4 2.8 18.3 33.27 9.98 
4815 1 6.9 0 2.1 3 13 23.64 7.09 
4816 3 6.7 0 0.1 O 9.8 17.82 5.39 
4817 2 0 2 41 0.2 83 15.09 4.53 
4818 3 7 2.2 5 3 20.2 36.73 11.02 
4819 2 0 0 4 0 6 10.91 3.27 
4820 0 0.2 0 0 1 1.2 2.18 0.65 
4821 2 0 0 0 6.8 8.8 16 4.8 
4822 0 10 3.1 4 8 25.1 45.64 13.69 
4823 1 9.3 0 1 15 12.8 23.27 6.98 
4824 5 8.3 41 6.1 63 29.8 54.18 16.25 
4825 4 6 2.1 4 0 16.1 29.27 8.78 
4826 2 0.2 0 6.1 3.5 11.8 21.45 6.44 
4827 1 0 0 4 2 7 12.73 3.82 
4829 0 10.5 O 0 2 12.5 22.73 6.82 
4830 0 2 0 0 0 2 3.64 1.09 
4831 0 13.8 9 1 0 23.8 43.27 12.98 
4832 2 0 0 4 0 6 10.91 3.27 
4833 0 6.5 2 6 2 16.5 30 9 


Table 10 Final marks (continued) 
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ID Q1 Q2 Q8 Qs Q5 Total (55) Per cent Scaled (30) 


4835 1 6.5 2 4 0 13.5 24.55 7.36 
4836 3 0 4.1 51 3.3 15.5 28.18 8.45 
4837 1 2 2.2 4 1 10.2 18.55 5.56 
4838 0 0.2 0 0 1 1.2 2.18 0.65 
4839 5 125 11 21 5 25.7 46.73 14.02 
4840 3 12.9 4 3.1 2 25 45.45 13.64 
4841 4 8.8 0 A 3.8 20.6 37.45 11.24 
4842 6 13.9 3.2 4 0 27.1 49.27 14.78 
4843 3 13.5 6.2 4 6.8 33.5 60.91 18.27 
4844 7 1.7 0 0 1 9.7 17.64 5.29 
4845 2 5.5 24 0.1 0 10 18.18 5.45 
4846 0 0 0 0.1 0 0.1 0.18 0.05 
4847 0 0.5 0 0 0 0.5 0.91 0.27 
4848 1 10.1 11 3 3 18.2 33.09 9.93 
4849 2 10.1 2.2 51 48 24.2 44 13.2 
4850 3 7.3 0 A 5 19.3 35.09 10.53 


Table 10 Final eramination marks, Business Mathematics, 2005-6, con- 
tributing 30 per cent towards overall points 


ID Q1 Q2 Q8 Qs Q5 Total (55) Per cent Scaled (30) 
4851 2 0 0 4 0 6 10.91 3.27 
4852 2 0 0 4 4 10 18.18 5.45 
4854 6 0 0 6 3 15 27.27 8.18 
4855 2.2 8 1 4 4 19.2 34.91 10.47 
4856 6 0 4 4 0 14 25.45 7.64 
4857 4.1 2.2 10 6 2 24.3 44.18 13.25 
4858 0 6.8 9 6.1 2 23.9 43.45 13.04 
4859 3 6.2 0 6.1 O 15.3 27.82 8.35 
4860 4.5 135 3.2 4.1 O 25.3 46 13.8 
4861 2 6.1 0 4 3 15.1 27.45 8.24 
4862 2 9.9 0 4 4.8 20.7 37.64 11.29 
4863 3 126 28 7 3.9 29.3 53.27 15.98 
4865 4 0 0 4 2.8 10.8 19.64 5.89 
4866 3 3 2 0 0 8 14.55 4.36 
4867 0 0 0 4 0 4 7.27 2.18 
4868 0 0.2 1 0 0 1.2 2.18 0.65 
4870 0O 4,2 0 0 0 4.2 7.64 2.29 
4871 4 8.7 3.1 4 5 24.8 45.09 13.53 
4872 1 0 0 7 0 8 14.55 4.36 


Table 10 (continued) Final ecam marks. 


From this we scale the total down to 30. The mean resulted is thus 7.74, 
the median 7.23, minimum 0.05, maximum 18.27, and the standard deviation 
4.46. 
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ID 
4661 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
A811 
4813 
4814 
4815 
4816 
4817 
4818 


Score 
5.67 
6.44 
9.49 
4.42 
11.24 
10.09 
10.42 
2.62 
4.64 
10.09 
3.71 
5.45 
9.98 
7.09 
5.39 
4.53 
11.02 


44 
16 


ID 
4819 
4820 
4821 
4822 
4823 
4824 
4825 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 
4836 
4837 


ID 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 
4848 
4849 
4850 
4851 
4852 
4854 
4855 


Score 
0.65 
14.02 
13.64 
11.24 
14.78 
18.27 
5.29 
5.45 
0.05 
0.27 
9.93 
13.2 
10.53 
3.27 
5.45 
8.18 
10.47 


Score 
7.64 
13.25 
13.04 
8.35 
13.8 
8.24 
11.29 
15.98 
5.89 
4.36 
2.18 
0.65 
2.29 
13.53 
4.36 


Table 11 Final examination marks and ranks, Business Mathematics, 2005-6, totalling 30 per 


cent 
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Final exam score distribution, Business Maths, 2005-6 
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Figure 20 Distribution of students’ final examination’s scores. 


Total exam score 


All two exams and three quizzes total to 80 per cent. Of this the average score 
is at 31.63, the median 32.03, minimum 11.02, maximum 51.81 and standard 
deviation 8.85. Then the scores and their ranks are shown in Table 12 and 
the distribution of scores in Figure 21. 
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ID 
4661 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
4811 
4813 
4814 
4815 
4816 
4817 
4818 


Score 
26.77 
34.77 
35.56 
36.22 
33.8 

34.66 
34.92 
17.05 
16.6 

36.09 
31.51 
26.75 
38.58 
38.49 
29.98 
16.46 
43.55 


ID 
4819 
4820 
4821 
4822 
4823 
4824 
4825 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 
4836 
4837 


Score 
24.64 
11.02 
14.9 

37.46 
36.32 
45.92 
22.08 
39.6 

36.68 
37.28 
20.06 
31.82 
24.77 
28.73 
25.73 
30.19 
30.7 


ID 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 
4848 
4849 
4850 
4851 
4852 
4854 
4855 


Score 
24.22 
32.25 
41.47 
44,24 
35.95 
51.81 
20.92 
19.22 
24.49 
23.24 
42.23 
27.6 

41.19 
19.51 
31.09 
38.22 
33.74 


Table 12 Total exam score, Business Mathematics, 2005-6, 80 per cent of total 


Score 
30.17 
45.19 
41.77 
40.15 
13.8 

29.5 

38.36 
44.65 
31.16 
28.46 
27.35 
26.69 
32.76 
40.49 
36.3 
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Total exam score distribution, Business Maths, 2005-6 
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number of students 


Figure 21 Distribution of total scores from all three quizzes and two 
exams. 


Homework 


For the marks for homeworks the mean is 8.08, median 8.00, minimum 7.00, 
maximum 9.00 and standard deviation 0.44. 
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Attendance 


The scores for attendance have as their average 9.33, median 9.5, minimum 
5.00, maximum 10.00 and standard deviation 1.02. Scores and their ranks are 
given in Table 13. 
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ID 
4661 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
4811 
4813 
4814 
4815 
4816 
4817 
4818 


R 
8 
5 
2 
2 
1 
2 
4 
2 
2 
1 
3 
4 
2 
1 
2 
1 
1 


ID 
4819 
4820 
4821 
4822 
4823 
4824 
4825 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 
4836 
4837 


Table 13 Score and rank of attendance. 


WRNNNNFHENHF NFR Wr Why 


ID 
4856 
4857 
4858 
4859 
4860 
4861 
4862 
4863 
4865 
4866 
4867 
4868 
4870 
4871 
4872 
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ID Midterm Quiz1 Quiz2 Quiz’ Final AW & Atd 


4661 5.9 0 5.2 10 5.67 12 
4801 8.13 1.6 8.6 10 6.44 16 
4802 7.87 2.1 6.1 10 9.49 17.5 
4803 9.1 4.5 8.2 10 4.42 17.5 
4805 4.97 1.6 6 10 11.24 19 
4806 6.17 2.4 6 10 10.09 17.5 
4807 5 1.6 7.9 10 10.42 16.5 
4808 4.13 1.4 6.9 2 2.62 18.5 
4809 5.07 1.5 3.4 2 4.64 17.5 
4810 5.4 2.4 8.2 10 10.09 18 
4811 6.6 4 7.2 10 3.71 17 
4813 7.1 4 8.2 2 5.45 15.5 
4814 9.4 1 8.2 10 9.98 17.5 
4815 10.8 2.4 8.2 10 7.09 18 
4816 8.43 = 6.2 10 5.35 17.5 
4817 3.63 2.3 A 2 4.53 18 
4818 8.23 6.3 8 10 11.02 18 
4819 5.27 1.9 4.2 10 3.27 18 
4820 3.37 1 A 2 0.65 17.5 
4821 4.3 1.8 4 = 4.8 17.5 
4822 4.27 1.6 7.9 10 13.69 16.5 
4823 10.03 1.6 7.7 10 6.98 16.5 
4824 9.87 1.6 8.2 10 16.25 17.5 
4825 6.7 1.6 5 7 8.78 14.5 
4826 10.77 4.2 8.2 10 6.44 18.5 
4827 10.77 3.9 8.2 10 3.82 18 
4829 10.07 2.2 8.2 10 6.82 18 
4830 3.77 ik 4.2 10 1.09 17 


Table 14 All results compiled into a single table, Business Mathematics, 
2005-6 
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ID Midterm Quiz1 Quiz2 Quiz3 Final AW & Atd 


4831 4.63 4 8.2 2 12.98 17 
4832 A 1.5 6 10 3.27 17.5 
4833 4.33 1.4 A 10 9 17.5 
4835 4.27 1.3 7.8 5 7.36 18 
4836 4.83 2.2 5.7 9 8.45 18 
4837 10.13 1 A 10 5.56 16.5 
4838 10.37 4 7.2 2 0.65 17 
4839 10.53 2 3.7 2 14.02 18 
4840 11.03 1.6 5.2 10 13.64 LY, 
4841 12.7 4.1 6.2 10 11.24 18 
4842 8.87 1.9 5.4 5 14.78 18 
4843 12.53 2.8 8.2 10 18.27 19 
4844 3.83 2.8 - 9 5.29 16.5 
4845 5.77 1.8 4.2 2 5.45 18 
4846 6.23 1 7.2 10 0.05 17.5 
4847 10.77 1.2 7 A 0.27 19 
4848 12 2.1 8.2 10 9.93 19 
4849 - 1.7 7.7 5 13.2 17.5 
4850 10.47 2 8.2 10 10.53 17.5 
4851 8.53 2.6 5.1 = 3.27 17.5 
4852 8.43 1.7 5.9 10 5.45 17.5 
4854 7.93 4.3 8.2 10 8.18 18 
4855 9.47 1 7.8 5 10.47 17 
4856 8.33 4 8.2 2 7.64 17.5 
4857 9.33 4.4 8.2 10 13.25 18 
4858 10.03 1.7 7 10 13.04 19 
4859 9.6 4 8.2 10 8.35 18 
4860 7 oe - 7 13.8 12 
4861 8.77 1.6 6.9 4 8.24 19 
4862 7.87 1.6 7.6 10 11.29 15.5 
4863 7.07 2.9 8.2 10 15.98 18 
4865 8.77 2.4 5.1 9 5.89 17 
4866 8.1 1.8 4.2 10 4.36 18 
4867 7.27 2.7 5.2 10 2.18 18 
4868 7.13 1.6 7.3 10 0.65 17.5 
4870 8.77 4 7.7 10 2.29 18 
4871 10.37 1.9 4.7 10 13.53 18 
4872 9.73 4 8.2 10 4.36 18 


Table 14 (continued) All results compiled into a single table, Business 
Mathematics, 2005-6. 


The total score has the mean 49.04, median 49.53, minimum 25.8, maximum 
70.81 and standard deviation 9.31. The score and ranking for every student 
are shown in Table 15. 
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ID 
4661 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
A811 
4813 
4814 
4815 
4816 
4817 
4818 


Table 15 Total score and ranking, Business Mathematics, 2005-6. 


Score 
38.77 
50.77 
53.06 
53.72 
52.8 

52.16 
51.42 
35.55 
34.1 

54.09 
48.51 
42.25 
56.08 
56.49 
47.48 
34.46 
61.55 


ID 
4819 
4820 
4821 
4822 
4823 
4824 
4825 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 
4836 
4837 


Score 
42.64 
28.52 
32.4 

53.96 
52.82 
63.42 
36.58 
58.1 

54.68 
55.28 
37.06 
48.82 
42.27 
46.23 
43.73 
48.19 
47.2 


ID 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 
4848 
4849 
4850 
4851 
4852 
4854 
4855 


Score 
41.22 
50.25 
58.47 
62.24 
53.95 
70.81 
37.42 
37.22 
41.99 
42.24 
61.23 
45.1 

58.69 
37.01 
48.59 
56.22 
50.74 


Score 
47.67 
63.19 
60.77 
58.15 
25.8 

48.5 

53.86 
62.65 
48.16 
46.46 
45.35 
44.19 
50.76 
58.49 
54.3 
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ID 
4661 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
A811 
4813 
4814 
4815 
4816 
4817 
4818 


Score 
38.77 
50.77 
53.06 
53.72 
52.80 
52.16 
51.42 
35.55 
34.10 
54.09 
48.51 
42.25 
56.08 
56.49 
52.75 
34.46 
61.55 


ID 
4819 
4820 
4821 
4822 
4823 
4824 
4825 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 
4836 
4837 


ID 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 
4848 
4849 
4850 
4851 
4852 
4854 
4855 


Score 
41.22 
50.25 
58.47 
62.24 
53.95 
70.81 
41.58 
37.22 
41.99 
42.24 
61.23 
56.38 
58.69 
41.12 
48.59 
56.22 
50.74 


Table 16 Score and rank after adjustment for absence from tests. 


Score 
47.67 
63.19 
60.77 
58.15 
51.60 
48.50 
53.86 
62.65 
48.16 
46.46 
45.35 
44.19 
50.76 
58.49 
54.30 
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Total exam score distribution, Business Maths, 2005-6 
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Figure 22 Distribution of the total score of tests, attendance and home- 
work. 


From the boundaries of each group portrayed in Figure 22 we arrive at the 
grading scheme in Table 17. 


Range Grade 
66.81, 100] A 
59.95, 66.81) Br 
49.66, 59.95) B 
39.38, 49.66) cr 
31.38, 39.38) C 
28, 31.38) Dr 
0, 28) (N/A) 


Table 17 Grading scheme derived from Figure 22. 


And then the grades, together with total scores, are shown in Table 18. 
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ID 
4661 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
4811 
4813 
4814 
4815 
4816 
4817 
4818 
4819 
4820 
4821 
4822 
4823 


Table 18 Total score and grade, Business Mathematics, 2005-6. 


Total 
43.07 
50.77 
53.06 
53.72 
52.80 
52.16 
51.42 
35.55 
34.10 
54.09 
48.51 
42.25 
56.08 
56.49 
52.75 


Grade 
(ens 


n 
T 


| 
T 


n 
T 
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ID 
4824 
4825 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 
4836 
4837 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 


Total 
63.42 
40.65 
58.10 
54.68 
55.28 
37.06 
48.82 
42.27 
46.23 
43.73 
48.19 
47.20 
41.22 
50.25 
58.47 
62.24 
53.95 
70.81 
41.58 
37.22 
41.99 
42.24 


Grade 
Bt 
ct 


ct 
ct 


ID 
4848 
4849 
4850 
4851 
4852 
4854 
4855 
4856 
4857 
4858 
4859 
4860 
4861 
4862 
4863 
4865 
4866 
4867 
4868 
4870 
4871 
4872 


Total 
61.23 
56.38 
58.69 
41.12 
48.59 
56.22 
50.74 
47.67 
63.19 
60.77 
58.15 
51.60 
48.50 
53.86 
62.65 
48.16 
46.46 
45.35 
44,19 
50.76 
58.49 
54.30 


Grade 
Bt 
B 
B 
ct 
ct 


Here student 4860 who had been absent from too many exams was given a 
C* instead of a B. Also, she said she would come to Quiz 3, but did not. 
That makes her marks fairly 51.60(0.9) = 46.44, that is a CT. 
Student 4661 had not come to Quiz 1. According to our procedure, therefore, 


the total of his score becomes 38.77/0.9 = 43.07, and a C* at that. 
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Judge not, that ye be not judged. 
Matthew 7:1 
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Opposite is the back cover of the book Percolation within percolation and 
Voronoi Tessellation by Kit Tiyapan, 2003. The pictures are tessellated with 
respect to the golden ratio. 


The smallest picture is the creek just below Kinder Downfall. The next- 
smallest one was taken on top of the plateau at Kinder Scout, at its rim where 
the cliffs drop down until they meet Kinder Reservoir, and the smallest one 
larger than it one of the locks on a canal we pass on our way to Hayfield. 


Kinder is a german word that means ‘children’. So the name of the place 
should more correctly be written ‘Kinderscout’. But then again there are 
Kinder Downfall and Kinder Low. To do the same thing everywhere would 
probably result in something seemingly out of place sitting in an English 
context. Therefore the name is normally written ‘Kinder Scout’. 


Kit Tyabandha 
Bangkok, April 2006 
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